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Preface 


For  those  who  have  had  the  opportunity  to  view  real 
staring  infrared  data  may  have  observed  instances  when  a 
target  blob  (aircraft)  moving  across  the  focal  plane  seemed 
to  possess  "target-like"  qualities.  These  glimpses  of  target 
characteristics  may  not  be  obvious  in  single  frames  of  data, 
but  can  be  inferred  from  the  motion  of  the  target.  The 
purpose  of  this  study  is  to  develop  processing  methods  that 
exploit  the  dynamic  nature  of  the  target  by  recovering  enough 
information  to  identify  the  target.  Unfortunately,  actual 
staring  sensor  data  was  not  available  to  test  the  algorithms 
developed  here. 
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Abstract 

A  unique  method  to  extract  and  classify  aircraft  targets 
from  low-resolution  staring  infrared  data  is  presented  which 
uses  sequences  of  sensor  frames  to  obtain  target  signature 
samples.  The  classifier  is  based  on  a  target  feature 
description  consisting  of  X  and  Y  energy  projections  of  the 
target.  The  target  information  is  first  extracted  from 
sensor  data  through  sampling  windows  consisting  of  contiguous 
rows  and  columns  of  detectors  in  the  path  of  the  expected 
target.  The  samples  are  then  reduced  to  remove  degradations 
from  noise  and  sensor  optics.  The  classifier  assumes  that 
approximate  target  position  and  velocity  are  available  from  a 
target  detection  and  track  subsystem.  Tests  were  performed 
using  a  software  sensor  simulation  that  includes  optical  blur 
and  sensor  motion  models.  The  classifier  is  shown  through 
simulation  to  successfully  extract  and  classify  various 
transport  and  fighter  class  aircraft  targets  with  sensor 
footprint  resolution  of  11  to  45  meters. 


x 


the  in-phase  case  and  given  as 


=  L  (11) 

where : 

=  the  resolution  of  the  in-phase  target  samples 
in  meters 

The  resolution  of  the  target  in  the  previous  example  crossing 
a  9  meter  footprint  would  then  be  9/4  or  2.25  meters. 
However,  in  the  more  likely  case,  if  the  target  is  "out  of 
phase"  with  the  sensor,  further  resolution  is  possible.  The 
amount  of  increased  resolution  is  fixed  for  a  rational  Fp 
and  is  given  as 


where : 


R 


o 


REM  ( ) 


(12) 


the  resolution  of  the  out-of-phase  target 
samples  in  meters 

the  fractional  part  of  the  enclosed  expression 


and  the  number  of  frames  Fn  required  to  resolve  the  target  is 
given  as: 


REM  (L/Rj_) 


(13) 


For  example,  consider  a  one-dimensiona  1  case  with  a  sensor 
having  a  footprint  of  20  meters,  and  a  frame  rate  of  40  Hz  in 
which  a  target  with  a  velocity  of  250  meters/sec  passes 
through  a  single  DFOV.  The  sample  width  is  then  250/40  or 
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to  the  detector  boundary: 


td  =  _k_  (9) 

vt 

where : 

=  target  dwell  time 
Vt  =  is  target  velocity 

When  the  target  dwell  time  is  related  to  the  system  frame 
rate,  an  important  parameter  which  describes  the  number  of 
times  a  point  on  the  target  will  be  sampled  as  it  moves 
across  the  detector  footprint  is  generated.  This  parameter 
is  designated  Fp  and  is  given  by, 

Fp  =  Fr  td  (10) 

where : 

F_  =  the  number  of  frames  during  the  dwell  time  of 
the  target 

Ff  =  the  system  frame  rate  (Hz) 

Thus,  a  target  which  has  a  dwell  time  of  0.1  seconds  passing 
through  an  DFOV  with  a  frame  rate  of  40  Hz  can  expect  to  have 
40*0.1  or  4  frame  samples  per  detector. 

The  ultimate  resolution  of  the  acquired  data  is  not 
fixed  by  the  sample  rate,  but  by  the  relative  phasing  of  the 
target  velocity  vector  to  to  the  DFOV.  For  the  previous 
example,  the  target  may  be  considered  "in  phase"  when  the 
target  crosses  the  footprint  in  a  integer  number  of  samples 
(Fp  is  an  integer).  Thus,  the  resolution  is  then  fixed  for 
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Figure  2-3.  Example  of  Simulated  Optical  Blur  Function. 


The  combined  effect  of  equations  (6)  and  (7)  on  equation 
(2)  gives  the  following  expression  describing  the  spatial 
transformations  to  the  incident  energy  on  the  detector, 


si,j=  JJ  pDFOVHd  (f  x'f  Ho  (f  x'f  y>exP  t  2j  (f  xx  +  f  yY>  ]  xdfy  (8) 
where : 

— i ,  j  =  the  output  of  detector  at  coordinates  (i,j) 

PqPOV  =  the  spectrum  of  PqFOV 

The  temporal  target  response  of  a  mosaic  sensor  depends 
on  the  target  speed  and  the  crossing  geometry.  The  time  in 
which  the  target  takes  to  cross  ^-he  footprint  of  a  single 
detector  is  referred  to  as  the  target  dwell  time;  the  dwell 
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Fourier  transform  is  given  by  (Ref  4), 


0* 


G(fX/fy)  =  JJ  g(x,y)  exp[-j2„  (fxx+fyy)  dxdy  (4) 
and  the  corresponding  inverse, 

g(x,y)  =  JJ  G(fx,fy)  exp [  j2*  (fxx+fyy)  dfxdfy  (5) 

The  detector  is  modeled  as  a  square  aperture  (of  area  =  L2 ) 
with  uniform  response  whose  transfer  function  is  given  by, 

Hd(fx,fy)=  sin (  Lfy  )  sin(  Lfy  )  (6) 

L^x  Lfy 

and  the  optical  transfer  function  which  is  approximated 
by  a  two-dimensional  gaussian  is  given  by, 

Hc(fx,fy)  =  exp  [-2  fx2  *2]exp[-2  fy2  „  2]  (7) 

where : 

°  -  the  approximate  width  of  a  detector  for  a 

matched  optic  configuration 

Figure  2-3  shows  the  gaussian  point-spread  function  (PSF) 
used  in  the  computer  simulation  along  with  the  relative  shape 
of  the  point-spread  function  compared  to  the  detector.  This 
function  is  convolved  with  input  scenes  to  simulate  the 
blurring  of  diffraction-limited  optics.  In  effect,  the  PSF 
acts  as  a  low -pass  filter  attenuating  high  spatial 
frequencies  present  in  input  scenes.  Removal  of  the  PSF  will 
play  a  significant  part  in  the  recovery  of  the  target  signal 
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of  data,  with  the  tiire  between  samples  referred  to  as  the 
frame  time.  The  mapping  of  energy  from  the  scene  onto  the 
focal  plane  can  be  model  in  terms  of  the  scene  radiant 
intensity  le-  The  scene  is  modeled  in  2-D  coordinates 
(u,v),  where  the  energy  transfer  to  the  sensor  fron,  point 
(u,v)  is 


P ( u , v )  =  Ie(u,v)Tan 


(2) 


where : 

P(u,v)  =  power  in  watts  at  the  sensor  from  point  (u,v) 

I  (u,v)=  radiant  intensity  in  watts/ steradian  from  a 
point  (u,v)  in  object  space 
n  =  represents  the  solid  angle  where  n  =  A0/R^ 

Ta  =  is  the  transmittance  of  the  atmosphere  which 

will  not  be  modeled  in  subsequent  calculations 

It  is  convenient  to  express  the  two-diniensional  nature  of  the 
energy  transfer  in  terms  of  a  number  of  point  intensities 
since  the  simulation  will  model  each  point  in  the  image  as  a 
integration  of  a  infinite  number  of  points  in  a  scene. 

PDFOV  =  If  P(u,v)  dxdy  (3) 

DFOV 

The  transfer  functions  of  both  the  optics  and  detector 
are  assumed  to  be  linear  and  shift  invariant  such  that  the 
detector  signal  can  be  modeled  in  the  Fourier  domain  as  a 
multiplication  of  the  2-D  Fourier  transform  of  the  signal 
with  the  transfer  function  of  the  optics  and  the  transfer 
function  of  the  detector  aperture.  The  two-dimensional 


Energy  radiated  from  the  target  and  the  background  is  imaged 
onto  a  mosaic  array  of  detectors  referred  to  as  a  focal 
plane.  Each  detector  receives  energy  from  an  area  subtended 
by  the  detector  f ie ld-of-view  (DFOV)  as  shown  in  Figure  2-2. 
The  spatial  energy  distribution  on  the  surface  of  the 
detector  can  be  considered  a  convolution  with  the  transfer 
functions  of  the  sensor  system  and  the  incident  target  and 
background  energy.  Similar  models  have  been  described  in 
(Ref  1,15,18),  however  the  target  under  investigation  is 
usually  represented  as  a  point  source  or  a  gaussian 
distribution.  Since  the  resolution  of  the  sensors  supported 
by  this  effort  range  from  10  to  100  meters,  it  is  important 
to  model  the  target  as  a  two-dimensional  signature  having 
specific  shape  characteristics  which  will  be  described  later. 

Several  sensor  specifications  have  a  direct  influence  on 
target  acquisition.  These  include  detector  geometry,  optics, 
range,  and  noise.  The  extent  of  the  DFOV  is  usually 
referred  to  by  the  area  in  the  scene  subtended  known  as  the 
detector  "footprint".  Thus,  the  footprint  is  represented  by 

L  =  R  sin  (©dfov^ 

where : 

L  =  width  of  footprint  in  meters  (area  =  ) 

F  =  is  the  range  (altitude) 

During  observation,  the  focal  plane  is  sampled  at  a  specified 
rate  to  generate  data.  Each  sample  is  referred  to  as  a  frame 
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The  system  effects  can  be  further  separated  into  sensor 
geometry  and  motion,  optical  and  detector  aperture  functions 
and  transmission,  and  background  noise  both  direct  and 
induced  by  sensor  instabilities.  Models  describing  sensor 
performance  and  associated  noise  characteristics  have  been  a 
topic  of  research  for  several  years  (Ref  1,7,8,12,16,18). 
The  intent  of  this  section  is  to  review  some  of  these  models 
as  they  apply  to  this  specific  problem,  and  to  support 
simulation  of  these  models  in  later  sections. 


Sensor  Scenario 

The  scenario  under  study  consists  of  a  high-altitude 
down-looking  staring  mosaic  sensor  and  a  target  in  cruise 
flight  passing  through  the  sensor  f ie ld-of-view  (SFOV). 


SENSOR 


TARGET 


BACKGROUND  — J  Fooifrmt 


Infrared 
Focal  Plano 

Syales  Transfer 
Function 


Detector 
Field -of -View 


TARGET  CLASSIFICATION  SCENARIO 


Figure  2-2.  Staring  Sensor  Scenario 


1 1  .  GENERAL  THEORY 


Introduction 

To  properly  analyze  the  signal  generated  from  a  staring 
mosaic  sensor,  it  is  necessary  to  model  the  essential 
character istics  of  the  sensor,  background,  and  target.  This 
analysis  of  the  interaction  between  the  target  signal  and 
various  degradation  effects  inherent  in  the  acquisition 
process  provide  a  basis  for  developing  an  approach  to 
counteract  these  effects.  Use  of  a  staring  infrared  sensor 
for  target  classification  presents  a  number  of  problems  which 
impact  the  quality  of  the  acquired  signal.  The  signal 
ultimately  used  for  target  correlation  must  pass  through 
several  stages  of  acquisition  each  of  which  contribute 
*  various  levels  of  degradation  to  the  original  signal.  These 

noise  contributions  are  generally  separated  into  those  caused 
by  the  scene  which  the  sensor  is  viewing,  and  the  sensor 
itself.  An  overview  of  this  process  is  shown  below  in  Figure 
2-1. 


Figure2-1.  Staring  Sensor  Data  Acquisition. 
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derived  and  implemented  to  obtain  a  performance 
characterization  relative  to  a  limited  number  of  target 
types.  Finally,  a  simulation  modeling  critical  aspects  of 
staring  sensor  signal  acquisition  process  is  developed  to 
generate  data  to  test  the  aforementioned  target  extraction 
scheme . 

Presentation  of  this  research  will  proceed  with  a 
discussion  of  general  theory  concerning  the  scenario 
dependent  sensor,  target,  and  background  considerations.  This 
is  followed  by  a  discussion  of  a  target  feature  set  derived 
from  image  sequences  and  later  applied  to  a  feature 
extraction  algorithm.  This  algorithm  is  then  implemented  and 
tested  in  a  computer  simulation  to  determine  the  feasibility 
of  low-resolution  target  classification.  Lastly,  the  results 
and  analysis  of  algorithm  performance  are  presented  with 
recommendations  for  continuing  study. 

Materials  and  Equipment 

All  thesis  research  was  performed  in  the  Signal 
Processing  Laboratory  at  AFIT  utilizing  the  Data  General 
Eclipse  S/250  and  Nova  computers  and  the  Octek  Image 
Analyzer  which  can  acquire  and  display  a  320x240  monochrome 
image  display  with  4  bit  resolution.  Target  and  background 
data  consisting  of  photographs  and  models  are  input  to  the 
Octek  digitizer  by  a  video  camera.  Programiming  was 
accomplished  using  FORTRAN  IV  and  FORTRAN  V  when  possible. 
Existing  support  software  is  utilized  as  available. 
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However,  it  is  expected  that  a  certain  amount  of  training  is 
required  to  optimize  any  feature  set  correlation  scheme  and 
that  a  practical  implementation  of  these  techniques  would 
utilize  actual  target  signatures. 

Specification  of  sensor  noise  contributions  will  be 
important  in  evaluating  the  classification  potential  of 
target  acquisition  techniques.  The  sensor  itself 
contributes  a  certain  amount  of  noise  to  the  image  in  the 
form  of  an  optical  transfer  and  detector  aperture  function 
(blurring),  jitter  and  drift  (movement  of  the  sensor),  and 
detector  variations.  It  is  assumed  that  contributions  due 
the  detector  electronics  will  be  background  limited  and 
therefore  neglected.  Likewise,  temporal  variations  due  to 
scintillation  and  other  sources  will  not  be  included. 
However,  contributions  due  to  spatial  gradients  in  the 
background  scene  are  included  to  a  limited  extent. 

Approach 

To  understand  the  impact  of  various  degradation 
mechanisms  inherent  in  the  staring  sensor,  the  target, 
background,  and  sensor  models  are  represented  initially  in  a 
simplified  manner  to  facilitate  identification  of  a  useful 
target  feature  set.  The  extraction  of  a  feature  set  from  the 
limited  amount  of  information  generated  by  a  target  is 
important  to  the  success  of  the  classification  scheme  and  is 
studied  extensively.  A  feature  extraction  method  is  then 


Algorithms  for  target  detection  are  not  included  due  to 
thesis  time  constraints.  Therefore,  any  implementation  of 
these  techniques  will  require  hand-off  information  from  a 
target  detection  scheme  concerning  target  position  and 
relative  direction. 

Assumptions 

Since  the  distance  between  the  sensor  and  the  target  is 
assumed  to  be  relatively  large  (on  the  order  of  50-100 
km),  the  target  may  be  represented  by  a  two- 
dimensional  projection.  An  aircraft  for  example,  would 
appear  as  a  flat  2-D  form  as  observed  from  space. 
Therefore,  all  target  classification  techniques  considered 
in  this  study  will  utilize  a  2-D  projection  of  the 
target  on  an  earth  background.  This  eliminates  the  need 
for  3-D  geometric  transformations  of  objects  thereby 
reducing  the  recognition  problem  to  a  single  aspect  angle 
of  the  target.  Given  this  simplification,  the  major 
considerations  for  a  recognition  methodology  will  be 
target  orientation  and  size  invariance. 

It  is  expected  that  the  target  signature  will  be 
a  known  quantity  such  that  features  such  as  size  and  shape 
may  be  utilized  to  identify  the  target.  Target  signatures 
used  in  this  study  are  intended  to  represent  a  first-order 
approximation  to  actual  signatures  which  are  unavailable  due 
to  security  classification.  It  is  acknowledged  that  this 
representation  could  affect  the  performance  analysis. 
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Scope 

This  effort  is  primarily  concerned  with  the  feasibility 
of  extracting  useful  target  information  from  low-resolution 
image  sequences  and  the  investigation  of  various  sensor  and 
scenario  dependent  parameters  which  impact  target 
classif ication.  Since  the  purpose  of  this  research  is  to 
support  the  Foreign  Technology  Division’s  (FTD)  need  to 
analyze  specific  sensor  data,  scenarios  used  in  this  study 
consist  of  a  high  altitude  IR  staring  sensor  with  a  nadir 
view  of  an  earth  background.  Other  geometries  are  not 
considered.  Scenes  are  generated  by  computer  simulation 
which  model  critical  effects  of  background  clutter,  target 
motion  and  sensor  optics.  Background  scenes  consist  of 
i Wb  visible-band  photographs  which  are  intended  to  provide  a 

first-order  approximation  of  an  infrared  cluttered 
background.  The  primary  objectives  of  this  effort  are: 

1.  Demonstrate  feasibility  of  classifying  various 
aircraft  types  from  low-resolution  data. 

2.  Investigate  various  scenario  dependent  parameters 
which  impact  the  target  recognition  process. 

3.  Perform  a  performance  analysis  to  ascertain  the 
potential  application  of  techniques  developed  herein  to  a 
practical  system  application. 

4  Develop  generalized  algorithms  which  may  be  applied 
to  specific  sensor  data  as  required  by  FTD. 
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Thus,  the  target  appears  as  a  single  dot  or  small  blob  and 
is  usually  characterized  and  identified  by  its  amplitude, 
direction,  and  sometimes  velocity.  However,  target  motion 
across  the  sensor  f ield-of-view  can  provide  enough  additional 
information  to  allow  target  classification.  Although  using 
image  sequences  to  enhance  resolution  is  not  new,  this 
application  provides  a  unique  approach  to  the  target 
classification  problem  and  may  serve  to  reduce  false  alarm 
rates  plaguing  current  system  performance. 

Problem 

The  problem  to  be  addressed  by  this  effort  is  the 
development  of  a  methodology  to  classify  low-resolution 
targets  based  on  features  extracted  from  temporal  variations 
in  sequences  of  infrared  imagery.  This  problem  can  be 
divided  into  four  parts: 

1.  Detection  and  location  of  potential  targets  in  a 
noisy  scene. 

2.  Extraction  of  target  information  by  observing 
target  translation  through  sequential  frames  of  data. 

3.  Transformation  of  acquired  data  into  a  useful  set 
of  features  unique  to  the  target. 

4.  Correlation  of  derived  features  with  those  of  known 


tarnets . 


LOW  RESOLUTION  TARGET  CLASSIFICATION 
FROM  A  STARING  INFRARED  SENSOR 


I .  Introduction 

Background 

A  high-altitude  or  space-based  infrared  sensor  using 
a  staring  mode  of  operation  offers  increased  resolution 
and  sensitivity  over  current  scanning  systems  for 
surveillance  missions.  A  staring  sensor  employs  millions 
of  detectors  compared  to  several  thousand  for  the  scanning 
system  and  can  look  or  "stare"  at  one  point  for 
relatively  long  periods  of  time,  thus,  collecting  more 
incident  energy  during  observation.  The  technology  for 
developing  the  detectors  and  associated  processing  for 
staring  systems  is  rapidly  maturing  and  prototype 
systems  are  scheduled  for  test  over  the  next  several 
years.  The  Foreign  Technology  Division  (FTD)  at  Wright- 
Patterson  AFB  has  a  requirement  to  analyze  data  provided 
by  these  prototype  staring  systems  and  to  develop  the 
associated  processing  techniques  which  will  automatically 
locate  and  identify  1 ow-reso  1  u t i on  targets  of  interest. 

Although  many  processing  techniques  have  been 
presented  in  the  literature  which  can  detect  a  target 
in  noisy  backgrounds,  these  techniques  consider  targets  of 
interest  to  be  "point"  targets.  Therefore,  information 
pertaining  to  the  target  form  or  type  is  largely  ignored. 
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6.25  meters  and  20/6.25  or  approximately  3  frames  would  be 
required  to  sample  the  data.  However,  if  more  samples  of 
data  were  used,  the  resolution  could  be  increased  to  RQ  = 
6.25  REM  (20/6.25)  or  1.25  meters  with  5  frames  of  data 
(illustrated  in  Figure  2-4.)  This  is  accomplished  by  summing 
the  samples  from  subsequent  detectors  with  that  of  the  first 
detector.  The  result  is  an  equally  spaced  sampled  energy 
profile  of  much  higher  resolution  then  samples  from  a  single 
detector.  Theoretically,  for  an  irrational  Fp,  the 
potential  resolution  would  be  infinite  for  an  infinite  number 
of  frames,  i.e.  lim  R_  =  0.  Practically  of  course  this  would 

F  '♦» 

N 

be  impossible  to  implement.  Therefore,  the  irrational  Fp  is 
aproximated  by  a  rational  value  with  little  error  induced  due 
to  the  low  frequency  energy  curve  which  the  data 
approximates.  The  application  of  this  resolution  enhancement 
is  a  key  factor  in  determining  a  useful  feature  extraction 
method  as  will  be  shown  in  the  next  chapter.  The  practical 
application  is  limited  however  due  to  sensor  noise  and  other 
signal  degradation  phenomena  depending  on  the  number  of 
detector  samples  to  be  reduced  to  a  single  sample. 
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SPECIFICATIONS: 


SENSOR  FOOTPRINT  -  20  METERS 
SENSOR  FRAME  RATE  -  40  Hz 
TARGET  VELOCITY  -  250  METERS/SEC 
SAMPLE  RESOLUTION  -  250/40  =6.25  METERS 


0  20  40  60  80  100  meters 


BY  COMBINING  THE  FIVE  DETECTOR  SAMPLES,  6. 25 -REM (20/6. 25)  OR 
1.25  METER  SAMPLE  RESOLUTION  CAN  BE  OBTAINED: 


DETECTOR  1 
♦ 

DETECTOR  2 
+ 

DETECTOR  3 
♦ 

DETECTOR  4 
♦ 

DETECTOR  5 


The  two-dimensional  scene  observed  from  the  staring 
mosaic  space  platform  can  contain  many  features  that  tend  to 
obscure  or  degrade  the  sensor's  target  detection  capability. 
When  the  target  is  weak  relative  to  the  background,  (such  as 
with  an  aircraft)  the  majority  of  the  incoming  signal  is 
considered  to  be  noise  (Ref  15:103).  Although  the  potential 
to  detect  weak  targets  against  structured  background  is  a 
primary  advantage  for  staring  IR  sensors,  the  background  can 
not  be  eliminated  completely  due  to  sensor  drift,  jitter, 
clutter,  cloud  movement  and  electronic  noise. 

Small  instabilities  in  the  sensor  can  lead  to 
background  spatial  variations  being  interpreted  as  temporal 
change  resulting  in  an  AC  signal  from  the  detector. 
Therefore,  the  earth  background  is  considered  a  noise  source. 
Generally,  signal  processing  must  be  performed  to  raise  the 
level  of  the  signal  relative  to  the  level  of  the  noise. 
The  type  of  noise  or  background  suppression  technique  used 
is  dependent  upon  the  nature  of  the  sensor  noise  sources. 

There  are  four  major  sources  of  noise  in  a  staring 
mosaic  sensor.  These  include  photon  fluctuations  of  incoming 
radiation  from  both  targets  and  background,  temporal 
variations  in  background  radiance,  spatial  variations  due  to 
movement  between  the  sensor  f ie 1 d-of -view  (SFOV)  and  the 
background,  and  finally,  noise  from  the  detector  and 
associated  electronics.  These  noise  sources  are  discussed  in 


detail  in  Pohlman  (1976)  and  Pollock  (1981).  Typically, 
sensors  are  designed  to  minimize  noise  sources  while 
enhancing  the  target  signal.  However,  it  is  extremely 
difficult  to  design  against  all  background  situations. 

A  prominent  noise  source  which  every  detection  system 
must  deal  with  is  clutter.  According  to  Pollock 
(1980),  "Clutter  may  be  defined  as  any  signal  which  must 
be  processed  to  some  degree  to  determine  if  it  is  a 
target  or  a  non-target”.  Sources  of  clutter  depend  upon  the 
sensor  and  viewing  conditions,  but  may  range  from  clouds  and 
glint  from  ocean  specular  reflection,  to  structured  patterns 
of  corn  fields  and  industrial  heat  sources  such  as  smoke 
stacks.  Schemes  for  suppressing  background  clutter  include 
adaptive  thresholding  based  upon  the  statistical  nature  of 
the  background  and  temporal  filtering  to  subtract  slow 
moving  backgrounds  from  fast  moving  targets.  A  technique  for 
suppressing  clutter  due  to  motion  of  the  sensor  (jitter)  is 
presented  in  Mazaika  (1982)  in  which  he  develops  a  model  for 
clutter  leakage  due  to  gradients  in  the  background  radiance. 
Jitter  is  a  significant  effect  in  many  sensor  designs  and  may 
severely  limit  the  detection  capability  of  the  sensor. 
Measurement  of  the  influence  of  temporal  and  spatial 
variations  is  usually  accomplished  by  spectral  analysis  of 
the  background  signal. 

Another  important  effect  of  the  background  is  its 


relative  temperature  contrast  with  the  target. 


Target 


at 


information  integrated  over  the  detector  footprint  can  be 
lost  due  to  the  masking  of  background  radiance  with  that  of 
the  target  signature.  This  process  is  not  reversible.  To 
see  this,  consider  the  2-D  signal  imaged  onto  the  focal  plane 
without  sensor  transfer  effects.  The  signal  on  the  detector 
at  tinne  t  is  given  by, 

sij(t)  =  At ( t ) Et  +  [L2  -  AM(t) ]EBKG(i,  j)  (14) 

where 

Sij  (t)  =  the  signal  on  detector  (i,j)  of  the  focal  plane 
ATlt)  =  the  area  of  the  target  in  the  DFOV  at  time  t 
AM(t)  =  the  area  of  the  radiant  background  masked  by 
the  target  (varies  from  0  to  AT(t))  at  t 
E-j.  =  the  average  irradiance  of  the  target 

in  watts/meter2 

EBKG  =  average  irradiance  from  the  background 

in  watts/meter2  on  detector  (i,j) 

The  effect  is  illustrated  in  Figure  2-4.  Given  a  3  by 
3  array  of  detectors,  a  single  cloud  partially  obscures  the 
field  of  view  of  the  central  detector  (all  other  background 
is  negligible).  For  illustration,  the  cloud  has  a  uniform 
signal  of  level  5  watts/meter.  Using  Equation  (14),  the 
signal  due  to  the  background  cloud  is  5L2.  If  the  target 
(also  with  a  uniform  signature  of  5  watts/ meter)  passes 
through  the  center  detector  along  path  1,  the  target  signal 
will  appear  on  the  detector  rising  to  a  maximum  amplitude  of 
5At  +  5L2.  However,  if  the  target  takes  path  2,  the  relative 
contrast  is  zero  since  5AT+512-5AT  =  5L2.  Therefore,  no 
target  information  will  be  present  on  the  center  detector 
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using  path  2.  This  clutter  severely  degrades  effect  to  the 
target  signature  in  low  signal  to  noise  ratio  scenarios. 


Figure  2-5.  Signature  Degradation  Due  to  Background  Contrast 


Target  Model 

The  model  for  the  target  signature  should  be 
representative  of  expected  target  signatures  in  the  8-12 
micron  band.  The  target  signature  for  this  region  is  largely 
due  to  the  target  thermal  temperature  and  resulting  blackbody 
emittance.  From  a  high-altitude  sensor  however,  the  actual 
signature  depends  on  the  radiance  contrast  between  the  target 
and  its  immediate  background.  Therefore,  the  spatial  extent 
of  the  target  can  be  represented  as  a  2-D  silhouette  whose 
shape  is  invariant.  Similar  approaches  to  representing 
thermal  targets  are  presented  in  several  studies  analyzing 
stationary  tank  and  ship  targets  (Ref  2,5,6,11,19).  For  this 
analysis,  the  target  motion  is  restricted  to  2-D  translation 
and  rotation  in  a  plane  parallel  to  the  2-D  image  plane.  For 
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the  purposes  of  classification,  it  is  assumed  that  the  target 
signature  is  a  known  quantity.  No  attempt  has  been  made  to 
study  variations  in  the  target  signature  and  it  is  assumed 
that  models  are  available  (or  can  be  developed)  for  a 
practical  implementation.  For  a  first  order  approximation, 
the  target  signature  results  from  the  contrast  between  the 
target  and  background.  This  can  include  both  positive  and 
negative  temperature  contrasts  which  can  vary  as  the  target 
moves  across  a  highly  structured  background.  The  target  is 
modeled  here  as  a  relatively  uniform  intensity  distribution 
over  the  target  shape  tapering  at  the  edges.  A  more  detailed 
model  would  account  for  the  the  difference  in  thermal 
signature  due  to  the  engine  temperature  and  the  adiabatic 

heating  of  the  aircraft  skin.  For  this  model,  the  engine  is 

represented  as  a  graybody  whose  radiance  is  proportional  to  a 
blackbody  at  engine  temperature  Te.  The  temperature  of  the 
aircraft  skin  Ts  is  given  by: 

Ts  =  t0[1  +  r (  7-  1)  M2]  (17) 

where : 

Ts:  Adiabatic  wall  temperature  in  Kelvin 

TQ:  Temperature  of  Flow 

r  :  Recovery  factor  (.85  -  1) 
y  :  Specific  heat  ratio 
M  :  Mach  number 

Both  temperatures  are  now  used  to  calculate  the  emmittance 


from  a  black  body  using  Plancks  law  Ref (9:35), 


where 


M  =  2*hc1 2  [exp(hc/AkT)  -1  ]_1  watts/m2 
“T5 


(15) 


M  =  the  blackbody  emmittance  in  watts/meter 
h  =  Plank's  constant  (6.6256x10“^^  watts/sec2) 
c  =  speed  of  light  (2.997925x10®  meters/sec) 

/  =  wavelength  (meters) 

k  =  Boltzmann's  constant  ( 1 . 38054xl0-23  watt  sec  K-1) 
T  =  Temperature  in  degrees  Kelvin  (either  Te  or  Ts) 


the  intensity  of  an  aircraft  engine  and  skin  can  be  modeled 
as  Lambertian, 


Ie  =  * M  Ae  watts/ steradian 

- 


(16) 


Thus,  the  actual  signature  of  a  target  may  be  inferred  from 
the  specific  features  of  the  target  such  as  the  number  of 
engines,  placement,  and  the  aircraft  structure  area.  A 
similar  method  to  represent  targets  is  presented  in  (Ref  2). 

The  aircraft  signatures  used  in  this  study  were  modeled 
with  very  little  detail.  Although  this  can  be  attributed  to 
the  resolution  of  the  hardware  used  in  the  simulation,  it  is 
assumed  to  be  representitive  of  a  typical  scenario.  This  can 
be  asserted  for  two  reasons: 

1.  At  long  ranges,  details  represented  by  small 
temperature  differences  are  unresolvable. 

2.  The  spatial  resolution  of  the  sensor  array 

tends  to  eliminate  high  spatial  frequencies  of  the  target. 
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The  targets  used  in  the  simulation  are  generated  by 
taking  pictures  of  plastic  scale  models  of  aircraft  using  a 
video  camera.  The  models,  which  are  painted  flat  black,  are 
placed  on  a  white  background  and  rotated  to  the  desired 
orientation  before  acquiring  the  picture.  After 
digitization,  the  graylevel  is  adjusted  as  necessary  to  match 
the  desired  signature  characteristics.  The  following  pages 
show  five  of  the  target  models  used  in  the  simulation.  Each 
target  is  shown  in  its  original  digitized  form,  and  in  its 
reduced  form  typical  of  the  targets  used  in  the  data  base  of 
the  simulation. 
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Figure  2-8.  DC-10  Target  Model. 


III.  FEATURE  SELECTION 


Introduction 

There  are  many  ways  to  represent  the  shape  and  intensity 
a  target.  As  expected,  the  features  chosen  to  represent  a 
target  have  a  direct  influence  on  the  potential  success  of 
the  classif ication  scheme.  The  ultimate  goal  of  any  feature 
selection  process  is  to  identify  those  features  which 
maximize  separability  of  the  expected  target  classes.  Other 
considerations  include  invariance  to  noise,  ease  of 
implementation,  and  computational  cost.  This  chapter  will 
focus  on  deriving  a  feature  set  which  will  use  to  the  maximum 
extent,  the  target  shape  information  extracted  from  low- 
resolution  imagery. 

Background 

Current  methods  to  represent  shape  and  intensity 
information  include  Fourier  boundary  descriptors  and 
invariant  moments  as  described  by  several  studies  (Ref 
2,5,19).  These  methods  use  features  such  as  height-to-width 
ratio,  perimeter-to-area,  normalized  intensity  moments,  and 
normalized  silhouette  moments.  However,  the  targets  used  in 
the  referenced  studies  were  of  relatively  higher  resolution 
then  those  considered  here.  An  accurate  description  of  a 
perceived  target  from  a  staring  sensor  can  be  at  best,  a 
"shaped  blob."  Thus,  the  scenes  considered  here  have  a 
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maximum  of  10  pixels  on  the  target,  with  a  more  likely 
resolution  of  3-4  pixels.  Even  prominent  target  features 
sue  hi  as  wings  and  tail  are  difficult  to  visualize  at  this 
resolution.  This  is  illustrated  below  in  Figure  10,  which 
shows  five  frames  of  a  B-52  at  high  resolution  along  side  a 
comparable  low-resolution  signature  as  would  be  observed  by 


Figure  3-1.  Target  Signature  Sequence  as 

Obtained  From  a  Staring  Sensor 


the  staring  sensor.  The  graylevel  variations  in  the  low- 
resolution  signature  are  due  to  small  changes  in  the  relative 
position  of  the  target  with  respect  to  the  focal  plane. 
While  this  resolution  may  seem  high  for  the  typical  staring 
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Redundant  Window  Sampling 

Aquisition  of  accurate  target  samples  is  obtained  in 
preprocessing  using  redundant  sampling  windows  with  limited 
size.  The  number  of  windows  required  in  each  direction  is 
arbitrary  and  depends  on  the  spatial,  temporal,  and  contrast 
of  background  clutter  in  the  local  area  subtended  by  the 
window.  For  homogeneous,  high-contrast  scenes  ,  two  to  three 
windows  may  be  sufficient.  For  backgrounds  used  in  this 
study,  five  to  ten  windows  in  each  direction  were  used.  An 
example  of  contrast  deficiencies  induced  by  high  spatially 
variant  scenes  is  shown  in  Figure  4-3.  In  the  illustration,  a 
portion  of  the  target  crossing  detector  2  of  detector  row  2 


0E7ECT0RS: 

i  /  2  /3  y  4  y  5  /  $  •' 

't 

f - 1  1 

ROM  i:  oi  [ 


SIGNAL 

FROM 

DETECTOR 

WINDOWS 

T 

ony  ?.  n  !  i 

J - '  1 — 1 

1T 

(NORMALIZED) 

.1  1 - 

ROU  3:  0 


i  r 

- 1 _ 1 - 

ROU  4:  01  ’ - 1 

Figure  4-3.  Background-Induced  Signal  Variations 

From  Parallel  Windows 

will  be  largely  lost  due  to  poor  relative  contrast.  However, 
much  of  the  same  portion  of  the  target  may  pass  through  the 
DFOV  of  detector  3,  window  3  where  greater  contrast  may  be 


Figure  4-2.  Phase  Relationship  of  Data  From  Parallel  Windows 


length  as  the  velocity  estimate  converges.  As  shown  above,  t 
is  equal  to  the  number  of  frame  samples  divided  by  the  frame 
rate  FR.  Therefore,  the  velocity  is  given  as, 


VTx  =  _k_  = 


VTy 


and  the  final  vector  as. 


1 

sin8y 


(20) 


VT  =  IVTx2  +  VTy2]l/2 

y  (21) 

Dir  =  tan^VTy  j 

In  low  noise  cases,  the  profiles  obtained  in  this  step 
may  be  adequate  for  preprocessing.  However,  since  background 
clutter  generally  corrupts  the  target  vector  profiles,  a 
single  sample  may  not  be  sufficient.  Also,  several  windows 
could  be  combined  to  obtain  an  increased  number  of  target 
samples.  Thus,  a  more  extensive  extraction  scheme  is  needed. 
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The  sampling  window  concept  differs  from  the 
mathematical  development  presented  in  Chapter  3.  As  opposed 
to  solving  a  set  of  linear  equations,  the  target  signal  is 
sampled  through  a  logical  window  as  a  convolved  sum. 
Therefore,  the  vectors  extracted  from  focal  plane  data  will 
be  deconvolved  to  solve  for  the  target  feature  vector.  This 
method  is  more  robust  than  the  simultaneous  equation  approach 
which  suffers  from  ill-conditioned  matrices  when  noise  is 
present . 

Since  hand-off  information  from  target  detection  is 
limited  to  a  relative  position  and  direction  vector  (computed 
from  target  track),  it  is  necessary  to  collect  preliminary 
information  concerning  target  velocity  and  direction  to  make 
the  extraction  process  more  efficient.  Thus,  preliminary 
windows  are  first  set-up  as  shown  in  Figure  4-1  based  on  the 
Largest  Expected  Target  (LET)  that  the  algorithm  is  expected 
to  process.  As  the  target  passes  through  the  sample  windows, 
the  profile  is  extracted  from  each  window  generating  a  phase 
relationship  between  parallel  windows  based  on  the  target 
velocity  and  crossing  geometry.  This  relationship  is 
illustrated  in  Figure  4-2,  and  allows  calculation  of  the 
target  velocity  vector  given  data  from  each  axis.  In  the 
event  of  severe  noise,  this  process  can  be  iterated  to  an 


accurate  velocity-direction  estimate  by  reducing  the  window 


’Tr 


Sampling  Windows 

Efficient  extraction  of  the  the  target  vectors  is 
accomplished  through  use  of  feature  sampling  windows.  These 
logical  windows  consist  of  columns  or  rows  of  selected 
detectors  (taking  advantage  of  focal  plane  geometry) 
continuously  sampling  the  target  profile  as  it  passes  through 
the  window  field  of  view.  In  this  way,  the  two-dimensional 
target  is  naturally  decomposed  into  orthogonal  vectors 
representing  the  target  feature  vectors  in  a  convolved  form. 
An  illustration  of  this  technique  is  shown  below  in  Figure  4- 
1,  where  two  pairs  of  windows  have  been  set-up  over  the 
length  of  the  focal  plane  (in  this  case  18  by  18  detectors) 
just  ahead  of  the  expected  target  vector. 
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Figure4-1.  Preliminary  Window  Configuration. 
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IV.  FEATURE  EXTRACTION 


Introduction 

This  chapter  presents  an  algorithm  developed  to 
accomplish  the  feature  extraction  process.  These  features 
consist  of  the  X-Y  vectors  described  in  chapter  3  which  are 
later  compared  to  a  target  database  for  classification. 
There  are  several  guidelines  closely  adhered  to  in  the  course 
of  this  development.  First  is  to  limit  the  quantity  and 
precision  of  the  required  a  priori  target  information.  It  is 
assumed  that  a  detection  scheme  exists  in  which  an  initial 
target  position  or  target  track  record  is  available. 
Therefore,  the  feature  extraction  algorithm  will  assume  this 
hand-off  is  information  available  with  some  limited  accuracy. 
Since  the  opportunity  to  observe  a  fast  moving  target  may  be 
small,  and  also  to  reduce  the  noise  effects  of  drift  and 
jitter,  it  is  desirable  to  acquire  the  target  information  in 
a  minimal  observation  time.  This  observation  time  will  vary, 
however,  depending  on  target,  sensor,  and  background  signal- 
to-noise  effects.  Also,  contributing  to  a  reduction  in  noise 
effects  is  the  sampling  of  a  smaller  portion  of  the  focal 
plane  during  target  extraction.  This  results  in  an  increase 
in  the  relative  s igna 1 - to-noi se  ratio  of  the  target  data. 
Lastly  is  the  desire  to  accomplish  the  feature  extraction 
process  with  a  minimal  amount  of  processing  resources  so  that 
implementation  of  these  algorithms  can  execute  in  real-time. 
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values.  Illustrated  in  Figure  3-5  are  two  aircraft  from  the 
previous  figure.  A  single  frame  of  target  data  is  shown 
where  the  signal  from  each  detector  is  represented  by  a 
intensity  value  in  the  Z-axis.  Below  the  single  frames  of 
data  are  the  final  target  representations  after  processing  of 
the  feature  vectors.  The  B-52  is  shown  with  a  25  meter 
footprint,  and  F-15  with  a  9  meter  footprint.  In  each 
diagram,  the  aircraft  nose  and  wing  shape  can  be  seen. 
Although  this  may  first  appear  as  a  crude  target  image,  the 
information  contained  in  these  generalized  signature 
representations  is  adequate  to  distinguish  between  various 
target  classes  over  a  range  of  scenarios. 
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range  since  targets  such  as  aircraft  are  sy  rr  met  r  i  ca  1  ale; 


one  axis.  Also,  the  memory  required  to  store  the  vec'ci  : 
(approximately  100  words  per  vector  at  2  bytes  per  point)  in 
small  so  that  a  target  could  be  fully  represented  in  les; 
than  10k  words  of  memory. 

A  three-dimensional  generalized  representation  of  the 
feature  vectors  can  be  generated  by  multiplying  the 
orthogonal  vectors  to  form  a  2-D  spatial  array  of  intensity 


Therefore,  the  ta  get  degree  characteristics  represented  by 
the  feature  vectors  are  directly  proportional  to  the 
resolution  of  the  vector.  This  also  implies  that  the 
resolution  of  the  target  vector  (or  the  number  of  points  in 
the  feature  set,)  affects  the  separability  of  the  target 
classes.  Shown  below  are  X-axis  vectors  of  three  targets  a 
B-52,  DC-10,  and  a  F-15.  Each  target  is  oriented  at  45 
degrees  and  is  sampled  at  a  one  meter  resolution. 


Figure  3-4.  Plots  of  Target  Feature  Vectors 

for  Three  Targets. 


An  obvious  disadvantage  of  this  feature  set  is  the 
variance  of  target  perspective  to  the  orthogonal  sampling 
axis.  To  fully  describe  a  potential  target,  a  target 
database  must  maintain  profiles  for  each  target  at  all 
possible  target  or ientations.  Fortunately,  most  targets  can 
be  represented  by  profiles  calculated  in  a  0  to  45  degree 


HIGH 

RESOLUTION 

SCENE 


LOU 

RESOLUTION 

SCENE 


SUM 

FRAMES 


SOLVE 

FEATURE 

VECTORS 


FEATURE 

VECTORS: 


2-D 

SCENE 

(MULTIPLY 

VECTORS) 


0 

0 

0 

0 

0 

0 

2 

2 

0 

_ 

0 

0 

0 

0 

0 

0 

4 

0 

0 

FRAME  3 


0 

0 

0 

0 

CM 

0 

0 

2 

0 

X-VECTORS  Y-VECTORS 


FRAME  1  0,0,4  0,0,4 
FRAME  2  0,2,2  0,0,4 
FRAME  3  0,4,0  0,2,2 


X1+X2+X3  =  0 

X2+X3+X4  =  0 

X3+X4+X5  =  0 

X4+X5+X6  =  0 

X5+X6+X7  =•  2 


XG+X7+X8  =  4 
X7+X8+X9  =  4 
X8+X9  =  2 
X9  =  O 


X1-X6=0,X7=2,XS=2,X9=0 


Y1+Y2+Y3  = 

Y2+Y3+Y4  = 

Y3+Y4+Y5  = 

Y4+Y5+YG  = 

Y5+YG+Y7  = 
YG+Y7+Y8  = 
Y7+Y8+Y9  = 
Y8+Y9  = 
Y9  = 

Y1-Y7=0,Y8=2,Y9=2 


X  (0,0, 0,0, 0,0, 2, 2,0) 


Y(0, 0,0, 0,0, 0,0, 2, 2) 


X9  X8 

X7 

XG  X5 

X4 

X3  X2  XI 

Y1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Y2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Y3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Y4 

0 

O 

0 

r0 

0 

0 

0 

0 

0 

Y5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

YG 

0 

0 

0 

0 

0 

0 

_0 

0 

0 

Y7 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Y8 

0 

4 

4 

0 

0 

0 

0 

0 

0 

Y9 

0 

4 

4 

0 

0 

0 

0 

0 

0 

ORIGINAL 

TARGET 


U 

Li 

u 

Li 
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detectors  in  four  frames  (the  fourth  frame  is  a  repeat  of  the 
first.)  The  detector  response  for  each  frame  is  shown  below 
the  high  resolution  frames.  By  applying  equation  18,  the 
data  are  summed  into  vectors.  These  vectors  form  general 
equations  which  are  solved  simultaneously.  The  resulting  X 
and  Y  vectors  are  the  descriptors  for  the  target.  A  two- 
dimensional  general  solution  is  obtained  by  multiplying  the 
resulting  vectors. 

The  two  vectors  X  and  Y  derived  from  the  above  possess 
many  desirable  attributes  to  classify  targets.  Both  vectors 
contain  information  on  target  size  such  as  length  and  wing¬ 
span  (within  the  resolution  outlined  in  Chapter  2).  Also 
included  are  the  shape  of  the  aircraft  structure  as  defined 
by  temperature  contrasts.  While  these  vectors  are  only 
generalized  solutions  to  specific  targets,  they  also 
represent  projections  of  the  target  form  which  will  be  shown 
to  be  separable  in  representative  signa 1 -to-noi se  cases. 
These  feature  vectors  will  be  the  basis  for  feature 
extraction  algorithms  described  in  later  chapters. 

Visualization  of  these  feature  vectors  as  energy 
profiles  of  target  signatures  reinforces  the  usefulness  of 
the  feature  set  to  describe  characteristics  of  the  target. 
The  individual  data  points  contained  in  a  feature  vector  are 
actually  an  integration  of  the  target  signature  over  the  area 
defined  by  the  resolution  width  of  the  sampling  process. 


value  of  the  subpixels  to  an  element  in  either  the  X  or  Y 
vectors.  By  solving  these  equations  simultaneously, 
solutions  for  the  subpixel  values  can  be  obtained.  This 
linear  process  reduces  the  X^  and  Yj  vectors  to  single 
vectors  X  ( 1 )  and  Y(l)  (1  is  the  length  of  the  vector  equal  to 
NFp.)  These  vectors  are  referred  to  in  subsequent 
discussions  as  "feature  vectors".  The  result  of  this 
manipulation  is  to  increase  the  resolution  of  the  extracted 
target  signature  by  combining  several  frames  of  data  into  a 
higher-resolution  one-dimensional  vector.  The  equation  for 
the  subpixel  solution  generates  either  X  or  Y  feature  vectors 
and  is  given  as: 


hki  0 

hkl  h12 
hkl  h12 


hkl  h12  h22 


0  0 


hNl  hN2 


H^l) 
I  Hi (2) 


Hx  (k) 
H2U) 
H2  (2) 

HN(k) 


where , 


Hj (k)  : 


is  the  value  of  pixel  j  of  frame  k  (either  X^ 
or  Yj) 

is  the  subpixel  i  of  pixel  j. 


A  simplified  application  of  the  previous  equations  is 
illustrated  in  Figure  3-3.  In  this  ideal  example,  the  target 
consists  of  a  square  subdivided  into  four  areas  each  with  a 
value  of  one.  The  square  moves  across  a  3  x  3  array  of 


Figure3-2. 


Summation  of  Focal  Plane  Data  to  Obtain 
Sampled  Data 

Frames  of  focal  plane  data  are  therefore  reduced  to 
vectors  and  can  be  represented  by, 

N 

Xik  =  2  Sij(k) 
j  =  l 

and,  (18) 

N 

Yjk  =  ^  sij (k) 

i  =  l 

By  acquiring  k  frames  of  data,  an  equal  number  of 
and  Yj  vectors  are  produced.  Each  element  in  the  vector 
represents  an  energy  projection  the  width  of  a  single 
detector  at  frame  k.  Since  the  target  moves  across  the 
detector  in  Fp  frame  samples,  each  of  these  samples  represent 
a  subpixel  resolution  R^  (Reference  Equation  (11)). 
Therefore,  each  vector  forms  a  general  equation  relating  the 


focal  plane.  Therefore,  tradeoffs  may  be  necessary  between 
the  sampling  time  of  the  target,  and  constraints  imposed  by 
the  system  and  scenario  dependent  noise  and  data 
requirements . 

It  would  be  desirable  to  resolve  the  target  spatially 
such  that  a  human  observer  could  recognize  the  target  by 
viewing  a  display  of  the  enhanced  target  signature.  However, 
resolving  the  target  in  a  two-dimensional  spatial  plane  is 
generally  not  possible  due  to  the  under-determined  system 
that  results  from  data  acquisistion.  To  see  this,  consider  a 
N  by  N  detector  focal  plane  sampled  at  a  rate  FR,  with  the 
number  of  samples  per  pixel  given  by  equation  (10)  as  Fp» 
Each  sample  Fp  can  then  be  thought  of  as  resolving  the  pixel 
to  Fp  square  subpixels.  To  spatially  resolve  the  target  on 
the  focal  plane,  N2Fp2  samples  would  then  be  required,  thus, 
each  detector  would  require  Fp^  samples.  However,  only  Fp 
unique  samples  are  generated.  A  modified  spatial 
representation  is  therefore  suggested  that  sacrifices  some 
spatial  features  for  the  increased  resolution  offered  by  the 
multiple  frames  of  data.  This  representation  can  be  realized 
through  summation  of  each  frame  of  data  orthogonally  into  two 
vectors  with  length  NFp,  representing  the  energy  projection 
along  the  the  focal  plane  X-Y  axis  for  that  frame  as 
illustrated  in  Figure  3-2.  In  this  example,  the  signals 
from  detectors  of  a  4x4  focal  plane  array  are  summed  during 
each  frame  of  sensor  data. 
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sensor,  one  pixel  and  subpixel  resolutions  are  also  possible 
with  large  sensor  footprints. 

As  the  previous  example  illustrates,  there  is  generally 
very  little  shape  information  apparent  in  single  frames  of 
data.  Therefore,  the  target  feature  descriptors  must 
utilize  as  much  information  as  possible  from  a  shaped-blob 
translating  across  the  SFOV.  This  information  is  generally 
in  the  form  of  a  differential  grayscale  based  on  signal 
interaction  with  the  focal  geometry.  Analysis  of  the 
interaction  between  the  target  signal  and  the  focal  plane 
geometry  is  a  key  factor  in  deriving  a  useful  feature 
descriptor . 

Feature  Derivation 

In  order  to  use  shape  information  to  classify  a  low- 
resolution  translating  target,  the  information  must  be 
"gathered"  over  several  frames  of  data.  In  this  manner, 
target  shape  information  can  be  combined  to  obtain  a 
signature  with  more  information.  The  target,  in  essence, 
will  represent  a  higher  resolution  shape  than  may  be  obtained 
from  a  single  frame  of  data.  The  amount  of  increased 
resolution  is  limited,  however,  by  the  temporal,  spatial  and 
noise  characteristics  of  a  given  scenario.  The  temporal 
limits  are  governed  by  equations  (11)  and  (12)  in  chapter  II, 
and  the  spatial  limits  by  the  geometry  of  the  detector 
footprint  and  the  resulting  column-row  configuration  of  the 
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achieved.  Exploitation  of  this  process  is  accomplished  in 
the  preprocessing  phase  of  the  feature  extraction  discussed 
in  the  following  section. 


The  position  and  length  of  the  window  boundaries  are 
based  on  the  target  velocity,  crossing  geometry,  and  the  LET. 
By  using  these  factors,  each  window  size  is  limited  to  that 
required  to  sample  the  target.  Also,  window  positioning  can 
be  set-up  to  sample  the  target  for  a  minimum  observation  time 
allowing  for  LET,  and  hand-off  error. 


sampling  window  parameters: 


Wx  =  LET  + 

sinGFp 

Wv  =  LET  + 

*  cbsBFp 

and  the  window  intercept: 


We 


We 


dx  =  LET 

2FpCos0 

dy  =  LET 

2Fpsin0 


(22) 


(23) 


A  typical  configuration  for  redundant,  limited  sized  window 
sampling  is  shown  below  in  Figure  4-5. 
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Figure  4-5.  Redundant  Window  Sampling. 


The  use  of  window  sampling  does  have  limitations. 
Targets  moving  parallel  to  either  focal  plane  axis  will  not 


be  sampled  in  that  direction  since  there  is  little  or  no 
velocity  component  in  the  perpendicular  direction.  The 
target  feature  set  is  then  reduced  in  half  limiting  potential 
cl assif icat ion.  Fortunately,  this  occurs  over  a  limited 
range  target  directions  depending  on  the  desirable  lengthof 
the  sampling  window.  For  targets  at  grazing  angles  and 
limited  backgrounds,  the  low  velocity  component  may  be  an 
advantage  in  this  case  since  longer  detector  sampling  times 
are  available. 

Vector  Processing 

Since  window  sampling  consists  of  summing  signals  from 
several  detectors  during  a  time  interval,  the  effects  of 
drift,  background  spatial  variations,  and  system  induced 
degradations  corrupt  the  original  signal  such  that 
preprocessing  is  necessary  before  classification  is  possible. 
The  preprocessing  used  to  "clean"  the  feature  vectors 
consists  of  two  major  steps;  reduction  of  the  data  to  a 
single  set  of  orthogonal  vectors,  and  deconvolution/filtering 
of  known  system  effects  such  as  the  OTF  and  detector  aperture 
functions.  The  vector  processing  proceedure  described  in  the 
following  text  is  illustrated  in  Figure  4-6. 

After  the  target  is  extracted  through  all  of  the  windows 
initially  allocated,  the  background  signal  is  then  subtracted 
from  the  window  data.  If  the  sensor  is  stationary,  this  is  a 
relatively  simple  procedure.  If  the  sensor  is  moving,  an 
estimate  of  the  drift-induced  signal  on  each  window  sample  is 
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derived.  This  is  accomplished  by  sampling  the  drift  signal 
both  before  and  after  target  passage.  If  there  is  no  drift, 
this  signal  will  be  DC.  With  sensor  movement,  the  signal 
will  be  varying  according  to  the  net  change  in  the  background 
between  samples.  For  slow  variations,  a  least  squares 
technique  may  be  used  to  estimate  the  background  signal 
during  target  passage.  However,  backgrounds  with  high  spatial 
variations  or  high  sensor  drift  may  severely  distort  the 
target  signature  so  that  any  useful  data  may  be  impossible  to 
acquire.  Unfortunately,  the  drift  signal  is  approximately 
similar  in  frequency  content  to  the  target  signal. 
Therefore,  the  target  signal  may  be  difficult  to  extract  in 
severe  drift-induced  background  noise. 

Once  the  target  signal  is  separated  from  the  background, 
the  windows  are  logically  arranged  to  obtain  the  optimum 
target  resolution.  Windows  are  grouped  into  sets  according 
to  the  phasing  of  the  target  with  respect  to  the  sensor  as 
described  in  chapter  II.  For  an  in-phase  target,  each  window 
would  be  a  separate  set  with  the  number  of  sample  windows 
equal  to  the  number  of  sets.  For  an  out-of-phase  target 
which,  for  example,  crosses  two  windows  in  an  integer  number 
of  samples,  each  set  is  then  composed  of  two  windows  with  the 
number  of  sets  equal  to  the  total  windows  divided  ty  two. 
The  information  as  to  the  arrangement  of  the  sets  is  also 
used  in  the  deconvolution  process  as  will  be  discussed  in  a 
later  section.  To  reduce  the  logical  window  sets  to  a  single 
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set,  it  is  desirable  to  obtain  the  best  target  samples  from 
each  of  the  sets.  This  is  done  by  reference  to  the  noise 
history  of  each  window.  Those  windows  with  low  drift  or  DC 
noise  are  are  averaged  with  corresponding  windows  in  other 
sets.  In  this  way,  biases  due  to  poor  target-background 
contrast  are  averaged  out.  The  remaining  window  set  is 
actually  a  number  of  vectors  to  be  combined  into  a  single 
vector.  The  sequencing  of  the  vector  (from  low  coefficients 
to  high)  is  a  critical  aspect  of  the  vector  recombination. 
There  are  several  ways  to  order  the  vectors.  One  method  is 
to  assume  the  longest  vector  is  first  with  gradually  shorter 
vectors  following  subsequent.  After  ordering,  the  resulting 
vector  is  the  generated  by  a  round-robin  extraction  from  each 
vector  in  the  window  set.  Another  method  is  to  order  the 
vectors  according  to  the  value  of  the  first  coefficient  of 
each  window  vector  in  the  set.  A  smaller  coefficient 
indicates  the  energy  in  the  vector  may  be  spread  out  over  a 
wider  area.  Once  the  vector  is  generated  and  normalized  it 
is  ready  for  deconvolution. 

The  deconvolution  process  attempts  to  remove  both 
the  finite  window  sampling  blur  and  the  optical  blur.  When 
the  detectors  comprising  the  windows  are  sampled,  the 
detector  aperture  forces  finite  sampling  window  effects  onto 
the  data.  The  width  of  the  detector  window  aperture  is 
dependent  on  the  number  of  samples  (or  frames)  in  which  a 
point  on  the  target  is  visible  in  the  detector  field  of  view. 
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Since  the  length  of  the  window  is  represented  digitally,  it 
must  be  an  integer  value.  Thus,  pixel  samples  which  are 
fractional  use  the  window  set  approach  outlined  above  to 
generate  an  integer  length  sampling  window. 

The  optical  transfer  function  (OTF)  is  modeled  here  as  a 
two-dimensional  Gaussian  function.  This  function  is 
separable  as  X  and  Y  functions  complimenting  the  windowing 
process.  With  the  OTF  modeled  as  a  one-dimensional  Gaussian, 
the  transfer  function  of  both  the  sampling  window  and  the 
OTF  is  represented  by  the  convolution  of  the  OTF  and  the 
sampling  window.  This  result  is  the  deconvolution  window 
used  to  recover  the  target  signal.  The  deconvolution  is 
performed  in  the  Fourier  domain  for  simplification  and 
at  consists  of  the  division  of  the  spectrum  of  the  target  signal 

by  the  spectrum  of  the  deconvolution  window.  This  result  is 
then  low-pass  filtered  in  the  Fourier  domain  to  remove  high 
frequency  noise  induced  from  external  sources  and  anomalies 
from  the  feature  extraction  process  itself.  The  calculated 
target  spectrum  is  then  inverse  transformed  to  obtain  a  time- 
domain  target  signal  representation.  A  flow  chart  of  the 
feature  vector  process  is  shown  in  Figure  4-7. 

The  feature  extraction  process  has  several  advantages. 
First,  the  data  can  be  acquired  in  real-time  since  no  pre¬ 
processing  is  necessary  to  collect  the  data.  Second, 
extraction  of  the  data  naturally  decomposes  two-dimensional 
functions  into  their  X  and  Y  components.  This  greatly 
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simplifies  the  deconvolution  process.  Third,  memory 
requirements  are  minimal  even  with  large  numbers  of  windows. 
Finally,  depending  on  the  scenario,  the  necessary  sampling  or 
observation  time  may  be  10  to  40  frames.  This  equates  to  250 
milliseconds  to  approximately  1  second  for  slowly  moving 
targets . 

Target  Classification 

The  target  feature  vectors  extracted  from  sensor  data 
are  compared  to  a  database  of  target  signatures  for 
identification.  Access  to  the  database  is  controlled  by  the 
target  orientation  (determined  from  observed  target  motion 
relative  to  the  X-Y  focal  plane  geometry.)  The  success  of 
the  classification  scheme  is  dependent  on  the  method  by  which 
the  potential  target  data  is  compared  with  that  of  the 
database.  This  study  tests  four  methods  to  compare  data. 
All  of  the  techniques  use  some  sort  of  normalization  and/or 
interpolation  to  perform  the  comparison.  For  this 
discussion,  model  data  will  refer  to  signatures  from  the 
target  database,  and  processed  data  will  refer  to  the  feature 
vectors  of  a  potential  target  extracted  from  sensor  data. 

The  target  database  is  constructed  from  data  taken  from 
the  five  model  aircraft  discussed  in  chapter  II.  Each  of  the 
targets  are  digitized  at  a  high  resolution  and  at  discrete 
orientations.  The  target  data  is  then  summed  into  rows  and 


columns  to  form  model  feature  vectors  and  saved  in  a  database 
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file.  To  perform  target  classification,  the  input  process 
data  is  compared  to  each  of  the  database  models  according  to 
the  processed  target  orientation.  The  model  which  has  the 
minimum  deviation  from  the  processed  data  is  judged  to  be  the 
processed  target.  Examples  of  the  target  data  base  are  shown 
in  Appendix  A.  Each  target  X-Y  projection  is  shown  at  ten 
orientations  from  0  to  90  degrees. 

The  methods  used  to  compare  the  model  and  processed 
data  evolved  from  the  nature  of  the  data  itself.  Recall  that 
each  point  in  the  feature  vector  represents  an  integration 
over  a  area  specified  as  the  resolution  of  the  data.  To 
compare  the  high  resolution  database  models  to  the  low- 
resolution  processed  data,  normalization  is  required  in  both 
magnitude  and  length  for  an  accurate  comparison.  The 
techniques  used  to  accomplish  the  comparison  are  a  normalized 
polynomial  interpolation,  a  matched  integration,  and  a 
integrated  linear  interpolation. 

The  Polynomial  Interpolation  method  uses  the  Newton 
form  of  the  interpolating  polynomial  to  obtain  a  polynomial 
fit  of  of  model  data  whose  order  is  equal  to  the  length  of 
the  model  vector.  The  method  bases  its  comparison  of  the 
processed  data  to  the  model  data  entirely  on  the  shape  of  the 
curve.  (All  data  is  normalized  in  magnitude  and  length  prior 
to  interpolation.)  Disadvantages  of  this  method  are  the 
erratic  nature  of  the  polynomial  at  higher  orders,  and  the 
elimination  of  length  as  a  feature. 
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The  Matched  Integration  method  performs  a  comparison  by 
integrating  the  model  data  to  match  as  closely  as  possible, 
the  resolution  of  the  processed  data.  To  obtain  the  proper 
integration  interval,  the  method  tries  to  match  the  first 
point  of  the  normalized  processed  data  with  successive 
integrations  of  the  model  data  until  the  first  points  are 
relatively  the  same.  Next,  the  remaining  model  data  is 
integrated  to  match  (as  closely  as  possible)  the  resolution 
of  the  processed  data  and  the  data  is  then  compared  point-by¬ 
point.  The  difference  between  the  model  and  processed  data 
indicates  the  degree  of  dissimilarity. 

The  Integrated  Linear  Interpolation  method  combines  the 
best  features  of  the  two  previous  methods.  The  model  data 
is  first  integrated  to  match  the  resolution  of  the  processed 
data.  This  consists  of  adding  model  data  points  to  match  a 
single  processed  data  point  and  repeating.  For  example, 
given  a  model  database  with  one  meter  resolution  and  a 
processed  target  which  is  extracted  at  a  2.25  meter 
resolution,  the  model  data  is  integrated  to  generate  a  2 
meter  model  target  (every  two  points  are  added  together.)  To 
correct  for  the  resolution  spacing  difference  between  the 
model  and  processed  targets  (2.25  verses  2),  the  processed 
data  is  linearly  interpolated  at  points  defined  by  the  model 
data.  Therefore,  the  processed  target  will  have  data  points 
at  2.25,  4.5,  7.75,  etc.  and  the  model  target  will  have  data 
points  at  2,  4,  6,  etc.  meters.  To  find  the  distance 
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between  targets  for  example,  the  processed  data  is 
interpolated  between  points  2.25  and  4.5,  at  point  4  given  by 


the  model  data.  Again,  the  difference  between  poi 

(mean-square  error)  indicates  the  closeness  of  the  data, 
overview  of  this  process  is  shown  in  Figure  4-8. 
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Figure  4-8.  Overview  of  Linear  Integrated  Interpolation 


V.  Software  Simulation  Description 


Introduction 

To  simulate  data  produced  by  a  staring  infrared  sensor, 
and  to  provide  a  software  testbed  for  classification 
algorithm  development  and  test,  a  computer  simulation  was 
developed  that  incorporates  key  features  of  staring  sensor 
performance.  This  chapter  presents  a  description  of  the  end- 
to-end  software  simulation  used  to  implement  and  test 
algorithms  discussed  in  previous  chapters.  The  simulation  is 
designed  to  allow  user-tai 1 oring  of  scenario  dependent 
parameters  and  is  useful  for  simulating  various  types  of 
targets  and  backgrounds  both  statically  and  dynamically.  The 
simulation  also  provides  flexibility  in  analyzing  algorithm 
performance  at  various  points  in  the  data  processing  path.  An 
overview  of  the  end-to-end  simulation  is  shown  below. 


Target  -> 
Model 


Backgnd 
Mode  1 


Figure  5-1.  Software  Simulation  Overview 
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VI.  Test  Results 


Introduction 

Tests  were  conducted  using  the  software  simulation  to 
determine  the  feasibility  of  target  classif ication  in  a  low- 
resolution  scenario.  It  is  not  expected  that  the  test  cases 
presented  here  will  fully  characterize  the  extraction 
process,  rather,  it  is  the  intent  to  obtain  a  general 
description  of  the  algorithm  performance  by  varying  certain 
target  and  sensor  parameters.  Given  this  philosophy,  limited 
testing  was  accomplished  by  gener-'  .  ng  scenes  whose 
parameters  were  defined  to  test  import,  it  aspects  of  the 
target  classifier.  This  chapter  will  discuss  the  test 
conditions  and  resulting  data  with  respect  to  the  performance 
of  thetarget  classifier. 

Scene  Generation 

For  each  test  case,  a  scene  sequence  is  generated 
consisting  of  20  to  50  frames  of  sensor  data  depending  on  the 
target  speed,  sensor  footprint,  and  extent  of  the  focal 
plane.  By  sequential  movement  of  the  target  (and  sometimes 
sensor)  and  application  of  system  transfer  functions,  scenes 
are  generated  according  to  user  specified  parameters. 
The  quality  of  feature  vector  extracted  from  sensor  data  is 
dependent  on  the  prevalent  noise  and  background  conditions, 
the  sensor,  and  motion  of  the  target.  It  is  therefore 
desirable  to  characterize  the  classifier  by  varying 
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closer  match  of  the  respective  target.  The  model  database 
target  are  indicated  on  the  horizontal  axis.  Since  this  is  a 
optimum  situation  (no  background  clutter),  these  results  can 
be  used  as  a  measure  for  the  tests  described  in  Chapter  VI. 


are  used  to  compare  data,  the  results  from  routines  INTERP 
and  INTEGRATE  are  included  for  reference  purposes  only.  Both 
these  routines  are  superseded  by  the  LINTEGRATE  routine  which 
performs  an  Integrated  Linear  Interpolation  (see  Chapter  4, 
page  54)  to  compare  data.  This  method  is  less  sensitive  to 
scale  and  energy  variations,  and  provides  a  more  reliable 
comparison  than  either  of  the  previous  methods.  The 
LINTEGRATE  method  was  used  to  generate  the  data  reported  in 
Chapter  VI. 

The  target  database  used  by  COMPARE  is  created  by  a 
utility  software  program  call  TARGEVAL  for  TARGet  EVALuator. 
This  program  takes  a  2-D  target  file  generated  from  TGEN  and 
calculates  the  associated  target  feature  vectors  for  the 
target.  The  target  feature  vectors  are  then  stored  according 
to  target  orientation.  This  process  is  repeated  for  each 
target  orientation  from  0  to  45  degrees  in  5  degree 
increments.  The  targets  used  for  the  data  base  have  a  1 
meter  pixel  resolution  and  are  based  on  high  resolution 
target  images  similar  to  larger  target  images  shown  in 
Figures  2-5  through  2-7. 

To  demonstrate  the  separability  of  targets  in  the  target 
database,  five  scenes  were  generated  using  each  of  the 
targets  without  a  background  (zero  energy.)  The  results  are 
shown  in  Figure  5-5.  The  vertical  axis  indicates  the 
deviation  of  the  extracted  target  to  each  of  the  model 
targets  in  the  database.  A  lower  deviation  indicates  a 
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number  samples  per  window  set)  and  generating  the  associated 
window  spectrum.  The  deconvolution  is  accomplished  by 
dividing  the  spectrum  of  the  input  vector  with  the  spectrum 
of  the  window.  Following  deconvolution,  DECON  filters  high 
frequency  components  of  the  target  spectrum  by  zeroing  the 
high  frequency  terms  in  both  the  real  and  imaginary  target 
data.  The  target  spectrum  is  then  inverse  transformed  and 
saved  in  a  file  for  target  classification.  1 

Target  Classifier 

The  target  classifier  is  implemented  in  a  program 
called  COMPARE  with  three  subroutines  which  perform 
comparisons  according  to  the  three  methods  described  in 
chapter  VI.  The  program  begins  by  reading  a  processed  target 
file  created  by  PROFILE.  After  the  data  is  normalized, 

COMPARE  calls  the  three  subroutines  INTERP,  INTEGRATE,  and 
L INTEGRATE  to  compare  the  input  data  to  each  of  the  model 
target  types  contained  in  the  target  database.  The  model 
database  is  accessed  according  to  the  input  target 
orientation.  The  results  of  each  comparison  are  stored  in  a 
table  and  presented  to  the  user  at  the  conclusion  of  the 
program. 

The  results  of  COMPARE  are  given  for  each  target  type 
according  to  the  deviation  of  the  input  target  to  each  target 
in  the  database.  COMPARE  uses  both  X  and  Y  target  feature 
vectors  if  present  in  the  input  data.  Although  three  methods 
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or  columns  of  detectors)  are  summed  and  placed  in  a  buffer. 
Since  this  is  the  only  operation  required  during  scene  data 
processing,  the  feature  extraction  process  continues  at  high 
speed  with  little  processing  overhead.  It  is  important  that 
the  windows  are  setup  far  enough  in  advance  of  the  target  to 
obtain  samples  of  the  background  prior  to  target  passage. 
These  background  samples  are  also  saved  in  a  buffer  to  allow 
calculation  of  D.C.  and  drift  induced  background  components 
for  removal  during  vector  processing. 

Vector  processing  is  accomplished  by  the  profile 
subroutine  VPROC.  Currently  VPROC  separates  only  D.C.  noise 
components  from  the  windowed  data.  After  removal  of  noise 
components,  windows  are  combined  and  averaged  according  to 
the  sampling  resolution  to  be  achieved  (see  chapter  IV.)  The 
data  is  then  reduced  to  a  single  feature  vector  and  is  ready 
for  deconvolution. 

The  sensor  induced  effects  are  removed  with  a 
inverse  filter  generated  by  prior  knowledge  of  the  OTF  and 
the  sampling  window.  The  coefficients  of  the  OTF  function 
are  generated  from  the  routine  GCOEFF  (also  resident  in  the 
scene  simulator)  and  passed  to  the  routine  DECON  which 
constructs  the  deconvolution  window  and  manages  the 
deconvolution  process  through  calls  to  a  one-d i mens  ions  1 
fourier  transform  routine  called  "FOUREA".  The  deconvolution 
window  is  constructing  by  convolving  the  OTF  function  with 
the  rectangular  sampling  window  (whose  length  is  equal  to  the 
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detectors  are  possible.  For  simulations  used  to  test  the 
extraction  algorithms,  10x10  and  20x20  detector  focal  planes 
were  used  to  insure  sufficient  subpixel  resolution. 

A  typical  output  of  the  scene  simulator  is  shown  in 
Figures  5-3  and  5-4.  The  high  resolution  scene  after  target 
masking  is  shown  beside  a  low-resolution  sensor  focal  plane 
output  illustrating  the  averaging  effect  of  the  detector 
array.  The  resolution  of  the  high  resolution  scene  is  2 
meters  per  pixel  and  the  low  resolution  pixel  is  10  meters. 
The  next  figure  illustrates  the  same  scene  after  addition  of 
the  gaussian  optics  blurring  function. 

Feature  Extractor 

The  feature  extraction  algorithm  described  in  chapter  4 
is  implemented  in  a  program  named  "PROFILE".  PROFILE  can 
process  targets  moving  in  any  direction  (0  -  90  degrees)  with 
minimal  hand-off  information  from  the  target  detection 
scheme.  PROFILE  currently  only  accepts  data  generated  form 
the  sensor  simulation  SCENE,  however,  any  data  which  is 
preprocessed  to  include  a  file  header  may  be  used.  There  is 
no  user-required  information  required  to  run  the  program 
since  all  needed  information  is  contained  in  the  scene  file 
header . 

PROFILE  begins  by  calculating  the  target  velocity 
vector  and  the  sensor  dependent  parameters  needed  for  window 
setup.  After  window  initialization,  data  is  processed  frame 
by  frame.  The  X  and  Y  windows  (consisting  of  contiguousrows 


target  data  file.  Thus,  a  target  with  a  pixel  resolution  of 
2  meters  input  into  the  scene  simulator  which  is  setup  with  a 
subpixel  resolution  of  15x15  defines  the  sensor  footprint  (L) 
to  be  30  meters.  The  background  used  in  this  study  is 
considered  to  have  no  scale  limitations.  Since  the  sensor 
footprint  consists  of  an  integer  number  of  subpixels  and  is 
therefore  not  continuous,  the  selection  of  target  speed  and 
sensor  motion  is  necessarily  limited  to  those  imposed  by  the 
subpixel  resolution. 

The  background  scene  used  extensively  in  this  study 
consists  of  a  photograph  of  a  cloud  layer  as  viewed  from 
above.  This  scene  is  especially  suitable  for  sensor 
evaluation  since  it  contains  both  high  and  low  spatial 
frequencies  in  an  anisotropic  pattern.  Also,  a  target  of  any 
size  can  be  superimposed  without  regard  to  the  scale  of  the 
background. 

Each  scene  file  generating  by  SCENE  contains  a  file 
header  describing  the  conditions  under  which  the  scene  file 
was  generated.  The  number  of  frames  generated  by  the  scene 
simulation  depends  on  the  relative  dwell  time  of  the  target 
and  is  usually  specified  so  that  the  target  crosses  the 
entire  sensor  field  of  view  during  the  simulation  run.  Focal 
plane  sizes  are  user  selectable  , however,  the  the  number  of 
detectors  and  associated  subpixels  are  limited  to  the 
absolute  size  of  the  background  scene  after  accounting  for 
sensor  motion.  Thus,  focal  plane  sizes  from  10x10  to  100x100 


62 


through  a  window  whose  size  is  defined  by  the  focal  plane 
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detector  and  sub-detector  (pixel)  dimensions.  The  window 
location  may  be  fixed  or  variable  depending  of  the  user- 
selected  drift  rate  parameter.  Sensor-induced  jitter  may  be 
modeled  in  a  similar  manner  although  is  not  implemented  in 
the  current  version  of  SCENE.  The  target  is  then  masked  onto 
the  background  according  to  its  current  location  in  the  scene 
sequence.  If  optics  are  included  in  the  simulation,  the  OTF 
or  poii.  t-spread  function  is  calculated  by  convolving  a 
averaging  function  whose  coefficients  are  determined  by  a 
two-dimensional  gaussian  function  generator  implemented  in  a 
subroutine  called  GCOEFF.  The  OTF  function  is  generated 
according  to  the  pre-defined  subpixel  dimensions  and  the 
desired  blur  radius.  For  tests  described  in  chapter  six,  a 
blur  radius  or  variance  of  the  gaussian  function  is 
approximately  equal  to  the  detector  size  for  a  matched-optics 
configuration.  Next,  the  high  resolution  scene  is  integrated 
over  the  logical  focal  plane  dimensions  to  model  the  detector 
averaging  function.  The  mapping  of  the  input  high  resolution 
data  to  the  low  resolution  focal  plane  data  is  accomplished 
within  the  constraints  imposed  by  the  hardware.  Each  low 
resolution  pixel  is  composed  of  an  integer  array  of  high 
resolution  subpixels.  The  pixel  dimensions  are  defined  in 
odd  subpixel  sizes  such  as  5x5,  7x7,  15x15,  etc.  to  allow 
easy  computation  of  the  optical  PSF.  The  resolution  of  the 
subpixel  is  determined  from  the  resolution  of  the  input 
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version  which  sacr if  ices  display  capabilities  for  increased 
processing  speed.  User  input  data  is  shown  in  Table  5-1. 

After  the  user  enters  preliminary  information 
initializing  the  sensor  system  and  target  parameters,  SCENE 
calculates  sensor  and  target  position  and  motion  based  on 
hardware  pixel  resolution.  The  program  then  enters  a  loop  to 
generate  frame  sequences  using  the  flow  of  operations 
outlined  below  in  Figure  5-2. 


Figure  5-2.  SCENE  Data-Generation  Flow-Chart. 

The  program  reads  the  background  data  into  a  buffer 


(Ref  10).  This  program  allows  storage  of  large  data  arrays 
accessible  during  program  execution  that  could  not  otherwise 
fit  into  core  memory  (64K). 

Sensor  Simulation 

The  purpose  of  the  sensor  simulator  is  to  generate 
scenes  consisting  of  sequences  or  frames  of  data  representing 
the  output  of  a  staring  mosaic  sensor.  This  process  is 
essentially  a  integration  of  the  input  data  (background  and 
target)  into  low-resolution  data  as  seen  by  an  array  of 
infrared  detectors.  The  simulation  includes  key  features 
such  as  specif ication  of  sensor  geometry,  optics,  frame  rate, 
and  sensor  motion.  The  sensor  simulation  is  called  "SCENE" 
and  is  implemented  in  two  versions  for  use  on  different 
computers  available  in  the  lab;  a  Nova  version  which  displays 
the  generated  scene  data  on  a  video  monitor,  and  an  Eclipse 


Table  5-1. 

SCENE  User  Input  Data 


User  Input  Files 

Program  Parameters 

-  Target  (From  TGEN) 

-  Focal  Plane  Size 

-  Background  (Video  File) 

-  Subpixel  Size 

-  Sensor  Field-of-View 

User  Input  Data 

(Window  on  Background 

-  Frame  Rate 

Scene) 

-  Drift  Rate 

-  Optical  Transfer  Function 

Output  Data 

-  Blur  Radius 

-  Filename 

-  Target  Velocity 

-  Number  of  Frames  to 

-  Target  Direction 

Generate 

-  Signature  (Graylevel) 

-  Initial  Position 

target  pixel  resolution  to  be  compatible  with  the  sensor 
simulation.  This  reduced  target  represents  the  high 
resolution  target  as  viewed  by  the  staring  sensor.  The 
preprocessing  program  is  called  "TGEN"  or  Target  GENerator. 
Before  using  the  program,  a  picture  of  the  target  model  is 
acquired  with  a  TV  camera  and  stored  in  memory  of  the  Octek 
image  analyzer.  The  field-of-view  of  the  earn. era  is 
determined  by  the  scale  of  the  plane  model  and  is  represented 
in  meters.  Thus,  the  user  inputs  two  parameters  to  TGEN,  the 
camera  field-of-view  (meters),  and  the  desired  pixel 
resolution  (meters).  TGEN  will  then  reduce  and  separate  the 
target  to  create  the  target  file.  The  target  file  format  is 
comprised  of  a  file  header  and  target  data.  The  header 
information  consists  of  the  number  of  pixels  that  define  the 
target,  the  direction  the  target  is  facing  (0-90  degrees), 
and  the  pixel  resolution.  Data  is  stored  by  X,  Y  and 
intensity  coordinates.  The  files  created  by  TGEN  are  named 
according  to  the  plane  type,  direction,  and  resolution  to 
access  the  target  database.  The  program  uses  the  Octek  to 
process  the  target  and  is  not  readily  transportable. 

Backgrounds  used  in  the  sensor  simulation  are  digitized 
with  the  Octek  and  preprocessed  to  a  256  x  256  x  16  format 
for  greater  dynamic  range  in  the  scene  simulator.  The 
backgrounds  are  then  stored  as  video  files  for  access  by  the 
sensor  simulator.  Manipulation  of  video  files  is 
accomplished  using  a  virtual  memory  utility  called  PICBUF 


The  simulation  consists  of  five  major  components;  the 
tar ge t / background  preprocessor,  the  sensor  simulation,  the 
feature  extractor,  the  target  classifier,  and  the  utility 
software.  Tables  and  listings  of  software  described  in  this 
chapter  are  contained  in  Appendix  B.  The  utility  software  is 
dedicated  to  performing  specific  simulation  performance 
measurements  and  maintenance  of  the  simulation.  This 
simulation  configuration  provides  a  flexible  medium  for 
algorithm  evaluation  whose  results  will  be  discussed  in 
Chapter  VI. 

Target/ Background  Database 

The  targets  and  backgrounds  used  by  the  sensor 
simulator  are  generated  using  a  video  camera  digitizer.  The 
target  database  consists  of  five  aircraft,  a  B-52,  F-15,  KC- 
135,  DC-10,  and  a  MIG-21.  The  plane  models  are  1:144  scale 
with  the  exception  of  the  MIG-21  which  is  1:72  scale.  With 
the  known  scale,  the  actual  size  can  then  be  represented  in 
the  simulation.  The  planes  are  painted  flat-black  to 
eliminate  variations  in  intensity  of  the  plane  surface  during 
digitization.  Also,  the  plane  resembles  a  relatively  uniform 
silhouette  at  this  point  which  is  easily  modified  through 
software  to  represent  a  target  signature  of  interest.  The 
planes  are  digitized  in  a  256x256x4  pixel  frame.  This  raw 
image  is  then  preprocessed  to  generate  the  target  file. 

The  target  preprocessing  consists  of  reducing  the 


parameters  which  define  each  of  these  components. 

The  sensor  footprint  is  determined  by  the  resolution  of 
the  target  and  the  number  of  subpixels  per  detector.  Since 
the  overall  size  of  the  background  scene  is  limited  (256  x 
256),  the  resulting  focal  plane  size  is  likewise  constrained. 
The  choice  of  sensor  footprint  is  therefore  bounded  by  these 
restrictions  and  resulted  in  four  sensor  footprints  chosen  as 
sensor  parameters  for  testing.  The  table  below  shows  these 
footprints  along  with  resulting  focal  plane  and  subpixel 
dimensions . 


Table  6-1. 
Sensor  Footprints 


Footprint 

Focal  Plane 

Subpixel 

Size 

Dimensions 

Dimensions 

10  meters 

20  x  20  detectors 

5x5  pixels 

20 

10  x  10 

9x9 

35 

10  x  10 

15  x  15 

45 

10  x  10 

21  x  21 

Other  sensor  parameters  were  selected  based  on  the 
typical  sensor  of  interest  by  FTD  and  by  those  typically  used 
in  other  studies  (Ref  1,18).  Table  6-2  shows  sensor 
specif ications  used  to  define  the  sensor  model  for  the  sensor 
simulation. 
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Detectors 


Sample  Rate 


-  Uniform  Response  -  Frame  Rate  40  Hz 

-  8-12  micron 

-  Background  Limited 

Optics  Configuration 

-  Gaussian  Point  Spread  -  Staring  Mosaic 

-  Matched  Optic/Detector  -  Nadir  View 

-  100%  Transmission 

Table  6-2.  Sensor  Specifications. 

A  target  is  specified  for  a  given  test  case  by  its  type, 
direction,  velocity  and  intensity. 

Target  Velocities  were  chosen  according  to  subpixel 
resolution  and  vary  over  a  range  of  125  to  700  meters  per 
second.  Although  several  of  the  actual  targets  represented 
in  the  simulation  are  incapable  of  speeds  in  this  range,  the 
velocities  are  for  characterization  purpose  only.  The 
relative  phasing  of  the  target  is  controlled  by  the  choice  of 
target  velocity.  Most  tests  used  both  in-phase  and  out-of¬ 
phase  velocities  for  characterization. 

Target  Direction  is  important  in  specifying  the  relative 
SFOV  crossing  geometry  and  the  number  and  length  of  the 
feature  vectors.  For  target  directions  from  15-75  degrees, 
two  vectors  are  usually  produced.  For  directions 
approximately  parallel  to  the  detector  geometry,  only  one 
vector  is  produced.  Intermediate  directions  within  the  15-75 
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degree  range  result  in  changes  in  the  X-Y  velocity  components 
and  can  be  represented  by  either  the  velocity  or  direction 
parameters.  Therefore,  only  two  directions  were  used  in  the 
tests,  0  and  45  degrees.  The  0  degree  direction  generates  a 
single  feature  vector,  and  should  result  in  a  decrease  in  the 
classification  ability  of  the  algorithms.  The  45  degree 
direction  generates  two  vectors  each  of  which  has  the  same 
feature  resolution. 

Target  Signal  Intensity  is  specified  by  thresholding  the 
target  signature  at  a  certain  graylevel.  This  results  in  a 
relative  increase  in  the  signal  from  the  target,  and  an 
increase  in  the  contrast  between  the  target  and  background. 
The  contrast  between  the  target  and  background  is  represented 
by  a  contrast  ratio  defined  by  the  signal  due  to  the 
background  over  the  local  area  in  which  the  target  is 
observed,  and  the  target  signal  itself.  This  value  is 
analogous  to  signal-to-noise  ratio  and  is  intended  to  examine 
the  performance  of  the  target  classifier  over  a  range  of 
target  and  background  signal  levels. 

Test  Cases 

Trials  were  conducted  at  each  footprint  over  a  range  of 
velocities  and  directions.  The  objectives  of  the  test  cases 
are  oriented  toward  examining  target  separability  with 
respect  to  different  target-sensor  scenarios.  Specifically, 
trials  were  run  to: 
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1. 


Determine  the  separation  of  target  classes  with 


decreasing  information  extracted  for  the  target. 

2.  Determine  the  classifier  sensitivity  to  orientation 
and  velocity  of  targets. 

3.  Determine  classifier  performance  with  decreasing 
contrast  or  reduced  target  signal  information. 

4.  Examine  the  effect  of  increased  sensor  footprint  on 
target  extractor  performance. 

5.  Examine  sensor  motion  influence  on  the  feature 
extraction  process. 

To  satisfy  these  test  objectives,  three  categories  of 
tests  were  performed.  The  first  group  of  tests  examined  the 
effect  of  target  motion  on  classifier  performance  relative  to 
the  SFOV  at  varying  velocities  and  direction.  This  was 
repeated  for  each  target  and  at  the  four  footprints  listed  in 
Table  6-1.  The  second  group  of  tests  consisted  of  the  five 
targets  flown  with  an  increasing  contrast  ratio.  Each  target 
traversed  the  SFOV  at  45  degrees  and  at  250  meter s / second. 
The  third  test  involved  drifting  of  the  sensor  during  target 
observation.  This  was  performed  with  the  same  target 
parameters  as  Group  II.  All  tests  used  the  cloud  background 
discussed  in  Chapter  5  except  Group  II  which  uses  a  single 
graylevel  background  for  several  tests. 
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Table  6-3. 
Summary  of  Tests 

■_ 

Test 

Target 

Vel 

Dir 

Contrast 

Footprint 

m 

► 

Group 

I 

All 

125 

to 

700  mps 

0,45  deg 

HI 

11  meters 

20 

35 

45 

• .  - 

Group 

II 

All 

250 

45 

VARIABLE 

20 

m 

Group 

III 

B-52 

250 

45 

HI 

20 

ra  id. 
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The  pass/fail  criteria  for  all  of  the  tests  are  based 
on  the  target  classifier  results.  The  target  classifier 
generates  output  for  each  test  consisting  of  the  relative 
deviation  of  the  feature  set  under  investigation  and  the 
target  database.  For  targets  that  are  successfully 
extracted,  it  is  desirable  to  have  a  measure  of  the 
"confidence"  of  the  target  estimate.  The  measure  used  here 
consists  of  the  interval  between  the  correct  target  choice 
and  the  next  closest  choice  as  specified  by  the  classifier. 
Therefore,  an  extracted  target  which  compares  similarly  with 
two  targets  in  the  database  should  have  a  low  confidence 
estimate,  and  a  target  that  has  no  deviation  from  the 
corresponding  database  target  while  maintaining  separation 
with  other  model  targets  in  the  database,  should  have  a  high 
confidence.  A  mathematical  figure  of  merit  called  a 
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Confidence  Distance  is  then  defined  as  the  separation  between 
the  correct  target  choice  and  the  next  possible  choice  given 
as  a  percentage: 


Confidence  Distance(%)  = 


i 


-[Correct  Target  Dev 
Closest  Target  Dev 


00 


For  Confidence  Distances  (CD) : 

<  1  :  Target  Classified  Incorrectly 
=  0  :  Target  does  not  separate 
>  1  :  Target  Classified  Correctly 


Analysis  of  the  relative  "correctness"  of  the  overall  target 
extraction/classification  will  be  in  terms  of  the  Confidence 
Distance  above.  Since  there  are  actually  three  comparison 
methods  used  in  the  classification  (Reference  Chapter  4),  the 
method  which  gave  the  most  consistent  results  is  used  for  the 
Confidence  Distance  measurement.  This  method  is  the 
Linear  Integrated  Interpolation  which  generated  all  of  the 
test  results  reported  here.  The  other  methods  did  produce 
greater  target  separability  for  some  tests,  however,  the 
Linear  Integrated  Interpolation  method  was  by  far  the  most 
reliable. 


Results  of  Group  Tests 

Approximately  175  trials  were  conducted  (limited  by 
computer  time)  with  various  sensor  and  target  parameters 
defined  for  the  Group  I  tests.  In  all  but  four  of  the  tests 
the  target  was  successfully  extracted  and  classified  from  the 
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simulated  data.  The  four  tests  that  failed  were  high  speed 
tests  (greater  than  600  meters/second)  which  resulted  in  only 
a  small  number  of  points  defining  the  feature  vector  of  the 
extracted  target.  However,  even  with  an  incorrect  plane 
classification,  it  was  still  possible  to  separate  plane 
classes.  For  example,  a  F-15  may  be  incorrectly  classified 
as  a  MIG-21,  or  the  KC-135  as  a  B-52.  There  were  no  tests 
which  the  separability  of  the  target  classes  (transport 
versus  fighter)  was  not  maintained.  With  respect  to 
individual  aircraft,  the  larger  aircraft  separated  with  a 
higher  confidence  than  smaller  aircraft.  This  remained  true 
at  single  pixel  resolutions  for  each  of  the  aircraft.  This  is 
mainly  due  to  differences  in  wing  structure  between  the  large 
aircraft,  and  the  larger  number  of  samples  that  define  the 
large  aircraft.  The  results  of  a  trial  using  a  20  meter 
footprint  and  out-of-phase  target  velocities  of  250  meters/ 
second  is  depicted  in  Figure  6-1.  The  figure  can  be  compared 
with  Figure  5-5  on  page  68.  Note  the  decrease  in  relative 
separation  distance  especially  on  the  fighter  aircraft.  This 
is  due  to  degradation  effects  of  the  background. 

At  higher  target  velocities,  the  resolution  of  target 
feature  vector  is  limited  to  several  data  points.  Therefore, 
the  upper  bound  for  the  utility  of  the  target  classifier 
defined  by  the  target  speed  with  respect  to  the  sensor. 
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sampling  rate.  For  these  tests,  target  speeds  above  600 
meters/ second  resulted  in  a  rolloff  of  the  classification 


confidence  interval.  These  results  are  shown  in  Figure  6-2, 
where  the  five  aircraft  generated  Confidence  Distances  over 
the  velocity  range  indicated.  At  higher  velocities,  the 
relative  phasing  of  the  target  becomes  much  more  important. 
For  out-of-phase  targets,  the  combination  of  many  frames  of 
data  results  in  higher  resolution  feature  vectors  and 
therefore,  an  improved  Confidence  Distance.  This  can  be  seen 


as  a  peak  in  the  aircraft  Confidence  Distances  at  460 
meters/second  in  Figure  6-2.  The  five  aircraft  are  shown  over 
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Figure  6-2.  Rolloff  of  The  Confidence  Distances  of 
Five  Aircraft  with  Increasing  Velocity 


a  range  of  velocities  sampled  with  a  20  meter  footprint.  The 
effectiveness  of  the  technique  diminishes  as  the  target  moves 
toward  an  in-phase  condition.  Therefore,  an  in-phase  target 
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with  a  velocity  of  600  m/s  may  be  incorrectly  classified, 
while  correctly  classified  when  out-of-phase  at  650  m/s. 

The  target  classifier  was  quite  effective  in  comparing 
the  extracted  signature  to  the  database  model  samples.  The 
Linear  Integrated  Interpolation  method  provided  the  best 
results  in  comparing  the  data.  The  key  features  of  the 
target  which  allow  target  separation  are  the  wing  structure 
and  target  length.  The  absolute  energy  in  the  extracted 
signatures  was  not  a  factor  in  separating  the  targets  after 
the  signature  increased  above  a  certain  contrast  threshold 
as  discussed  in  the  next  section.  Although  it  may  seem  that 
leng tha 1  one  may  be  sufficient  to  separate  the  aircraft  used 
here,  successful  classification  of  extracted  targets  with 
resolutions  of  5  meters  demonstrate  that  wing  shape  is  a 
dominate  factor.  Specifically,  length  indicates  the 
difference  between  target  classes  (B-52  versus  a  F-15),  and 
wing  -shape  differentiates  between  intraclass  aircraft  (B-52 
versus  a  KC-135). 

From  the  analysis  performed  in  Chapter  II,  the  target 
velocity  versus  the  system  sampling  function  define  the 
resolution  of  the  sampled  data.  This  would  indicate  that  the 
sensor  footprint  resolution  is  not  a  major  factor  in  defining 
the  performance  of  the  classifier.  Therefore,  the 
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performance  of  the  classifier  with  a  target  at  a  10  meter 
footprint  should  not  be  different  than  the  performance  with 
the  same  target  at  50  meters.  However,  the  system  transfer 
functions  and  the  associated  vector  processing  seem  to  limit 
extrapolation  of  the  test  results  to  higher  footprint  sizes. 
Data  from  the  tests  indicate  that  there  is  a  rolloff  in  the 
confidence  intervals  for  all  targets  as  the  footprint  size 
increases.  An  example  of  this  is  shown  in  Figure  6-3.  This 
suggests  that  the  increased  losses  from  sensor  transfer 
functions  at  higher  footprints  becomes  more  difficult  to 
remove.  At  targets  resolutions  below  one  pixel,  the  blur  or 


Figure  6-3.  Rolloff  of  Confidence  Distance  Versus  Footprint 
For  a B-52  (velocity  =  250  meter s/ second  against 
a  cloud  background) 
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low  pass  filter  effect  of  the  system  reduces  the  information 
content  of  the  target  signature  to  a  point  were  deconvolved 
becomes  ineffective  in  recovering  the  target  signal.  In 
fact,  for  sensors  with  a  large  blur  radius  (DFOV<  02),  the 
deconvolution  process  tends  to  distort  the  recovered  time- 
domain  feature  vector.  Another  problem  with  extending  this 
research  to  higher  footprints  is  the  limitation  of  dynamic 
range  imposed  by  the  system.  The  incremental  energy 
difference  from  target  motion  across  larger  footprints  is 
much  less  than  that  of  smaller  footprints  with  similar 
targets.  Therefore,  the  target  information  is  more  difficult 
to  extract  and  more  susceptible  to  temporal  noise  in  the 
respective  DFOV.  Therefore,  from  the  results  of  these  tests, 
an  upper  limit  to  the  target  classifier  scene  data  is  on  the 
order  of  a  single  pixel  on  a  target. 

The  results  of  the  Group  I  tests  indicate  that  the 
targets  are  separable  within  the  velocities,  directions,  and 
resolutions  tested.  The  performance  of  the  classifier  met 
expectations  with  relatively  low  error  rates.  The 
performance  is  bounded  for  decreasing  input  data  resolutions 
and  high  target  velocities. 

Results  of  Group  II  Tests 

The  Group  II  tests  measure  the  relative  contrast 
between  the  target  and  background  to  determine  the  effect  of 
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background  clutter  on  target  extraction.  Two  backgrounds  are 
used  in  this  test,  the  standard  cloud  scene,  and  a 
homogeneous  background  with  a  constant  graylevel.  The 
spatial  distribution  of  the  energy  in  the  background  cloud 
scene  is  anisotropic,  so  that  both  high  and  low  target 
contrasts  occur  during  target  observation.  This  background 
energy  distribution  is  particularly  suitable  to  demonstrate 
the  feature  extractors  ability  to  obtain  a  good  target 
feature  representation  in  a  variable  contrast  scenario.  The 
other  background  is  created  by  filling  memory  with  a  constant 
graylevel  value  of  7  (on  a  scale  from  0  to  15;  0  is  Black,  15 
is  White.) 

The  first  tests  determined  the  ability  of  the  feature 
extraction  method  to  obtain  good  samples  from  targets  with 
different  energy  contrasts.  By  increasing  the  target  energy 
threshold,  an  improvement  in  average  target/background 
contrast  is  achieved.  Tests  were  conducted  using  the 
variable  contrast  cloud  scene  and  the  single  contrast  gray 
scene.  Figure  6-4  shows  the  results  of  two  aircraft  against 
the  cloud  background.  With  reference  to  Figure  6-4,  a  lower 
target  deviation  indicates  a  closer  comparison  with  the 
respective  database  target.  The  target  energy  refers  to  the 
relative  contrast  between  the  target  and  background.  The 
data  shows  a  relatively  constant  deviation  for  targets  over 
the  contrast  range.  These  test  results  illustrate  the 
classifier’s  ability  to  separate  aircraft  in  a  variable 
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Figure  6-4.  Effects  of  Variable  Target  Energy  on  Target 

Extraction  Performance  (Top:KC-135,  Bot:  F-15) 
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contrast  background.  This  analysis  indicates  that  the 
feature  extraction  method  using  redundant  window  sampling, 
can  consistently  obtain  a  good  sample  of  the  target  signature 
given  the  opportunity  to  view  portions  of  target  at  high 
contrast.  Therefore,  this  method  takes  advantage  of  good 
contrast  in  selected  windows  to  generate  a  higher  quality 
feature  vector  then  could  otherwise  be  extracted  through  a 
single  window.  This  is  illustrated  in  Figure  6-5,  where  the 
confidence  intervals  for  the  scene  sequence  are  shown  given 
multi-window  versus  single  window  sampling  using  the  cloud 
background.  Note  the  Confidence  Distance  is  much  less  for 
the  B-52  sampled  thru  a  single  window. 

The  remaining  tests  used  the  single  gray-level 
background  versus  the  cloud  background.  A  KC-135  using  the 
same  scene  parameters  as  depicted  in  Figure  6-4,  was 

simulated  against  the  gray  background.  The  results  are  shown 
in  Figure  6-6. 

A  comparison  between  a  gradually  increasing  KC-135 
signature  against  both  backgrounds  (from  Figures  6-4  and  6-6) 
is  shown  in  Figure  6-7.  As  shown,  the  KC-135  with 
relatively  lower  target  energy  can  be  extracted  with  a  higher 
Confidence  Interval  when  viewed  against  a  variable  contrast 
background. 
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Figure  6-5.  Comparison  of  Single  versus  Multiple  Window 
Target  Extraction  (B-52  at  250  m/s) 
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Figure  6-7. 


Comparison  by  Confidence  Distance  of  a  KC-135 
Extracted  From  a  Cloud  versus  a  Constant 
Background 
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Group  III  Tests  -  Effects  of  Sensor  Motion 

The  influence  of  sensor  motion  on  target  extraction  is 
test  using  a  20  meter  footprint  and  a  sensor  drift  rate  of  90 
meters  per  second  (or  approximately  11%  DFOV  drift  per 
frame).  Figure  6-8  shows  a  comparison  of  data  extracted 
from  parallel  sampling  windows.  The  target  is  a  B-52  with  a 
velocity  of  250  meters/second.  The  first  set  of  curves  a  B- 
52  signature  with  no  background,  the  other  set  of  curves  show 
a  B-52  against  a  cloud  background.  Figure  6-9  shows  the  same 
target  flying  over  the  same  background  scene  while  the 
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Target  Extracted 
(No  Background) 


Target  Extracted 
Cloud  Background 
(No  Drift) 


Figure  6-8.  B-52  Signature  Extracted  from  Four  Parallel 

Windows  (Left:  No  Background,  Right:  Cloud 
Background  [no  drift] ) 


sensor  is  drifted.  The  severe  degradation  in  the  drifted 
data  is  apparent  and  the  classifier  failed  to  extractthe 
signature  successfully.  Unfortunately,  the  low  frequency 
changes  in  the  background  when  drifted  are  of  the  same 


>1*1 


Target  Extracted 
Cloud  Background 
(Drifted) _ 


Figure  6-9.  Effect  of  Sensor  Drift  on  Target  Signal 
Extracted  From  Parallel  Sampling  Windows 
(B-52  Against  a  Cloud  Background) 


The  current  target  extraction  algorithm  demonstrates 
severe  weakness  in  separating  a  moving  background  from  a 


moving  target.  The  effects  of  sensor  motion  require  further 
study  before  a  practical  application  of  the  feature 
extraction  method  can  be  implemented. 
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VII.  Conclusion  And  Recommendations 


<3* 


Conclusion 

The  methods  developed  in  this  study  to  extract  and 
classify  targets  from  simulated  low -resolution  staring 
infrared  sensor  data,  are  shown  to  be  successful  with  various 
sensor  and  target  scenarios.  The  classifier  uses  a  feature 
set  based  on  one-dimensional  feature  vectors  extracted  from 
sensor  data.  The  feature  vectors  consist  of  orthogonol 
projections  of  the  target  spatial  signature.  The  target 
information  is  extracted  through  logical  sampling  windows 
consisting  of  contiguous  rows  or  columns  of  detectors.  The 
windows  are  set-up  prior  to  target  passage  based  on  minimal 
hand-off  information  from  a  target  detection  scheme. 
Increased  resolution  is  accomplished  by  combining  information 
from  adjacent  windows  thereby  taking  advantage  of  the 
target's  relative  motion  with  respect  to  the  sensor. 

Results  of  performance  tests  indicate  that  aircraft  are 
separable  and  can  be  classified  with  high  confidence  using 
the  one-dimensional  feature  extraction  method.  Aircraft  were 
successfully  classified  with  single  pixel  resolutions. 
Extrapolation  of  the  results  of  this  study  are  limited  by 
the  sensor  optics  and  sampling  resolution.  The  loss  of 
information  from  the  sensor  optics  limit  the  size  of  the 
sensor  footprint  which  the  classifier  can  be  applied  (for  a 
matched  optic  configuration.)  Further  limitations  to 
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classifier  performance  depend  on  the  target  velocity  and 
sensor  sampling  rate.  The  loss  of  target  information  due  to 
poor  target  contrast  can  be  offset  by  multiple  looks  if  the 
background  is  anisotropic.  Sensor  motion  generates  severe 
degradations  in  the  extracted  data  to  the  extent  that 
classification  is  not  possible  without  preprocessing. 

The  algorithms  developed  in  this  study  can  be 
implemented  in  real-time  since  little  overhead  processing  is 
required.  Target  detection  handoff  information  is  limited  to 
target  position  and  direction.  The  extraction  method  can 
correct  for  errors  in  handoff  information  once  the  extraction 
process  is  initiated. 

Recommendations 

Since  this  effort  consists  of  a  target  classification 
feasibility  study  which  encompassed  a  broad  range  of  topics, 
the  depth  of  research  was  necessarily  limited.  Only  limited 
analysis  and  testing  could  be  accomplished  due  to  time 
constraints.  Therefore,  recommendations  for  follow-on  study 
consist  of  further  analysis  and  testing  with  regard  to 
understanding  the  physical  mechanisms  involved  in  loss  of 
target  information  in  a  staring  sensor  scenario,  and,  the 
methods  involved  in  countering  these  mechanisms  to  retrieve 
useful  target  information.  One  such  method  was  presented 
here,  developed  in  a  "first-order"  simplified  sensor-model 
simulation.  However,  extending  this  study  effort  will 
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require  enhanced  models  supported  by  solid  analysis  and 
testing.  Therefore,  specific  recommendations  are  oriented 
toward  understanding  the  staring  sensor  scenario,  and 
development  and  test  of  enhanced  extraction  methods  based  on 
the  improved  sensor  model: 

1.  Further  research  on  thermal  properties  of  typical 
infrared  targets  from  space  would  support  the  assumption  that 
a  finite  database  could  be  generated  for  desirable  targets. 
Since  target  graylevel  can  vary  as  a  function  of  ambient 
operating  conditions,  analysis  of  target  signature  variations 
under  different  viewing  conditions  could  establish  if  such  a 
database  is  possible. 

2.  As  shown  in  this  study,  sensor  motion  relative  to 
the  background  can  contribute  severe  distortion  in  the 
extracted  target  signal.  By  sampling  the  background  around 
the  target  sampling  window,  and  with  knowledge  of  the  sensor 
motion,  it  may  be  possible  to  estimate  the  motion  induced 
signal  variation.  Further  study  is  therefore  required  to 
determine  the  feasibility  of  separating  background  from 
target  signals. 

3.  The  sensor  model  used  in  this  study  approximated 
the  transfer  function  for  a  staring  infrared  sensor.  An 
improved  sensor  model  could  yield  more  conclusive  results. 
This  model  should  include  an  improved  detector  model 
to  account  for  non-uniformity,  finite  detector  response  time, 
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detector  integration  time,  and  detector  noise  constraints. 

4.  To  compliment  three,  the  background  scenes  viewed 
by  the  sensor  should  include  greater  dynamic  range  and  a 
variety  of  thermal  possibilities  (ocean,  mountains, 
farmland.)  Obtaining  real  infrared  data  could  enhance  the 
credibility  of  the  test  results. 

5.  This  study  used  the  one-dimensional  feature  vector 
to  describe  and  separate  aircraft  from  IR  data.  However, 
given  a  larger  number  of  target  possibilities,  ~  -  ture 
vector  separation  may  be  more  difficult.  Therefore, 
extensive  testing  using  a  large  number  of  aircraft  will  be 
required  to  determine  if  the  general  solutions  afforded  by 
the  one-dimensional  feature  vectors  described  in  this  study 
will  separate  aircraft. 

6.  It  would  be  desirable  to  augment  the  target 
extractor  software  with  a  target  detection  algorithm  so  that 
data  from  different  sources  could  be  tested  automatically 
without  user  interaction. 
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Appendix  A 


Target  Database 

The  target  database  used  to  classify  extracted  targets 
is  depicted  in  Figures  A-l  through  A-10  on  the  following 
pages.  The  five  targets  are  represented  as  X  and  Y 
projections,  each  having  10  different  orientations  with 
respect  to  the  vertical  axis.  The  vertical  axis  (top  of  the 
page)  is  designated  as  zero  degrees  with  angles  increasing 
clockwise  to  the  vertical  axis.  The  diagram  below 
illustrates  the  convention  used  to  represent  the  data: 


X-Axis 


Y-Axis 


Each  figure  contains  10  signature  graphs  which  represent 
the  target  signature  projected  against  the  indicated  axis. 
The  signature  graph  uses  normalized  enoigv  (vertical  axis), 
verses  samples  of  the  target  defined  in  one  motor  increments 
(horizontal  axi s)  . 
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^a******************************************************************* 

C 

C  PROGRAM  COMPARE  VERSION  1.1 

C 

C  FUNCTION:  THIS  PROGRAM  COMPARES  A  PROCESSED  X-Y  VECTOR 

C  GENERATED  BY  THE  RESOLUTION  ENHANCEMENT  ROUTINES 

C  AGAINST  A  DATA  BASE  OF  KNOWN  TARGET  MODELS. 

C 

C  Macro:  FORT  COMPARE 

C  RLDR  COMPARE  INTERP  INTEGRATE  LINTEG  FORT. LB 

C 

C  LANGUAGE:  FORTRAN  IV 

C 

PARAMETER  ASIZE=100 

REAL  RXMAX , RYMAX , RXDI S , RYD I V 
INTEGER  NUMPLANE 

REAL  XPROC (ASIZE) ,YPROC (ASIZE) ,XP INDEX (ASIZE) ,YP INDEX (ASIZE) 
REAL  XMODEL (ASIZE)  , YMODEL ( AS  I ZE )  , XM I NDEX ( AS  I ZE )  ,YMINDEX (ASIZE) 
REAL  XPTEST (ASIZE) , YPTEST ( AS I ZE ) ,XPIN(ASIZE) ,YPIN(ASIZE) ,MPS 
DIMENSION  INFILE (20) ,MODFILE (7) 

INITIALIZE  VARIABLES 

RXMAX=0 
RYMAX =0 
XPMAX=0 
YPMAX=0 
RXDI V  =  0 
RYDI V=  0 
NUMPLANES  =  5 
I  TEST  =  0 

ACCEPT " RUN  TEST?  ",ITEST 

IF  ( I TEST . NE . 1 )  GOTO  14 
TYPE 

TYPE"  (1)  INTERPOLATION" 

TYPE"  (2)  INTEGRATION" 

TYPE"  (3)  LINEAR  INTERPOLATION" 

TYPE 

ACCEPT"WHICH  ONE?  ",ITEST 
C  GET  PROCESSED  DATA 

14  ACCEPT"DO  YOU  WANT  A  PRINTOUT  OF  RESULTS?  ",IPR 

TYPE  "ENTER  PROCESSED  DATA  FILENAME:  " 

READ (11,16)  INFILE ( 1 ) 

16  FORMAT ( S2  0 ) 
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LINTEG.FR 


S,A  Calculates  Deviation  Feature  Target 

Between  Target  Vectors  Vectors  Dev 

w/Linear  Interpolation 

INTEGRATE  S,A  Calculates  Deviation  Feature  Target 

Between  Target  Vectors  Vectors  Dev 

w/Normalized  Integration 

INTERP.FR  S , A  Calculates  Deviation  Feature  Target 

Between  Target  Vectors  Vectors  Dev 

w/Polynomial  Interpolate 


Directory  Of  Appendix  B 

Program  Page 


COMPARE.FR . 110 

DDA'A.FP . 118 

DECON.FR . 119 

F0KT4T05.FR . 123 

FOUREA.FR . 124 

FRAC.FR . 126 

FRDUMP.FR . 127 

GCOEFF.FR . 129 

GETDA.FR  . 129 

GETHDR.FR . 130 

INTEGRATE. FR  .  131 

INTERP.FR . 133 

INVPLOT.FR . 135 

LINTEG.FR . 138 

PIXFRAME.FR . 141 

PIXMOD.FR . 142 

PROFILE  .  FR . 143 

PSPREAD.FR . 154 

SCENE.FR . 156 

TARGEVAL.FR . 162 

TGEN.FR . 170 

VPROC.FR . 173 


m 


GCOEFF.FR 

S,  A 

Generates  2-D 

Gaussian  Coefficients 

2-D 

Array 

Target  Extraction 

Program 

Type 

Function 

Inputs 

Output 

PROFILE.FR 

M ,  A 

Extracts  Target 

Signature  Sensor  Data 

IR  Data 

Target 

Feature 

Vectors 

decon.fr 

S ,  A 

Filters  Feature  Vectors 

In  Fourier  Domain 

X-Y 

Vectors 

Feature 

Vectors 

vprroc.fr 

S  /  A 

Reduces  Sampled  Window 
Data  to  X-Y  Vectors 

Window 

Data 

X-Y 

Vectors 

plotio.fr 

S  ,  A 

Plots  Vector  Data  Via 
Calls  to  CALC0MP 

— 

fourea.fr 

S ,  A 

Computes  1-D  Fourier 
Transform 

Time 

Freq 

Freq 

Time 

gethdr.fr 

S  t  A 

Reads  Header  From 

Sensor  Data  File 

— 

Header 

Data 

getda.fr 

S ,  A 

Reads  A  Frame  of  Sensor 
Data  From  File 

— 

IR  Data 
Frame 

ddata.fr 

5,  A 

Displays  Frames  of 
Low-Res  Data  On  Octek 

IR  Data 
Frame 

“ 

pixframe.fr 

S ,  A 

Calculates  Target 
Velocity  in  Pixels/Frame 

Target 

Parameters 

frac.fr 

S  ,  A 

Calculates  Out-of-Phase 
Resolution 

Target 

Parameters 

frdump.fr 

M ,  U 

Dumps  IR  Sensor  Data 
Created  By  SCENE  at 

Line  Printer 

SCENE 

Data 

Target  Classification 

Program 

Type 

Function 

Input 

Output 

compare.fr 

M,A 

Classifys  Extracted 
Targets  With  Database 

Feature 

Vectors 

Target 

Dev 
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Appendix  B 


Sof tware  Listings 

The  software  used  in  this  thesis  is  grouped  according  to 
application.  This  includes  target/background  generation, 
scene  simulation,  target  extraction,  and  target 
classification.  Within  these  categories  the  software  is 
typed  according  to  hierarchy  and  function.  An  application 
(A)  or  a  utility  (U)  program  may  be  a  main  routine  (M),  or  a 
subroutine  (S).  The  Software  is  written  in  FORTRAN  IV  and 
FORTRAN  V  for  use  on  the  Data  General  Nova  and  Eclipse 
computers  respectively. 


Tar get /Background 

Generation 

Program 

Type 

Function 

Input 

Output 

TGEN.FR 

M ,  A 

Generates  Target  Files 
Form  Video  Data 

— 

Target 

Files 

TARGEVAL.FR 

M,U 

Generates  Target  X-Y 
Vectors 

INVPLOT 

M ,  U 

Prints  Target  Database 

At  Line  Printer 

Target 

Files 

- 

P1XM0D.FR 

M,U 

Converts  Standard 

Video  Files  to  16  Bit 

256x256x4 

Video 

256x25 

xl6 

Scene  Generatio' 

Program 

Type 

Function 

Input 

Output 

SCENE.FR 

M ,  A 

Generates  Simulated 

Focal  Plane  Output 

Backgrnd 

Target 

Sensor 

Data 

pspread.fr 

S ,  A 

Convolves  2-D  Video 

With  2-D  Coefficents 

For  Point-Spread  Funct 

Video 

File  & 

PSF  Data 

Blured 

Scene 

Figure  A-10.  Y-Axis  MIG-21  Data  at  10  Degree  Increments 


Figure  A-8.  Y-Axis  F-15  Data  at  10  Degree  Increments' 


Figure  A-7.  X-Axis  F-15  Data  at  10  Degree  Increments 


( I**0IX> 


( itteix) 


< :*t0ix> 


n**0Tx) 


<  tt*0tX> 


Figure  A-5.  X-Axis  DC-10  Data  at  10  Degree  Increments 


Figure  A-2.  Y-Axis  B-52  Data  at  10  Degree 


Figure  A-l.  X-Axis  B-52  Data  at  10  Degree  Increments 


CALL  0PEN(4, INFILE, 1,IER) 

TYPE  "READING  PROCESSED  DATA 

READ (4,30)  1DIR 
30  FORMAT (A2) 

TYPE"DI R :  ",IDIR 
READ (4)  MPS 
TYPE"MPS :  ",MPS 
READ (4)  PRESX 
TYPE  "X  RES:  ", PRESX 
READ (4)  PRESY 
TYPE  "Y  RES:  ", PRESY 
READ (4)  IPX 

IF  (ITEST.NE.O)  TYPE"IPX:  ",IPX 
DO  40  1=1, IPX 

READ  (4)  XPROC(I) 

IF  (ITEST.NE.O)  TYPE  "XPROC ( " , I , " ) =  ",XPROC(I) 
IF  (XPROC  (I)  .GT.XPMAX)  XPMAX  =  XPROC ( I ) 

XPINDEX (I) =1 

40  CONTINUE 

READ (4)  IPY 

IF  (ITEST.NE.O)  TYPE"IPY:  ",IPY 
DO  50  J=1 , IPY 

READ (4)  YPROC(J) 

IF  (ITEST.NE.O)  TYPE  " YPROC ( " , J , " ) =  " ,YPROC(J) 
IF  (YPROC (J) .GT.YPMAX)  YPMAX= YPROC ( J) 

YPINDEX (J) =J 

CONTINUE 

CALL  CLOSE ( 4 , IER) 

C 

C  COMPARE  MODEL  AND  PROCESSED  DATA 

C 

DO  400  IP=1,NUMPLANES 

RXMAX=0 

RYMAX=0 

MODFILE (1) ="MD" 

MODFILE (2) =IDIR 

C  INITIALIZE  PROCESSED  DATA  TEST  ARRAYS 

DO  60  1=1, IPX 

XPTEST (I) =XPROC (I) 

XPIN ( I ) =XP INDEX ( I ) 

60  CONTINUE 

DO  70  J=1 , IPY 

YPTEST(J) =YPROC ( J) 

YPIN (J) =YPINDEX (J) 


m 

50 


111 


CONTINUE 


PICK  PLANE  MODEL 

GOTO  (110,120,130,140,150)  , IP 

MODFILE (3) = "B5 " 

MODFILE ( 4 ) = " 2 . " 

MODFILE (5) =”DB" 

MODFILE (6) =0 
GOTO  200 
MODFILE (3) ="KC" 

MODFILE (4) ="13" 

MODFILE (5) = " 5 . " 

MODFILE (6) = "DB" 

MODFILE (7) =0 
GOTO  200 
MODFILE (3) = " DC " 

MODFILE  (4) ="10" 

MODFILE (5) ="X. " 

MODFILE (6) = "DB" 

MODFILE (7) =0 
GOTO  200 
MODFILE (3) = "FI " 

MODFILE (4) = " 5 . " 

MODFILE (5) ="DB" 

MODFILE (6) =0 
GOTO  200 
MODFILE (3) ="MI" 

MODFILE (4) ="G2" 

MODFILE (5) ="1 . " 

MODFILE (6 ) ="DB" 

MODFILE (7) =0 

CALL  OPEN  (3, MODFILE, 1,IER) 

WRITE (10,201)  (MODFILE(I) ,1=1,4) 

FORMAT (IX, "JUST  OPENED  FILE:  " , 7A2) 

GET  MODEL  DATA 

READ (3)  RESOLU 

TYPE"MODEL  RESOLUTION:  ", RESOLU 
READ (3)  IMX 

IF  (ITEST.NE.0)  TYPE"IMX:  ",IMX 
DO  220  1=1, IMX 

READ (3)  XMODEL ( I ) 

IF  (ITEST.NE.0)  TYPE " XMODEL ( " , I , " ) =  ",XMODEL(I) 
IF  ( XMODEL ( I ) . GT . RXMAX )  RXMAX=XMODEL ( I ) 

XMINDEX (I) =1 

CONTINUE 
READ (3)  IMY 

IF  (ITEST.NE.0)  TYPE" IMY :  ",IMY 


DO  240  J= 1 , IMY 

READ ( 3 )  YMODEL(J) 

IF  (YMODEL(J) .GT. RYMAX)  RYMAX =YMODEL ( J) 
YMINDEX  (J) =J 

240  CONTINUE 

CALL  CLOSE ( 3 , IER) 

C  NORMALIZE  MODEL  DATA 

IF  ( I  TEST .  NE . 0 )  TYPE"NORMALI ZING  DATA.." 

IF  (ITEST.NE.O)  TYPE"RXMAX :  ",RXMAX,"  RYMAX :  " , RYMAX 
DO  250  1=1, IMX 

XMODEL ( I ) =XMODEL { I ) / RXMAX 

250  CONTINUE 

DO  26C  J= 1 , IMY 

YMODEL(J) =YMODEL ( J ) /RYMAX 

260  CONTINUE 


C  NORMALIZE  PROCESSED  DATA 


270 


272 


273 

274 

275 

280 


281 


IF  (ITEST.NE.O)  TYPE " XPMAX  = " , XPMAX , "  YPMAX= " , YPMAX 

XSCALE=FLOAT (IMX) /FLOAT (IPX) 

YSCALE=FLOAT( IMY) /FLOAT ( I PY ) 

IF  (ITEST.NE.O)  TYPE " XSCALE :  ",XSCALE,"  YSCALE : " , YSCALE 

DO  270  1=1, IPX 

XPTEST ( I ) =XPTEST ( I ) /XPMAX 
XPIN(I)=XPIN(I)*XSC.  \E 

CONTINUE 
DO  272  J= 1 , IPY 

YPTEST ( J ) =YPTEST ( J ) /YPMAX 
YPIN (J) =  YPIN (J) * YSCALE 

CONTINUE 

IF  ( (IPR.NE.l) .OR. (IP.NE.l) )  GOTO  280 
WRITE (12,900) 

WRITE (12,901)  ( INFILE (L) ,L=1,10) 

WRITE (12,902) 

WRITE (12,273) 

WRITE (12,274)  (XPTEST(L) ,L=1,IPX) 

WRITE (12,275) 

WRITE (12,274)  (YPTEST (L) ,L=1,IPY) 

FORMAT (IX, "DUMP  OF  X  PROCESSED  DATA") 

FORMAT (10F6. 2) 

FORMAT (IX, "DUMP  OF  Y  PROCESSED  DATA") 

WRITE (12,902) 

IF  (ITEST.EQ.0)  GOTO  285 
TYPE"DUMP  X  OF  PROCESSED  ARRAYS:  " 

WRITE (10,281)  (XPTEST ( K ) ,K=1,IPX) 

FORMAT (10F6. 2) 
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WRITE (10,281)  (XPIN(K) ,K=1,IPX) 

TYPE 

TYPE"DUMP  Y  OF  PROCESSED  ARRAYS:" 

WRITE (10,281)  (YPTEST(K) ,K=1,IPY) 

WRITE (10,281)  (YPIN(K) ,K=1,IPY) 

TYPE 

IF  (IPR.EQ.l)  WRITE (12,902) 

IF  (IPR.EQ.l)  WRITE (12 ,902) 

IF  (IPR.EQ.l)  WRITE (12 ,926) 

IF  (IPR.EQ.l)  WRITE (12,928)  (MODFILE (L) ,L=1 , 7) 

IF  (IPR.EQ.l)  WRITE (12 ,908)  RESOLU , RESOLU 
IF  (IPR.EQ.l)  WRITE (12,910)  IMX,IMY 
IF  (IPR.EQ.l)  WRITE (12 ,950) 

CALL  INTERP ( IPX , XPTEST , XPIN , IMX , XMODEL , XM INDEX , RDIV , I TEST , I  PR) 
RXDIV=RDIV 

IF  (IPR.EQ.l)  WRITE (12 ,952) 

CALL  INTERP ( IPY , YPTEST , YPIN , I MY , YMODEL, YM INDEX , RDI V , I TEST  ,  I  PR) 
RYDIV=RDIV 

DO  290  1=1, IPX 

XPTEST ( I ) =XPROC ( I ) /XPMAX 

CONTINUE 
DO  295  J=1 , IPY 

YPTEST ( J) =YPROC { J ) /YPMAX 

CONTINUE 

IF  (IPR.EQ.l)  WRITE (12,954) 

CALL  I NTEGRATE ( IMX , I PX , RESOLU , PRESX , XMODEL , XPTEST , RDEV , 

* ITEST , IPR) 

RIXDEV=RDEV 

IF  (IPR.EQ.l)  WRITE (12,956) 

CALL  INTEGRATE (IMY, IPY, RESOLU, PRESY, YMODEL, YPTEST, RDEV, 

* ITEST , IPR) 

RIYDEV=RDEV 

CALL  LINTEG ( IMX , IPX , RESOLU , PRESX , XMODEL , XPTEST , RDEV , RDEV 2 , 

* ITEST , I  PR) 

RI2XDEV=RDEV 

RLXDEV=RDEV2 

CALL  LINTEG ( IMY , I PY , RESOLU , PRESY , YMODEL , YPTEST , RDEV , RDEV 2 , 

* ITEST , I  PR) 

RI 2YDEV=RDEV 
RLYDEV=RDEV2 

RTOTAL= (RESOLU /PRESX ) *RXDIV+ (RESOLU/PRESY) *RYDIV 

RITOTAL=RIXDEV+RIYDEV 

RI 2TOTAL=RI 2XDEV+RI 2YDEV 

RLTOTAL=RLXDEV+RLYDEV 


GOTO  (310,320,330,340,350)  , IP 


BX=RXDIV 
BY  =  RYDI V 
BT=RTOTAL 
BIX=RIXDEV 
BIY=RI YDEV 
BIT=RITOTAL 
BI 2X=RI 2XDEV 
BI 2Y=RI 2YDEV 
BI 2T=RI 2T0TAL 
BLX=RLXDEV 
BLY=RLYDEV 
BLT=RLTOTAL 
GOTO  400 
RKX=RXDIV 
RKY=RYDIV 
RKT=RTOTAL 
RIKX=RIXDEV 
RIKY=RI YDEV 
RI KT=RITOTAL 
RI  2X=RI 2XDEV 
RI  2Y  =  RI 2 YDEV 
RI  2T  =  RI 2TOTAL 
RLX=RLXDEV 
RLY=RLYDEV 
RLT=RLTOTAL 
GOTO  400 
DX=RXD I V 
DY=RYDI V 
DT=RTOTAL 
D I  X  =  R I XDEV 
DI Y=RI YDEV 
D I T=  R I  TOTAL 
D I 2X  =  R I 2XDEV 
D I  2Y  =  R I  2 YDEV 
DI  2T  =  RI 2TOTAL 
DLX=  PLXDEV 
DLY=RLYDEV 
DLT=RLTOTAL 
GOTO  400 
FX=RXDI V 
FY  =  RYD I V 
FT=FTOTAL 
FI  X=RI XDEV 
F  IY  =  RI YLEV 
F  I  T=  R I  TOTAL 
FI  2X  =  RI 2XDEV 
F I  2Y  =  R I 2YDEV 
F  I  2T  =  RI 2TOTAL 
FLX  =  RLXDEV 
FLY=RLYDEV 
FLT=RLTOTAL 
GOTO  400 


RMX=RXDIV 
RMY=RYDI V 
RMT=RTOTAL 
RIMX=RIXDEV 
RIMY=RIYDEV 
RIMT=RITOTAL 
RI 2MX=RI 2XDEV 
RI 2MY=RI 2YDEV 
RI 2MT=R1 2T0TAL 
RLMX=RLXDEV 
RLMY=RLYDEV 
RLMT=RLTOTAL 

CONTINUE 

TYPE 

WRITE (10,500) 

WRITE (10,902) 

WRITE (10,510) 

WRITE (10,902) 

WRITE (10,520)  BX , BY , BT ,BIX,BIY,BIT 
WRITE (10,522)  RKX , RKY, RKT, RIKX , RI KY , RI KT 
WRITE (10, 524)  DX , DY , DT , DIX , DI Y , DIT 
WRITE (10,526)  FX , FY , FT , FIX , FI Y , FIT 
WRITE (10,528)  RMX , RMY , RMT , RIMX , RIMY , RIMT 
WRITE (10,902) 

WRITE (10,530) 

WRITE (10,902) 

WRITE (10,520)  BI2X,BI2Y,BI2T, BLX , BLY , BLT 
WRITE (10,522)  RI 2X , RI 2Y , RI 2T , RLX , RLY , RLT 
WRITE (10,524)  DI2X,DI2Y,DI2T, DLX , DLY , DLT 
WRITE (10,526)  FI2X,FI2Y,FI2T,FLX, FLY , FLT 
WRITE (10,528)  RI 2MX , RI 2MY , RI 2MT , RLMX , RLMY , RLMT 
TYPE 

FORMAT (IX,"  SUMMARY  OF  CORRELATION") 

FORMAT ( 10X , "POLY  INTERP" , 3 OX ," INTEGRATION" ) 

FORMAT (1X,"B52:  " , 3F10 . 4 , 5X , 3F1 0 . 4 ) 

FORMAT ( IX , "KC13  5 :  " , 3F10 . 4 , 5X , 3F1 0 . 4 ) 

FORMAT (IX, "DC10:  " , 3F10 . 4 , 5X , 3F10 . 4 ) 

FORMAT (IX, "FI  5:  " , 3F10 . 4 , 5X , 3F1 0 . 4 ) 

FORMAT (IX, "MIG21  " , 3F10 . 4 , 5X , 3F10 . 4 ) 

FORMAT (1 OX, "2ND  INTEGRATION" , 15X , "LINEAR  INTERP") 

CALL  CLOSE (3,1 ER) 

IF  (1PR.NE.1)  GOTO  999 

PRINT  RESULTS 

DO  800  1=1,24 

WRITE (12,902) 


CONTINUE 


WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 


(12,902) 

(12,902) 

(12.900) 

(12.901) 

(12.902) 
(12,904) 
(12,906) 
(12,908) 
(12,910) 
(12,902) 
(12,902) 
(12,912) 
(12,902) 
(12,914) 
(12,916) 
(12,918) 
(12,920) 
(12,922) 
(12,902) 
(12,902) 
(12,924) 
(12,902) 
(12,914) 
(12,916) 
(12,918) 
(12,920) 
(12,922) 
(12,902) 
(12,902) 
(12,958) 
(12,902) 
(12,914) 
(12,916) 
(12,918) 
(12,920) 
(12,922) 
(12,902) 
(12,902) 
(12,960) 
(12,902) 
(12,914) 
(12,916) 
(12,918) 
(12,920) 
(12,922) 


( INFILE ( I ) ,1=1,10) 

IDIR 

MPS 

PRESX , PRESY 
IPX , I PY 


BX , BY , BT 
RKX , RKY , RKT 
DX , DY , DT 
FX , FY , FT 
RMX , RMY , RMT 


BIX , BIY , BIT 
RI KX , RIKY , RIKT 
DIX , DI Y , DIT 
FIX,FIY,FIT 
RIMX , RIMY , RIMT 


BI2X,BI2Y,BI2T 
RI2X,RI2Y,RI2T 
DI2X,DI2Y,DI2T 
FI2X,FI2Y,FI2T 
RI 2MX , RI 2MY , RI 2MT 


BLX , BLY , BLT 
RLX , RLY , RLT 
DLX , DLY , DLT 
FLX ,FLY,FLT 
RLMX , RLMY , RLMT 


FORMAT (IX, "CORRELATION  DATA  FOR  FILE:  ") 

FORMAT (IX, 10A2) 

FORMAT ( " 

FORMAT (IX, "DIRECTION:  ",A2,"  DEGREES") 

FORMAT (IX, "SPEED:  ",F6.2,"  METERS/S") 

FORMAT (IX, "TARGET  RESOLUTION:  X:",F6.2,"  Y:",F6. 


910  FORMAT (IX, "DATA  POINTS:  X:",I3,"  Y:" ,13) 

912  FORMAT (IX, "DEVIATION  SUMMARY  (POLYNOMIAL  INTERPOLATION)”) 

914  FORMAT ( IX , "  B52:  ",3F8.4) 

916  FORMAT ( IX , "  KC135:  ",3F8.4) 

918  FORMAT ( IX , "  DC10:  ",3F8.4) 

920  FORMAT ( IX , "  F15:  ",3F8.4) 

922  FORMAT (IX,"  MIG21:  ",3F8.4) 

924  FORMAT (IX, "DEVIATION  SUMMARY  (INTEGRATION  METHOD)") 

926  FORMAT (IX, "MODEL  FILE:  ") 

928  FORMAT (IX, 7A2) 

950  FORMAT (IX, "DATA  FOR  INTERPOLATION  X-AXIS:"  ) 

952  FORMAT (IX, "DATA  FOR  INTERPOLATION  Y-AXIS:") 

954  FORMAT (IX, "DATA  FOR  INTEGRATION  X-AXIS:") 

956  FORMAT (IX," DATA  FOR  INTEGRATION  Y-AXIS:”) 

958  FORMAT (IX, "DEVIATION  SUMMARY  (2ND  INTEGRATION  METHOD)") 

960  FORMAT (IX, "DEVIATION  SUMMARY  (LINEAR  INTERPOLATION)") 


999  END 


C 

C  SUBROUTINE  DISPLAY  DATA 

/m  c 

“  C  CALLED  BY:  PROFILE. FR 

C 

C  FUNCTION:  DISPLAYS  LOW-RESOLUTION  SENSOR  DATA  ON  THE 

C  OCTEK  IMAGE  DISPLAY 

C 

C  LNAGUAGE :  FORTRAN  IV 

C 

£***★★**★*★****★*★*★★★*★***★*★■**★**★★★*****★****★★★★*************★*** 

SUBROUTINE  DDATA ( IAS , ISUBPIX , IFPX , IFPY , RLOW, ICT) 

INTEGER  XPOS , YPOS , IFPX , IFPY , ICT  ( 12  0 ) 

REAL  RLOW ( IAS , IAS ) ,MAX 

IX  =  38 
I  Y=20 

I SCALE = ISUBPIX 
MAX  =  0 

IF  (ISCALE.LT. 5)  ISCALE=5 
IF  ( I  SCALE .GT . 10 )  ISCALE=10 

DO  100  J= 1 , I FPY 

DO  100  1=1, IFPX 

IF  ( RLOW ( I , J ) . GT . MAX )  MAX = RLOW ( I , J) 

100  CONTINUE 
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DO  200  J= 1 , IFPY 

DO  200  1=1, I FPX 

RRES= ( RLOW { I , J ) *15/ MAX ) 

IGRAY= I FI X ( RRES ) 

XPOS=IX+ ( I -1 ) *ISCALE 
YPOS=IY+ (J-l) *1  SCALE 

CALL  PXF ILL { I CT , I  GRAY , XPOS , I SCALE , YPOS , I SCALE ) 

200  CONTINUE 

RETURN 
END 


SUBROUTINE  DECONVOLUTION 
CALLED  BY:  PROFILE. FR 

FUNCTION:  DECONVOLVES  INPUT  TARGET  VECTOR  WITH  RECTANGULAR 

OR  GAUSSIAN  WINDOW  AND  LOW  PASS  FILTER'S  VECTOR 
IN  THE  FOURIER  DOMAIN 

LANGUAGE:  FORTRAN  V 


SUBROUTINE  DECON ( XREAL , RREAL , IWIN , IPTS , I TEST , I PR , BLUR) 
REAL  XREAL (64) ,XIMAG(64) ,CONV(20) ,COEFF ( 10 , 10 ) 

REAL  RREAL (64) ,RIMAG(64) ,WREAL(64) ,WIMAG(64) 

C  CLEAR  ARRAYS Z 

DO  10  1=1,64 

XI MAG ( I ) =0 
RREAL (I) =0 
RIMAG(I) =0 
WREAL(I) =0 
WIMAG(I) =0 
IF  (I.GT.15)  GOTO  10 
CONV(I) =0 

10  CONTINUE 

INUM  =  3  2 

ITEMP= IWIN+ I PTS 
IF  ( ITEMP.GT. I TEMP)  INUM=64 
NFPNTS=INUM/ 4+1 
I SG=-1 

DO  100  1=1, INUM 
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RREAL ( I ) =XREAL ( I ) 

RIMAG ( I ) =XIMAG ( I ) 

100  CONTINUE 

TYPE 'CALCULATING  CONVOLVED  TARGET  TRANSFORM...' 
CALL  FOUREA ( XREAL , XIMAG , INUM, ISG) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

200 

202 

205 

C208 

C210 

C215 

C216 

C217 

218 

C219 

220 

222 

224 


IF  (ITEST.NE.l)  GOTO  290 
WRITE (10,200) 

WRITE (10,202)  (XREAL (I) ,1=1, INUM) 

WRITE (10,208) 

WRITE (10,202) 

TYPE 

IF  (IPR.NE.l) 

WRITE (12,216) 

WRITE (12,205) 

WRITE (12,217) 

WRITE (12,205) 

WRITE (12,218) 

WRITE (12,219) 

WRITE (12 , 200 ) 

WRITE (12,205) 

WRITE (12,208) 

WRITE (12,205) 

FORMAT ( IX ,' REAL  VALUES: 

FORMAT (1X,7F10.4) 

FORMAT (1X,6F11.4) 

FORMAT (IX, ' I MAG  VALUES: 

FORMAT (IX, 'MAGNITUDE: ' ) 

FORMAT (IX, 'PHASE:  ') 

FORMAT (IX, "INPUT  REAL 
FORMAT (IX, "INPUT  IMAG 
FORMAT  (  "  " ) 

FORMAT (IX,"  CONVOLVED  TARGET  FFT  OUTPUT") 

FORMAT (IX, "WINDOW  TRANSFORM:  (LENGTH= " , I  3 , " ) " ) 

FORMAT (IX, "TARGET-  TIME  DOMAIN  VECTOR") 

FORMAT (IX, "DECONVOLUTION  WINDOW:  (LENGTH:  ",I3, 


(XIMAG ( I ) ,1=1, INUM) 


GOTO  290 


( RREAL ( I ) ,1=1, INUM) 
(RIMAG (I) ,1=1, INUM) 


(XREAL (I) ,1=1, INUM) 

(XIMAG (I) , 1=1, INUM) 
') 


) 


VALUES: 
VALUES : 


>") 


GENERATE  DECONVOLUTION  WINDOW 


290  ILEN= IWIN 

IF  (BLUR.EQ. 0)  GOTO  498 
TYPE"USING  GAUSSIAN  WINDOW...." 

C  GAUSSIAN  WINDOW 

RMEAN= IWIN/ 2+1 

CALL  GCOEFF (BLUR , RMEAN , IWIN , COEFF ) 

IMEAN= IFIX ( RMEAN + . 5 ) 

RMAX=0 
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DO  485  1=1, IWIN 

IF (COEFF ( I , I MEAN ) . LE . RMAX )  GOTO  485 
RMAX  =COEFF ( I , IMEAN ) 

485  CONTINUE 

DO  490  I=1/IWIN 

CONV ( I ) =COEFF ( I , IMEAN) /RMAX 

490  CONTINUE 

IF  (IPR.NE.l)  GOTO  491 
WRITE (12,492) 

WRITE (12,202)  (CONV ( I ) , 1=1 , IWIN) 

492  FORMAT (IX, "GAUSSIAN  FUNCTION:  ") 


CONVOLVE  WINDOW 


491 


495 

496 


c 

498 


500 


C 

C 


ILEN=2*IWIN-1 
DO  495  1=1,1 LEN 

DO  495  J=1 , IWIN 

IF  ( ( I-J+l) .GT. IWIN)  GOTO  495 
IF  ( (I-J+l) .LT.l)  GOTO  495 
WREAL(I) =WREAL ( I ) +CONV ( I-J+l) 

IF  (IPR.EQ.l)  WRITE (12,496)  I,WREAL(I) 

CONTINUE 

FORMAT ( IX , "WREAL ( " , 1 2 , " ) =  "  ,  F8 . 3 ) 

GOTO  500 


RECTANGULAR  WINDOW 


DO  500  1=1, IWIN 

WREAL ( I ) = 1 


CONTINUE 

WRITE (10,224)  I LEN 

WRITE (10,202)  (WREAL ( I ) ,I=1,ILEN) 


ISG  — 1 

CALL  FOUREA (WREAL, WIMAG, INUM, ISG) 
IF  (IPR.NE.l)  GOTO  525 
WRITE (12,218) 

WRITE (12,220) 

WRITE (12,200) 

WRITE (12 ,205) 

WRITE (12,208) 

WRITE (12,205) 

WRITE (12,218) 


IWIN 

(WREAL ( I ) ,1=1, INUM) 
(WIMAG (I) ,1=1, INUM) 


DIVIDE  TRANSFORMS  TO  DECONVOLUTE  TARGET 


525  IF  (NFPNTS.EQ.O)  GOTO  550 

ICENTER= INUM/ 2+1 
ILO= ICENTER-NFPNTS/ 2 
IHI = ICENTER+NFPNTS/ 2 
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IF  (IPR.EQ.l)  WRITE (12,545)  NFPNTS , I CENTER, IHI  , ILO 
545  FORMAT (IX, "FILTERING:  "  ,  1 3  ,  "  PNTS ,  CENTER:  ",I3,"  HI:  ", 

*13,"  LO:  " , 1 3 ) 

550  DO  600  1=1, INUM 

RDI V= (WREAL(I) *WREAL ( I ) +WIMAG ( I ) *WIMAG  (I) ) 

IF  (ABS (RDIV) .GT. 0.00001)  GOTO  575 
RREAL (I ) =0 
RIMAG (I) =0 
GOTO  585 

575  RREAL (I ) = ( XREAL ( I ) *WREAL ( I ) +XIMAG (I ) *WIMAG (I) ) /RDIV 

RIMAG (I) = (WREAL(I) *XIMAG (I) -XREAL (I) *WIMAG(I) ) /RDIV 
585  IF  ( (I.GE.ILO) .AND. (I.LE.IHI) )  RREAL(I)=0 

IF  ( (I.GE.ILO) .AND. (I.LE.IHI) )  RIMAG(I)=0 

600  CONTINUE 

C  CHANGE  BACK  TO  TIME  DOMAIN 

IF  (ITEST.NE.l)  GOTO  611 
C  TYPE"TARGET  DECONVOLUTED  F-DOMAIN" 

C  TYPE" REAL:" 

C  WRITE (10,202)  (RREAL ( I ), 1=1 , INUM) 

C  TYPE" IMAG : " 

C  WRITE (10 , 202)  (RIMAG ( I ), 1=1 , INUM) 

C  TYPE 

IF  (IPR.NE.l)  GOTO  611 
WRITE (12,218) 

®  WRITE (12,205)  (RREAL ( I ) , I =1 , INUM) 

C  WRITE (12,205)  (RIMAG ( I ), 1=1 , INUM) 

WRITE (12,218) 

611  I SG=1 

CALL  FOUREA (RREAL, RIMAG, INUM, ISG) 

C  IF  (ITEST.NE.l)  GOTO  999 

C  TYPE 

C  WRITE (10,222) 

C  WRITE (10,200) 

C  WRITE (10,202)  (RREAL ( I ) , I =1 , INUM) 

C  WRITE (10,208) 

C  WRITE (10,202)  ( RIMAG ( I )  , I  =  1 ,  INUM ) 

IF  (IPR.NE.l)  GOTO  999 
WRITE (12 ,218) 

WRITE (12,222) 

WRITE (12 ,200) 

WRITE (12,205)  (RREAL (I) ,1=1, INUM) 

WRITE (12,208) 

WRITE ( 12 , 205)  (RIMAG ( I ) , 1=1 , INUM) 

999  RETURN 

END 
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c******************************************************************** 

c 

C  PROGRAM:  F0RT4T05  -  Utility  LANGUAGE:  FORTRAN  IV 

C 

C  FUNCTION:  CONVERTS  TARGET  FILES  CREATED  BY  FORTRAN  IV  PROGRAM 

C  TO  FORTRAN  V  FORMAT  COMPATABLE  WITH  THE  ECLIPSE 

C  VERSION  OF  SCENE. FR 

C 

DIMENSION  INF I LE (13) ,OUTFILE(13) 

TYPE"ENTER  INPUT  FORTRAN  4  DATA  FILE:  " 

READ (11,10)  INF ILE { 1 ) 

10  FORMAT (S13) 

TYPE"ENTER  OUTPUT  FORTRAN  FILENAME: 

READ (11,10)  OUTFILE  ( 1 ) 

CALL  OPEN (1 , INFILE, 1 , IER) 

IF  (IER.NE.l)  TYPE"ERROR  OPENING  INPUT  DATA  FILE” 

IF  (IER.NE.l)  STOP 

CALL  OPEN(2, OUTFILE, 3, IER) 

IF  (IER.NE.l)  TYPE "ERROR  OPENING  OUTPUT  DATA  FILE" 

IF  (IER.NE.l)  STOP 
READ ( 1 )  NPTS 
READ ( 1 )  FDIR 
READ ( 1 )  TRES 
READ ( 1 )  IXL 
READ ( 1 )  IYL 

TYPE"NUMBER  OF  PLANE  DATA  POINTS:  ",NPTS 
TYPE "DIRECT I ON  TARGET  IS  FACING  :  " , FDIR 
TYPE" RESOLUTION  OF  TARGET  PIXEL  :  " , TRES 
TYPE"TARGET  DEFINED  IN  WINDOX  X  :  ",IXL,"  Y:  " , IYL 
WRITE (2,101)  NPTS 
WRITE (2,102)  FDIR 
WRITE (2,102)  TRES 
WRITE (2,101)  IXL 
WRITE (2,101)  IYL 
DO  100  1=1, NPTS 

READ ( 1 )  ITARX 
READ ( 1 )  ITARY 
READ ( 1 )  ITARI 
WRITE (2,101)  ITARX 
WRITE (2,101)  ITARY 
WRITE (2,101)  ITARI 
TYPEWRITING  POINT:  ",I 

CONTINUE 

101  FORMAT ( IX , I  5 ) 

102  FORMAT ( 1 X , F 1 2 . 6 ) 

END 
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SUBROUTINE  FOUREA ( XREAL , X I MAG ,  N  , I SN ) 

DIMENSION  XREAL (N) ,XIMAG (N) 

C  FOUREA  -  VER.  1.0  -  15  JUL  84 

C  COMPUTES  1-D  FFT  USING  RADIX  2  DECIMATION  IN  TIME  ALGORITHM 

C  DATA  AND  RESULTS  ARE  PASSED  IN  REAL  ARRAYS  XREAL  AND  XI MAG 

C  ADAPTED  BY  D.  KING  FROM  IEEE  PROGRAMS  FOR  DIGITAL  SIGNAL 

C  PROCESSING,  SECTION  1.1  (P  1.1-5),  BY  C.M.  RADER. 

C  MODIFICATIONS  INCLUDE  CLEANUP  OF  CODE,  CHECKING  OF  INPUT 

C  PARAMETERS,  AND  USE  OF  REAL  RATHER  THAN  COMPLEX  ARITHMETIC. 

C  CHECK  ISN  ( 1 = INVERSE  FFT,  -l=FORWARD  FFT) 

ISI=-1 

IF  (ISN.EQ.l)  ISI=1 

C  CHECK  FOR  POWER  OF  2  UP  TO  15 

IF  (K.LT.2)  GO  TO  991 
NN  =  1 

DO  10  1=1,15 
NN=NN*2 

IF  (NN.EQ.N)  GO  TO  20 

10  CONTINUE 

GO  TO  992 
20  CONTINUE 

PI =  4 . *ATAN  ( 1 .  ) 

PI = FLOAT ( ISI ) *PI 
FN  =  N 

C  THIS  SECTION  BIT-REVERSES  THE  INPUT  DATA 

J=1 

DO  80  I  =  1 , N 
IF  (I-J)  30,40,40 
30  TEMPR= XREAL (J) 

TEMPI =X IMAG  ( J) 

XREAL (J) =XREAL  ( I ) 

XI MAG ( J) =X IMAG  ( I ) 

XREAL ( I ) =TEMPR 
X IMAG ( I ) =TEMP I 
40  M=N/2 

50  IF  (J-M)  70,70,60 

6C  J= J-M 

M=  ( K+ 1 ) / 2 
GO  TO  50 
70  J= J+M 

80  CONTINUE 

C  COMPUTE  BUTTERFLIES 
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c********************************************************************* 

c 

c 

C  SUBPROGRAM  LINTEG  LANGUID:  FORTRAN  V 

C 

C  CALLED  BY:  COMPARE. FR 

C 

C  FUNCTION:  COMPARES  EXTRACTED  TARGET  DATA  WITH  MODEL  DATA. 

C  PROGRAM  FIRST  INTEGRATES  MODEL  DATA  TO  MATCH 

C  EXTRACTED  DATA,  THEN  LINEAR  INTERPOLATES  TO 

C  CONPUTE  DEVIATION. 

C 

C  NOTE:  MAIN  PROGRAM  FOR  TARGET  VECTOR  CLASSIFICATION 

C 

c******************************************************************** 

SUBROUTINE  LINTEG (IMPTS,IPPTS, RMRES , RPRES , MODEL , PROC , RDEV , RTOTAL , IT , I P ) 

INTEGER  IMPTS  ; NUMBER  OF  MODEL  POINTS 

INTEGER  IPPTS  ; NUMBER  OF  POINTS  IN  TEST  CASE 

REAL  RMRES, RPRES , MODEL ( I MPTS ) ,PROC(IPPTS) 

REAL  RDEV, RDEV2 ,MTEST (100) ,PTEST(75) 

INTEGER  PNXTPNT , PCURPNT , MCURPNT 

IF  (IPPTS.GT.3)  GOTO  5 
RDEV=0 
RTOTAL=0 
GOTO  999 

C  GET  RESOLUTION  WIDTH 

5  INTWD= IF  IX ( RPRES / RMRES + . 5 ) 

IF  (IT.EQ.3)  TYPE" INTEGRATION  WINDOW:  " , INTWD 

C  INTEGRATE  MODEL  TO  MATCH  TEST  CASE 

RWAG=  0 
RDEV  =  0 
ICNT= 1 

10  CONTINUE 

RPT  =  0 

DO  20  K=i, INTWD 

IND  =  ( ( ICNT-1) *INTWD+K) 

IF  (IND. GT. IMPTS)  GOTO  20 
RPT=RPT+MODEL (IND) 

IF  (IT.EQ.3)  TYPE" ADDING  MODEL  POINT:  " , IND 

20  CONTINUE 

IF  (IT.EQ.3)  TYPE"MODEL  POINT  (", ICNT ,") =  ",RPT 
IF  (RPT.GT.RWAG)  RWAG=RPT 
MTEST (ICNT) =RPT 
ICNT= ICNT+ 1 
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OPEN  2 , INFILE, ERR=70 
WRITE (10,65)  ( INFILE ( J ) , J= 1 , 6 ) 

65  FORMAT (IX, "OPENING  FILE:  ",6A2) 

GOTO  90 

C  IF  ERROR  OPENING  FILE  THEN  HERE 

70  WRITE (10,75)  ( INFILE (J ), J=1 , 6 ) 

75  FORMAT (IX, "ERROR  OPENING  FILE:  " ,6A2) 

GOTO  200 

C  GET  FILE  HEADER  INFORMATION 

90  READ (2, 81)  RESOLD 

TYPE"FILE  RESOLUTION:  “ , RESOLU 
READ (2,82)  IXPNTS 
TYPE "LENGTH  OF  X-AXIS:  ", IXPNTS 
DO  95  J=l, IXPNTS 

READ  (2,82)  XDATA(J) 

95  CONTINUE 

READ (2, 82)  IYPNTS 

TYPE" LENGTH  OF  Y-AXIS:  ", IYPNTS 

DO  100  J  =  l, IYPNTS 

READ (2, 82)  YDATA(J) 

100  CONTINUE 

81  FORMAT (F12. 6) 

82  FORMAT (15) 

IF  ( I  AXIS . EQ. 2 )  GOTO  120 

110  CALL  PLOT10 (XDATA, IXPNTS, I, X0,Y0, RESOLU) 

GOTO  130 

120  CALL  PLOT10  (YDATA, IYPNTS, I ,X0,Y0 , RESOLU) 

130  CONTINUE 

CLOSE  2 

200  CONTINUE 

999  END 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
15 


20 

» 

30 


40 


50 


55 


GOTO  (1,2, 3, 4, 5, 6,7, 8, 9, 10), I 
INFILE  (2)  =  "00" 

GOTO  15 
INFILE (2) ="10” 

GOTO  15 
INFILE (2) ="20" 

GOTO  15 
INFILE (2) ="30" 

GOTO  15 
INFILE (2) ="40" 

GOTO  15 
INFILE (2) ="50" 

GOTO  15 
INFILE (2) ="60" 

GOTO  15 
INFILE (2) ="70" 

GOTO  15 

INFILE (2) ="80" 

GOTO  15 
INFILE (2) ="90" 

INFILE (1) ="MD" 

GOTO  (20,30,40,50,55) ,1 PLANE 

INFILE ( 3 ) = "B5 " 

INFILE  (4) ="2." 

INFILE (5) ="DF" 

INFILE (6) =0 
GOTO  60 
INFILE (3) ="KC" 
INFILE(4) ="13" 

INF ILE ( 5 ) = " 5 . " 

INFILE ( 6 ) = "DF" 

INFILE (7) =0 
GOTO  60 
INFILE (3) ="DC" 

INFILE  (4) ="10" 

INFILE ( 5 ) = "X . " 

INFILE ( 6 ) ="DF" 

INFILE (7) =0 
GOTO  60 
INFILE ( 3 ) = " F 1 " 

INFILE (4) ="5. " 
INFILE (5) ="DF" 
INFILE (6) =0 
GOTO  60 
INFILE (3 ) =" MI " 
INFILE ( 4 ) = "G2 " 
INFILE (5) ="1. " 
INFILE (6) ="DF" 
INFILE (7) =0 
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WRITE (12,510)  (RDATA(I) ,I=1,INUM) 


999 


RETURN 

END 


m 


FUNCTION:  PLOTS  TARGET  DATABASE  AT  LINE  PRINTER. 

DATABASE  CREATED  BY  PROGRAM  TGEN.FR.  TARGET  DATA 
HAS  A  -.DF  EXTENTI ON . 


MACRO: 


FORTRAN  INVPLOT 

RLDR  INVPLOT  PLOTIO  PLOT5.LB  @FLIB@ 


£********************************************************************* 

C 

C  PLANE  INVENTORY  PLOTTER  -  Utility  LANGUAGE:  FORTRAN  V 

C 

c 
c 
c 
c 
c 
c 
c 

£**★**★********’******★*★****★*****★★**★**★★*★**■*★**★★***★**★*******★* 

PARAMETER  ASIZE=50 

REAL  XDATA (100) ,YDATA (100) 

DIMENSION  INFILE (13) 

COMPILER  FREE 

C  SELECT  DATA  BASE 

TYPE 
TYPE" 

TYPE 


TARGET  DATA  BASE  PLOTTER" 


TYPE" 

1) 

B-52  " 

TYPE" 

2) 

KC-135" 

TYPE" 

3) 

DC-10" 

TYPE" 

4) 

F-15" 

TYPE" 

5) 

MIG-21" 

TYPE 

TYPE" 

-1) 

ABORT  " 

TYPE 

ACCEPT "SELECT 

ONE 

->  ",  I PLANE 

IF  (I  PLANE . EQ . 

-1) 

GOTO  999 

TYPE 

TYPE" 

AXIS  TO  PLOT: 

TYPE 

TYPE" 

1) 

X-AXIS  " 

TYPE" 

2) 

Y-AXIS  " 

TYPE 

ACCEPT "SELECT  AXIS  ->  ",IAXIS 
TYPE 

DO  200  1=1,10 
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9  CONTINUE 

F  (NP1MK  +  1) =FLAST 

10  CONTINUE 

C  CALCULATE  PN(X)  FOR  VARIOUS  VALUES  OF  X 

20  RDI V=0 

RLDEV=0 

I VAR= I F IX (FLOAT { INUM) / 4 . 0 ) 

IBNDRY  =  INUM- I VAR 

ILBND= I VAR+ 1 

DO  30  J=ILBND, IBNDRY 

PNOFX=0 

DO  29  1=1,  NP1 

c  IF  (IT.EQ.l)  TYPE'F ( ' ,1 , ' ) = ' ,F (I ) , '  Y(I)=',Y(I) 

PNOFX  =  F ( I )  +  ( J-Y  ( I ) ) *PNOFX 

29  CONTINUE 

IF  (IT.EQ.l)  WRITE (10,600)  J,PNOFX 
600  FORMAT ( ’NUMBER:  ’,14,’  F(X):  ’,F10.2) 

IF  (IT.EQ.l)  TYPE"MODEL  DATA ( " , J , " ) =  " , RDATA ( J) 


GUESS (J) =PNOFX 

C  CHECK  TO  SEE  IF  MID  POINT  (J-1)  IS  RADICAL 

C  IF  YES,  SUBTRACT  ITS  DEVIATION  FROM  TOTAL 

C 

ILOW= J-2 

IF  (ILOW.LT.l)  ILOW= 1 
IMID= J-l 

IF  ( IMID.LT. 1 )  IMID=1 
AVG= (GUESS (J) -GUESS (ILOW) ) /2 
IF  (ABS (GUESS (IMID) -AVG) .LT.0. 8)  GOTO  610 

RDEV=RDEV-RLDEV  , 

IF  (IT.EQ.l)  TYPE "DELETING  RDEV  OF  POINT:  " ,IMID  • 

6 ! 0  IF  ( (PNOFX.GT.l)  .OR.  (PNOFX.LT. 0) )  GOTO  30 

RTEMP=ABS (PNOFX- RDATA  ( J) ) *ABS (PNOFX-RDATA ( J) ) 

RDIV  =  RDIV+RTEMP  ,! 

RLDI V=RTEMP  j 

30  CONTINUE 

IF  (IT.EQ.l)  TYPE 'DEVIATION  FOR  PROCESSED  RDATA:  ,RDIV  , 

IF  (IPR.NE.l)  GOTO  999 

WRITE (12,505) 

505  FORMAT (IX, "INTERPOLATED  DATA") 

WRITE (12,510)  (GUESS  (I)  , 1  =  1, INUM) 

510  FORMAT (10F6. 2)  ; 

WRITE (12,515)  : 

515  FORMAT (IX, "MODEL  DATA") 
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303 

FORMAT ( IX , " INTEGRATION 

DATA 

(MODEL):  " ) 

304 

FORMAT ( IX , " INTEGRATION 

DATA 

(TEST) : ") 

305 

FORMAT (10F6. 2) 

999 

RETURN 

END 

C 

C  SUBROUTINE  INTERP:  MODIFIED  FROM  ' ELEMENTART  NUMERICAL  ANALYSIS' ,BY 
C  S.D  CONTE  AND  CARL  DE  BOOR  (1972) ,  PG  227. 

C 

C  CALCULATION  OF  THE  NEWTON  FORM  OF  THE  POLYNOMIAL  OF  DEGREE 

C  LESS  THAN  N,  WHICH  INTERPOLATES  AT  F (X)  AT  Y(I),  1=1,..., NPI. 

SUBROUTINE  INTERP ( IMNUM , F , Y , INUM , RDATA , DINDEX , RDIV , IT , I PR) 

INTEGER  INUM, NPI , IPNUM 

REAL  F ( IMNUM) ,Y(IMNUM) , PNOFX , X , DX , RDATA ( INUM) , DINDEX ( INUM) 
REAL  GUESS (100) , RLAST , RNEXT , RDEV 
X  =  0 
DX  =  0 

IF  (INUM.GT.3)  GOTO  5 
RDEV=0 
GOTO  999 

5  NPI = IMNUM 

C  CALCULATE  THE  DIVIDED  DIFFERENCES 

N=NP1-1 

IF  (N.EQ.O)  GOTO  20 

DO  10  K=1,N 

FLOATK=K 
NP1MK=NP1-K 
FLAST=F (1) 

DO  9  1=1 ,NP1MK 

IPK=I+K 

DY  =  Y (IPK) -Y  (I ) 

IF  (DY.EQ.O)  GOTO  8 

F ( I ) = ( F ( I + 1 ) -FLAST ) /DY 
FLAST=F (1+1) 

GOTO  9 

8  F ( I ) =F ( I +1 ) /FLOATK 
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GOTO  10 


IF  (ABS (RMAT-FPT) .LE.ABS (RRMAT-FPT) )  GOTO  30 
ICNT=ICNT-1 


INTEGRATE  MODEL  DATA  AND  NORMALIZE 


RPT=  0 

DO  40  1=1,  ICNT 

RPT=RPT+MODEL ( I ) 

CONTINUE 
MTEST (1) =RPT 
RMAX=RPT 

DO  50  1=2, IPPTS 
RPT  =  0 

DO  4  5  K  =  1 , INTWD 

INDEX=ICNT+K+ (1-2) *INTWD 
IF  (INDEX. GT.IMPTS)  GOTO  45 
RPT=RPT+MODEL ( INDEX ) 

CONTINUE 

IF  (RPT.GT.RMAX)  RMAX=RPT 
MTEST ( I ) =RPT 
IF  (IT.NE.2)  GOTO  50 
TYPE 

TYPE"  CALC  POINT:  ",I 

TYPE "ENERGY:  MODEL: " ,RPT , "  PROC : " , PROC ( I ) 

TYPE" RMAX :  " , RMAX 

CONTINUE 

NORMALIZE  MODEL  DATA 

DO  55  1=1, IPPTS 

MTEST ( I ) =MTEST ( I ) /RMAX 

CONTINUE 

CALCULATE  DEVIATION 
RDEV=0 

DO  200  1=1, IPPTS 

RDEV=RDEV+ ABS (MTEST ( I ) -PROC  ( I ) ) 

IF  (IT.EQ.2)  TYPE "DEVI AT I ON  FOR  POINT:", I," 

CONTINUE 

IF  (IP.NE.l)  GOTO  999 
WRITE (12,303) 

WRITE (12,305)  (MTEST (I) ,1=1, IPPTS) 

WRITE (12,304) 

WRITE (12 ,305)  (PROC(I) ,1=1, IPPTS) 


SUROUTINE  INTEGRATE 


VERSION  1.0 


CALLED  BY:  COMPARE. FR 

FUNCTION:  COMPARES  EXTRACTED  TARGET  VECTOR  WITH  DATABASE 

VECTOR 

NOTE:  THIS  PROGRAM  INCLUDED  FOR  COMPARISON  PURPOSES  ONLY 
NOT  A  RELIABLE  COMPARISON  TECHNIQUE 


SUBROUTINE  INTEGRATE ( IMPTS , IPPTS , RMRES , RPRES , MODEL, PROC , RDEV, IT , IP) 
INTEGER  IMPTS, IPPTS 

REAL  RMRES, RPRES, MODEL (IMPTS) ,PROC(IPPTS) , RDEV, MTEST (100) 

IF  (IPPTS. GT. 3)  GOTO  2 
RDEV  = 0 
GOTO  999 

GET  INTEGRATION  WIDTH 
INTWD=IFIX (RPRES/RMRES) 

IF  (IT.EQ.2)  TYPE" INTEGRATION  WINDOW:  " , INTWD 

INTEGRATE  MODEL  TO  WAG  AT  NORMALIZATION 
RWAG=0 

DO  5  1=1, IPPTS 
RPT  =  0 

DO  4  K=l, INTWD 

IND= ( I -1 ) * INTWD +K 
IF  (IND.GT. IMPTS)  GOTO  4 
RPT=RPT+MODEL ( IND) 

CONTINUE 

IF  (RPT.GT.RWAG)  RWAG=RPT 

CONTINUE 

IF  (IT.EQ.2)  TYPE "MAX  ELEMENT= " , RWAG , "  (WAG)" 

START  AT  LOW  END  AND  MATCH  ENERGY 

AS  CLOSE  AS  POSSIBLE  TO  FIRST  POINT  OF  PROCESS  DATA 

FPT=PROC (1) 

RMAT=0 

ICNT=0 

IF  (RMAT.GE.FPT)  GOTO  20 
RRMAT=RMAT 
ICNT= ICNT+ 1 

RMAT=  RMAT+ (MODEL (ICNT) /RWAG) 

IF  (IT.EQ.2)  TYPE"NEW  NORMALIZED  MODEL  PNT ( " , ICNT , " ) : 


REAL  PIX,  RLOW ( SAR , SAR) 
COMPILER  FREE 
I ERROR =0 

TYPE’ READING  DATA _ ' 

DO  250  1=1,  NBLKX 

DO  250  J=l,  NBLKY 


READ  (2)  PIX 
RLOW (J, I) =PIX 
TYPE  ’PIXEL' ,J,I. '=' .PIX 


CONTINUE 
GOTO  999 

RETURN 


SUBROUTINE  GET  HEADER  - (GETHDR) 


CALLED  BY:  PROFILE. FR 


LANGUAGE:  FORTRAN  V 


FUNCTION:  OPENS  TARGET  DATA  FILE 

READS  IN  FILE  DESCRIPTION  HEADER 


SUBROUTINE  GETHDR ( SAMP , DIR , MPS , FRATE , NBLKX , NBLKY , LOWRES ) 
INTEGER  SAMP, NBLKX  , NBLKY, FRATE 
REAL  LOWRES, DIR, MPS 
DIMENSION  INFILE  (20) 

TYPE’ ENTER  NAME  OF  LOWRES  TARGET  FILE:  ' 

TYPE 

READ (11,5)  INFILE ( 1 ) 

FORMAT (S20) 

CALL  OPEN (2 , INFILE) 

READ (2)  SAMP 

TYPE  'TARGET  FILE  HAS:  ', SAMP ,' FRAMES ' 

READ  (2)  DIR 

TYPE  'DIRECTION;  ',DIR,'  RAD' 

READ  (2)  MPS 

TYPE  'VELOCITY:  ',MPS,'  M/S' 

READ  (2)  FRATE 

TYPE  'FRAME  RATE:  ', FRATE,'  HZ' 

READ  (2)  NBLKX,  NBLKY 

TYPE  'WINDOW  SIZE:  ', NBLKX,’  BY  ', NBLKY 

READ  (2)  LOWRES 

TYPE  'DATA  RESOLUTION:  ', LOWRES,'  METERS  ' 

RETURN 

END 


I 


11./'. 


Q*  *******************************************************************  * 

c 

C  SUBROUTINEGaussianCOEFFicient  generator  LANGUAGE:  FORTRAN  V 
C 

C  CALLED  BY:  PROFILE . FR ,  SCENE. FR 

C 

C  FUNCTION:  RETURNS  A  COEFFICENT  ARRAY  CONTAINING  THE 

C  DISCRETE  VALUES  OF  A  2-D  GAUSSIAN  SPREAD  FUNCTION 

C 

C  VARIABLES:  SIG  -  STANDARD  DEVATION 

C  {EQUIVILENT  TO  (DETSIZE) **0 . 5 ) 

C  RM  MEAN 

C  NUM  -  WIDTH 

C  COEFF  -  RETURN  DATA 

C 


TYPE"THE  FILE  HAS",INUM,"  FRAMES" 

ACCEPT "HOW  MANY  FRAMES  TO  PRINT?  " , IN 

WRITE (12,400)  ( INFILE ( I ) ,1=1,5) 

WRITE (12,405) 

WRITE (12,410)  INUM 
WRITE (12,430)  ISIZE 
WRITE (12,445)  ISP 
WRITE (12,425)  IFRATE 
WRITE (12,440)  FPRINT 
WRITE (12,450)  SIG 
WRITE (12,455)  RDRIFT 
WRITE (12,405) 

WRITE (12,435)  ILXST,ILYST 
WRITE (12,420)  RMPS 
WRITE (12,415)  RD1R 
WRITE (12,470)  TRES 
WRITE (12,460)  ITGRAY 
WRITE (12,465)  PIXDIST 
WRITE (12,405) 

WRITE (12,405) 

400  FORMAT (IX,"  FRAME  DUMP  FOR  FILE:  ",5A2) 

405  FORMAT (80X) 

410  FORMAT (IX, "FILE  CONTAINS:  ",I12,"  FRAMES") 

415  FORMAT (IX, "TARGET  DIRECTION:  ",F12„6,"  RADIANS") 

420  FORMAT (IX, "TARGET  SPEED:  ",F12.6,"  M/S") 

425  FORMAT (IX, "FRAME  RATE:  ",I12,"  FRAMES/S") 

£>  430  FORMAT (IX, "FOCAL  PLANE  SIZE:  "  ,112,"  SQ  DETECTORS") 

435  FORMAT (IX, "TARGET  POSITION:  ",14,14,"  @  FRAME  1") 

440  FORMAT (IX, "FOOTPRINT:  ",F12.6,"  METERS") 

445  FORMAT (IX, "SIZE  OF  SUBPIXEL:  ",I12,"  SQ  " ) 

450  FORMAT (IX, "STD  DEV  OF  OTF :  ",F12.6) 

455  FORMAT (IX, "DRIFT  RATE:  ",F12.6,"  METERS/SEC") 

460  FORMAT (IX, "TARGET  THRESHOLD:  ",I12,"  GRAYLEVEL") 

465  FORMAT (IX, "PIXEL  VELOCITY:  ",F12.6,"  PIXELS/SEC" ) 

470  FORMAT (IX, "TARGET  RESOLUTION:  ",F12.6,"  METERS") 

DO  100  1=1, IN 

CALL  GETDA (IFPX, ISIZE, ISIZE, RLOW, I ERROR) 

WRITE (12,200)  I 
WRITE (12,205) 

DO  50  J=l, ISIZE 

WRITE (12,210)  ( RLOW ( K , J ) ,K=1, ISIZE) 


50 

CONTINUE 

WRITE (12,205) 

WRITE (12,205) 

100 

CONTINUE 

200 

FORMAT (IX, "DATA  FRAME 

205 

FORMAT (IX," 

210 

FORMAT ( IX, 20F5.1) 

END 

128 
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MFRAME=HFRAME 
GOTO  300 

200  IF  ( (1-HFRACT) .EQ.LFRACT)  GOTO  250 

MFRAME=LFRAME 
GOTO  300 

250  IF  (HFRAME . GT . LFRAME )  GOTO  275 

MFRAME=HFRAME 
GOTO  300 

275  MFRAME=LFRAME 

300  MPIX=IFIX ( FRACT*MFRAME+ . 5) 

AFRAC= FLOAT (MPIX) /FLOAT (MFRAME) 

IF  (INEG.EQ.l)  MPIX=MPIX* (-1) 
ERROR= ABS ( AFRAC-FRACT) 

RETURN 

END 


C 

C 

C 

C 

C 

C 

C 

C 


UTILITY  FRAME  DUMP  LANGUAGE:  FORTRAN  IV 

FUNCTION:  Dumps  Sensor  Data  Created  By  PROGRAM  SCENE. FR 

******************************************************************** 


PARAMETER  IFPX=20 , IFPY=20 

REAL  RLOW ( IFPX , IFPY) , RDIR, RMPS , FPRINT 
DIMENSION  INFILE (13) 

ACCEPT " ENTER  SCENE  DATA  FILENAME  ->  ", 
READ (11,11)  INFILE ( 1 ) 

11  FORMAT (S13) 

CALL  OPEN ( 2 , INFILE , 1 , IER) 

READ (2)  INUM 
READ ( 2 )  ISIZE 
READ (2)  ISP 
READ (2)  IFRATE 
READ (2)  SIG 
READ (2)  RDRIFT 
READ (2)  RLOWRES 
READ (2)  ILXST 
READ (2)  ILYST 
READ (2)  RMPS 
READ (2)  RDIR 
READ (2)  TRES 
READ (2)  ITGRAY 
READ ( 2 )  PIXDIST 
TYPE 
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q* ******************************************************************* 

c 

C  SUBROUTINE  FRAC  LANGUAGE:  FORTRAN  V 

C 

C  CALLED  BY:  PROFILE. FR 

C 

C  FUNCTION:  CONVERTS  THE  REMAINDER  OF  THE  SAMPLING  RATE 

C  ( #DETECTORS/FRAME)  TO  FRACTION.  THE  DENOMINATOR 

C  IS  USED  TO  DETERMINE  TARGET  PHASING 

C 

Q* ******************************************************************* 


SUBROUTINE  FRAC ( FRACT , MFRAME , MPIX , ERROR) 

REAL  FRACT , ERROR , FRACNM , REALNM , HFRACT , LFRACT 
REAL  NFRAC , NHFRAC 

INTEGER  MFRAME, MPIX, INTNUM , MAXNUM 
INTEGER  LFRAME , HFRAME 

C  INITIALIZE  VARIABLES 


INEG=0 

IF  (FRACT. GE.O)  GOTO  20 
FRACT=ABS (FRACT) 

INEG=1 

20  HFRACT=0 

LFRACT=1 
MAXNUM=5 

DO  100  1=1, MAXNUM 

C  GET  FRACTIONAL  PART  OF  NUMBER 

REALNM=FRACT* I 
INTNM=IFIX (REALNM) 

FRACNM =REALNM-INTNM 
C  CHECK  MINIMUM  LOW  FRACTION 

IF  ( FRACNM. GE. LFRACT)  GOTO  120 
LFRACT = FRACNM 
LFRAME= I 

C  CHECK  MINIMUM  LOW  FRACTION 

120  NFRAC=ABS ( 1-FRACNM) 

NHFRAC=ABS (1-HFRACT) 

IF  (NFRAC. GE. NHFRAC)  GOTO  100 
HFRACT = FRACNM 
HFRAME= I 

100  CONTINUE 

C  DETERMINE  SMALLER  FRACTION 

IF  ( (1-HFRACT) .GE. LFRACT)  GOTO  200 
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MMAX=1 

90  IF  (MKAX-N)  100,130,130 

100  ISTEP=2*MMAX 

DTHETA=PI /FLOAT (MMAX) 

DO  120  M= 1 , MMAX 
THETA=DTHETA*FLOAT (M-l) 

C=COS (THETA) 

S=SIN (THETA) 

DO  110  I =M ,N , ISTEP 
J= I +MMAX 

TEMPR=C*XREAL ( J) -S*XIMAG ( J) 

TEMPI =C*XIMAG (J) +S*XREAL( J) 

XREAL(J) =XREAL(I) -TEMPR 
XI MAG (J) =XIMAG (I) -TEMPI 
XREAL ( I ) =XREAL ( I ) +TEMPR 
XI MAG ( I ) =XIMAG ( I ) + TEMP I 
110  CONTINUE 

120  CONTINUE 

MMAX= I  STEP 
GO  TO  90 

130  IF  (ISI)  160,140,140 

C  FOR  INVERSE  TRANSFORM  (ISI=1)  -  MULT  BY  1/N 

140  DO  150  1=1, N 

XREAL ( I ) =XREAL ( I ) / FN 
<9u  X I  MAG  ( I )  =  X I  MAG  ( I )  /  FN 

150  CONTINUE 

160  RETURN 

C  ERRORS 

991  TYPE  "<7>N  =  ",N 

STOP  FOUREA01  -  N  MUST  BE  GREATER  THAN  1 

992  TYPE  " < 7 >N  =  ",N 

STOP  FOUREA02  -  N  MUST  BE  INTEGER  POWER  OF  2 


END 


IF  (IND.LT. IMPTS)  GOTO  10 

IF  (IT.EQ.3)  TYPE "MAX  MODEL  POINT:  " , RWAG 

1CNT= ICNT-1 
DO  50  1=1/ ICNT 

MTEST ( I ) =MTEST ( I ) / RWAG 

50  CONTINUE 

C  TEST  DEVIATION  FROM  NORMALIZED  MODEL 

RDEV=0 
1  =  1 

60  IF  ((I.GT.ICNT) .AND. (I.GT.IPPTS) )  GOTO  100 

IF  (I.GT.ICNT)  GOTO  65 
RMPNT=MTEST (I) 

GOTO  70 

65  RMPNT=0 

IF  (I.GT.IPPTS)  GOTO  75 
RPROC=PROC (I) 

GOTO  80 
RPROC=0 

RDEV=RDEV+ ABS ( RMPNT-RPROC ) 

IF  (IT.EQ.3)  TYPE"M  PNT :  " , RMPNT , "  P  POINT:  ",RPROC 
IF  (IT.EQ.3)  TYPE "DEV I AT I ON  OF  POINT:  "  ,  I ,  "  IS  " , RDEV 
IF  (IT.EQ.3)  TYPE 
1  =  1  +  1 

GOTO  60 

IF  (IP.NE.l)  GOTO  230 
WRITE (12,200)  ICNT 
WRITE (12,210)  (MTEST (I) ,1=1, ICNT) 

WRITE (12,220)  IPPTS 
WRITE (12 ,210)  (PROC(I) ,1=1, IPPTS) 

FORMAT (IX, "NEW  INTEGRATION  (MODEL):  (PNTS=  ”,I3) 

FORMAT (IX, 8F6. 2) 

FORMAT (IX, "NEW  INTEGRATION  (TEST  CASE):  (PNTS=  ”,I3) 

C  LINEAR  INTERPOLATE  MODEL  DATA  TO  MATCH  TEST  AND  CALC  RDEV 2 

230  RTOTAL=0  ; TOTAL  DEVIATION 

PCURPNT=1  ; POINT  TO  CURRENT  TEST  POINT 
PNXTPNT=2  ; POINT  TO  NEXT  TEST  POINT 
MCURPNT= 1  ; POINT  TO  CURRENT  MODEL  POINT 

MNXTPNT  =  2 

RDISTHI =RPRES  ; TEST  DATA  HIGH  POINT  (INCREMENTED  BY  X  WIDTH) 
RMMF.E  S  =  RMRES  *  INTWD 

IF  (IT.EQ.3)  TYPE"MODEL  RES:  " , RMMRES , "  PROC  RES:  ",RPRES 
240  RTEMP=  MCURPNT*RMMRES  ; CHECK  INTEGRATION  WIN 

IF  ( (RTEMP.GT. RDISTHI) .AND. (MCURPNT.GT. 1) )  GOTO  350 

IF  (IT.EQ.3)  TYPE 


70 


75 

80 


OL 

100 


200 

210 

220 
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IF  (IT.EQ.3)  TYPE"CURRENT  MODEL  POINT:  ", MCURPNT 

IF  (MCURPNT.LE. ICNT)  GOTO  250 
RMODL=0 
GOTO  260 

RMODL=MTEST (MCURPNT) 

IF  (PNXTPNT.LE. IPPTS)  GOTO  270  ; CHECK  IF  END  OF  DATA 

RPROC=0 
GOTO  280 

RHI  =  PROC  (PNXTPNT)  ; CALCULATE  SLOPE 

RLO=PROC (PCURPNT) 

IF  (PCURPNT. LT. 1)  RLO=0 

RM= (RHI-RLO) / RPRES 

IF  (IT.EQ.3)  TYPE"SLOPE :  " , RM 

B=RHI -RM*RPRES*PNXTPNT  ;CALC  INTERCEPT 

IF  (IT.EQ.3)  TYPE" INTERCEPT:  ",B 

PTEST (MCURPNT) =RM*MCURPNT*RMMRES+B 
RPROC=PTEST (MCURPNT) 

RDEV2=ABS (RPROC-RMODL) *ABS (RPROC-RMODL) 

RTOTAL=RTOTAL+RDEV2 

IF  (IT.NE.3)  GOTO  300 

TYPE" INTERP  BETWEEN  PROC  POINTS:  ", PCURPNT , PNXTPNT 
TYPE"PROC  POINT:  ",RPROC 
TYPE "MODEL  POINT:  " , RMODL 

TYPE "DEV  FOR  THIS  POINT:  ",RDEV2,"  TOTAL:  " , RTOTAL 

TYPE"RDISTHI :  ",RDISTHI 

TYPE 

MNXTPNT=MNXTPNT+1 

IF  ( (XTPNT.GT. IPPTS) .AND. (MNXTPNT.GT. ICNT) )  GOTO  400 
MCURPNT =MCURPNT+ 1 

GOTO  240 

INCREMENT  TO  NEXT  WINDOW 

PNXTPNT = PNXTPNT +1 
PCURPNT=PCURPNT+ 1 
RD I STH I = RPRES* PNXTPNT 
GOTO  240 

IF  (IP.NE.l)  GOTO  999 

WRITE (12,500)  ICNT 

WRITE (12,210)  (MTEST(I) ,1=1, ICNT) 

WRITE (12,510)  PCURPNT 

WRITE (12,210)  (PTEST ( I ) , I = 1 , MCURPNT ) 

FORMAT (IX, "LINEAR  INTERPOLAT  MODEL  DATA:  (",I3,"  POINTS)") 
FORMAT (IX, "LINEAR  INTERPOLAT  TEST  DATA:  (",I3,"  POINTS:") 


999 


RETURN 

END 


I  c* 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


VERSION  2  26  JULY  1984 

VERSION  1 


SUBPROGRAM  PIXFRAME  LANGUAGE:  FORTRAN  V 

FUNCTION:  COMPUTES  TARGET  VELOCITY  IN  PIXELS/FRAME 


************ 


********* 


K  X  K  X  W  K 


SUBROUTINE  PIXFRAME (DIR ,MPS , FRATE , LR, PIXFRX , PIXFRY) 

INTEGER  FRATE ,SHX, SHY 
REAL  MPS , LR , VELX , VELY , PI 
REAL  DIR , PIXFRX , PIXFRY 

PI=3. 14159265 

C  GET  COMPONENT  VELOCITIES 

VELX=SIN (DIR) *MPS 
VELY=COS (DIR) *MPS 

C  INSURE  CORRECT  SIGN  ON  VELOCITIES  (FOURTH  QUAD  ORIENTATION) 

IF  ( (DIR. GT. (PI/2) ) .AND. (DIR. LT. (1.5*PI) ) )  GOTO  100 
VELX=ABS ( VELX ) 

GOTO  150 

100  VELX=-ABS (VELX) 

150  IF  ( (DIR.GE.0) .AND. (DIR.LT.  PI))  GOTO  175 

VELY=ABS (VELY) 

GOTO  200 

175  VELY=-ABS (VELY) 

200  CONTINUE 

C  CALCULATE  LOWRES  PIXEL/FRAME  IN  X-Y  DIRECT 

PIXFRY=VELY/ (LR*FRATE ) 

PIXFRX=VELX/ (LR*FRATE) 


RETURN 

END 
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PROGRAM  PIXMOD 


LANGUAGE:  FORTRAN  V 


FUNCTION:  USES  PICBUF  TO  CONVERT  A  STANDARD  256x256x4 

VIDEO  FILE  TO  256x256x16  FOR  USE  BY  PROGRAM  SCENE. FR 

******************************************************* 

INTEGER  IBUF1 ( 300 ) ,  IBUF2(300) 

INTEGER  I ARRAY ( 256 ) 

DIMENSION  INFILE  (13),  OUTFILE(13) 

TYPE 

TYPE" *OCTEK  PICTURE  FILE  CONVERSION*" 

TYPE 

TYPE"ENTER  BACKGROUND  VIDEO  FILENAME:  " 

READ (11,10)  INFILE ( 1 ) 

FORMAT ( S20 ) 

TYPE"ENTER  OUTPUT  FILENAME:  " 

READ (11,10)  OUTFILE ( 1 ) 

TYPE"BU ILDING  BUFFERS..." 

CALL  PICFMT (IBUF1 , INFILE, IHDR) 

CALL  MAKB(IBUFl) 

CALL  PICIN(IBUF1, INFILE, IHDR) 

CALL  PFMT{IBUF2, 256, 256,16,1) 

CALL  MAKB ( IBUF2 ) 


DO  100  1=1,256 

TYPE"LINE :  ",I 

CALL  GROW ( IBUF1 ,1,1,256, I ARRAY) 
CALL  PROW (IBUF2, 1,1, 256, IARRAY) 

CONTINUE 

TYPE"SAVING  PICTURE..." 

IHDR=0 

CALL  PICOUT(IBUF2, OUTFILE, IHDR) 

CALL  RELB ( IBUF1 ) 

CALL  RELB ( I BUF2 ) 


non 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PROGRAM  PROFILE  LANGUAGE:  FORTRAN  V 

FOR  THE  ECLIPSE 


FUNCTION:  EXTRACTS  TARGET  SIGNATURE  FROM  STARING  SENSOR 

DATA. 


MACRO: 


FORTRAN  PROFILE 

RLDR  PROFILE  DECON  PIXFRAME  FRACVPROC  FOUREA 
GCOEFF  GETHDR  GETDA  F5PICBUF.LB  @FLIB@ 


PARAMETER  MAXWIN=10 
PARAMETER  MAXFR=50 
PARAMETER  LETAR=90 
PARAMETER  IFPMAX=20 


C  DECLARE  ARRAYS 


;  Number  of  Sample  Windows 
;  Max  Frames  Program  Can  Process 
;  Largest  Expected  Target 
;  Size  Of  Focal  Plane  (Detectors) 
;  N  x  N 


DIMENSION  INFILE  (13) 

INTEGER  FPLANE ( IFPMAX , IFPMAX) 

INTEGER  IWINPOSX1 (MAXWIN) , IWINPOSX2 (MAXWIN) , IWINPOSY1 (MAXWIN) 
INTEGER  IWINPOSY2 (MAXWIN) ,XP (MAXWIN) ,YP (MAXWIN) 

REAL  RLOW( IFPMAX, IFPMAX) ,RXSEQ (MAXWIN, MAXFR) 

REAL  XM (MAXWIN) ,YM (MAXWIN) , RYSEQ (MAXWIN, MAXFR) 

REAL  XWNOISE (MAXWIN) ,  YWNO I SE (MAXWIN) 

REAL  XAVG (MAXFR) ,YAVG (MAXFR) , RREAL (MAXFR) 


COMPILER  FREE 


CLEAR  MONITER  FOCAL  PLANE  DISPLAY  ARRAY 
DO  5  1=1, IFPMAX 

DO  5  J=l, IFPMAX 

FPLANE  (J, I) ="  " 

5  CONTINUE 

DO  7  1=1, MAXWIN 

DO  7  J=l, MAXFR 

RXSEQ (I , J) =0 
RYSEQ ( I , J) =0 
XAVG ( J) =0 
YAVG(J) =0 
RREAL (J) =0 

7  CONTINUE 

NBND=MAXWIN 

C  GET  INTERACTIVE  DATA 
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IWPRT=0 

ACCEPT"SEND  DATA  TO  PRINTER?  ",IPR 

ACCEPT "CHANGE  MAX  NUMBER  OF  WINDOWS?  :",IANS 

IF  (IANS.NE.l)  GOTO  8 

ACCEPT "ENTER  NUMBER  OF  WINDOWS:  " ,NBND 
ACCEPT "GRAPH  WINDOW  DATA  ?",IANS 
IF  (IANS.NE.l)  GOTO  10 

TYPE "GRAPH  DATA:  1)  BEFORE  PROCESSING" 

TYPE"  2)  AFTER  PROCESSING" 

ACCEPT"  ?" , IWPRT 

GET  SCENE  HEADER  INFORMATION 

ACCEPT "ENTER  SCENE  DATA  FILE  NAME:  " 

READ (11,14)  INFILE ( 1 ) 

FORMAT (SI 3) 

OPEN  2, INFILE 

READ (2)  INUM  ;  NUMBER  OF  FRAMES 

READ (2)  IFPS  ;  FOCAL  PLANE  SIZE 

READ (2)  ISP  ;  SUB  PIXEL  RES 

READ (2)  IFRATE  ;  IFRATE 

READ (2)  BLUR  ;  OTF  STD  DEV 

READ  (2)  DR  ;  DRIFT  RATE 

READ (2)  FPRINT  ;  FOOT  PRINT 

READ (2)  ICTX  ;  INITIAL  TARG  POSITION 

READ (2)  ICTY  ; 

READ (2)  RMPS 
READ (2)  RDIR 

READ (2)  TRES  ;  PLANE  RES 

READ (2)  ITH  ;  PLANE  THRESHOLD 

READ (2)  PXDST  ;  VEL  IN  PIX/ FRAME 

IFPX=IFPS 

IFPY=IFPS 

PRELIMINARY  VELOCITY  VECTOR  ACQUISITION 


DETERMINE  PHASE  RELATIONSHIP 
-  SETUP  WINDOW  CONSTRAINTS 


GET  TARGET  SHIFT 

CALL  PIXFRAME (RDIR , RMPS , IFRATE , FPRINT , SHX , SHY ) 
TYPE"P/F  X  :  "  ,  SHX  ,  "  Y:",SHY 
RESX=ABS (SHX*FPRINT) 

RESY=ABS (SHY*FPRINT) 

PIXFRX=ABS (1/SHX) 

PIXFRY=ABS (1/SHY) 


o  n  o 


TYPE"F/P  X : " , PIXFRX , "  Y:  n,PIXFRY 

IPIXFRX=IFIX (PIXFRX) 

IPIXFRY=IFIX (PIXFRY) 

FRACX=ABS (PIXFRX-IPIXFRX) 

FRACY  =  ABS (PIXFRY-IPIXFRY) 

TYPE"FRACX  :  ",FRACX,"  FRACY FRACY 

IF  ( (FRACX.GE.0.2) .AND. (FRACX.LE.0.8) )  GOTO  40 
IWPSX=1 
IWINTX=1 
GOTO  42 

40  CALL  FRAC (FRACX , IWPSX , IWINTX ,ERRORX ) 

42  IF  (  (FRACY. GF,. 0.2)  .AND.  (FRACY. LE. 0.8)  )  GOTO  44 

IWPSY=1 
IWINTY=1 
GOTO  46 

44  CALL  FRAC (FRACY, IWPSY, IWINTY,ERRORY) 

46  TYPE "FOR  X  DIR:  WIN/ SET= " , IWPSX , "  NUMERAT IWINTX 

TYPE"FOR  Y  DIR:  WIN/SET=" , IWPSY , "  NUMERAT IWINTY 

IWPSX= IFIX (IWPSX+.5) 

IWPSY=IFIX ( IWPSY+ . 5) 

NSETX=NBND/ IWPSX 
NSETY=NBND/ IWPSY 

TYPE"SETS  OF  WINDOWS:  X:  ",NSETX,"  Y:  " ,NSETY 


C  CALCULATE  WINDOW  FOR  DECONVOLUTION 

IWINCX=IFIX (ABS (1/SHX) *IWPSX+.5) 

IWINCY=IFIX (ABS (1/SHY) *IWPSY+. 5) 

TYPE "DECONVOLUTION  WINDOW:  X:",IWINCX,"  Y:",IWINCY 
RESX= (RESX/ IWPSX) 

RESY= (RESY/ IWPSY) 

TYPE "NEW  RES  X:  ",RESX,"  Y:  " , RESY 
TYPE 

DETERMINE  WINDOW  SIZES 

IF  (RDIR.EQ. 0 )  RDIR=0 . 01 
IF  (RDIR.GT. 1.5707)  RDIR=1.5707 

C  SET  INITIAL  WINDOW  SIZE  BY  DETERMINING  FLITE  INTERSECTION 

DX=FLOAT (LETAR) / (SIN (RDIR) *2*FPRINT) 

DY=FLOAT (LETAR) / (COS (RDIR) *2*FPRINT) 

IWINX  =  INT ( 2  *DX+ . 5 ) +2 
IWINY=INT(2*DY+.5) +2 


INSURE  WINDOW  LENGTH  ON  FOCAL  PLANE 

IF  (IWINX.LT. 1)  IWINX=1 
IF  (IWINY.LT. 1)  IWINY=1 
IF  (IWINX.GT. IFPX)  IWINX= IFPX 
IF  (IWINY.GT. IFPY)  IWINY=IFPY 


TYPE"WINDOW  SIZE 


X : " , IWINX , "  Y : " , IWINY 


DETERMINE  X-WINDOW  POSITIONS 


C  GET  FIRST  CENTER 

XDIR=TAN ( (3.14159265/2) -RDIR) *LETAR/ (2*FPRINT) 

IXCENT “ ICTX+ IFIX (FLOAT (LETAR) /  (2*FPRINT) +.5) 

IYCENT=ICTY-IFIX (XDIR) +1 

C  INSURE  CENTER  ON  FOCAL  PLANE 

IF  (IXCENT. LT.l)  IXCENT=1 
IF  (IYCENT.LT. 1)  1 YCENT= 1 
IF  (IXCENT. GT. IFPX)  IXCENT=IFPX 
IF  (IYCENT.GT. IFPY)  IYCENT=IFPY 

C  STORE  FIRST  CENTER 

IWINPOSX1 ( 1) =INT ( IXCENT+ .5) 

IWINPOSX2 ( 1 ) =  INT ( I YCENT+ .5) 

C  GET  REMAINING  CENTERS  RELATIVE  TO  FIRST  CENTER 

RXCENT= FLOAT ( IWINPOSX1  ( 1 ) ) 

RYCENT= FLOAT ( IWINPOSX2 ( 1 ) ) 

DO  100  1=2 ,NBND 

75  RXCENT=RXCENT+SIN (RDIR) 

RYCENT=RYCENT-COS (RDIR) 

IX= INT ( RXCENT+ . 5 ) 

IY=INT (RYCENT+ . 5 ) 

IF  (IX. GT. IFPX)  IX=IFPX 
IF  (IY.GT.IFPY)  I Y= IFPY 
IF  (IX. LT.l)  IX=1 
IF  (IY.LT.l)  IY=1 

IF  ( (IX.EQ. IWINPOSX1 (1-1) ) .AND. (IX.LT. IFPX) )  GOTO  75 
I WINPOSX1 (I ) =IX 
IWINPOSX2 (I) =IY 

100  CONTINUE 

DO  110  1=1 ,  NBND 

TYPE"X  WIN  COORD:  X : " , IWINPOSX1 ( I )  , "  Y : " , IWINPOSX2  ( I ) 
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110  CONTINUE 

£******************************************************************** 

C 

C  DETERMINE  Y-WINDOW  POSITIONS 

C 

£******************************************************************** 

C  GET  FIRST  Y-CENTER 

YDIR=LETAR/ ( 2 *FPRINT*TAN ( ( 3 . 14 15 9 2 6 5/ 2 ) -RDIR) ) 

IXCENT= ICTX+ IFIX (YDIR) -1 

IYCENT=ICTY-IFIX (FLOAT (LETAR) / (2*FPRINT) +.5) 

C  INSURE  CENTER  ON  FOCAL  PLANE 

IF  (IXCENT.LT. 1)  IXCENT=1 
IF  (IYCENT.LT. 1)  IYCENT=1 
IF  (IXCENT.GT. IFPX)  IXCENT=IFPX 
IF  (IYCENT.GT. IFPY)  I YCENT= IFPY 

C  STORE  Y-CENTER 

IWINPOSY1 (1) =IXCENT 
IWINPOSY2 ( 1 ) = INT ( I YCENT+ . 5 ) 

C  GET  REMAINING  Y-WINDOW  CENTERS 

RXCENT  =  FLOAT (IWINPOSY1  (1) ) 

RYCENT= FLOAT ( IWINPOSY2  ( 1 ) ) 

DO  200  1=2 ,  NBND 

175  RXCENT=RXCENT+SIN (RDIR) 

RYCENT=RYCENT-COS (RDIR) 

IX=INT (RXCENT+ . 5) 

IY=INT (RYCENT+ . 5) 

IF  (IY.EQ. IWINPOSY2 (1-1) )  GOTO  175 
IWINPOSY1  (I) =IX 
IWINPOSY2  (I) =IY 

200  CONTINUE 

DO  210  1=1, NBND 

TYPE"Y  WIN  COORD  X : " , IWINPOSY1 ( I ) , "  Y: " , IWINPOSY2 ( I ) 

210  CONTINUE 

Q  ******************************************************************* 

C 

C  LOOP  HERE  TO  ACQUIRE  AND  PROCESS  DATA 

C 

Q********* *********************************************************** 

DO  1000  IFR=1 , INUM 

C  GET  FRAME  OF  DATA  AND  DISPLAY 


o  o  o  o  o 


TYPE 'READING  DATA - ' 

DO  250  1=1,  IFPX 

DO  250  J=l,  I FPY 

READ  (2)  PIX 
RLOW ( J , I ) =PIX 

250  CONTINUE 

C 

C  SUM  X-WINDOW  DATA 

C 


TYPE"SUMMING  X-WINDOWS _ " 

DO  300  1=1 ,NBND 

IXEND= IWINPOSX2 (I ) +INT (LETAR/ (SIN (RDIR) *2*FPRINT) ) 

IF  ( IXEND.GT. I FPY)  IXEND=IFPY 

I X ST ART = I XEND- IWINX 

IF  ( IXSTART . GT . I FPY)  IXSTART=IFPY 

IF  (IXSTART. LT.l)  IXSTART=1 

DO  300  K=IXSTART, IXEND 
IF  (K.LT.l)  GOTO  300 
FPLANE ( IWINPOSX1 ( I ) ,K) ="XX" 

RXSEQ ( I , IFR) =RXSEQ ( I , IFR) +  RLOW ( IWINPOSX1 (I) ,K) 
IF  (IFR.NE.l)  GOTO  300 

XWNOISE ( I ) = RXSEQ (1,1) 


300  CONTINUE 

SUM  Y-WINDOW  DATA 


TYPE"SUMMING  Y-WINDOWS - " 

DO  400  1=1 ,NBND 

I YSTART= IWINPOSY 1 (I ) -INT (LETAR/ (COS (RDIR) *2*FPRINT) ) 

IF  (IYSTART.LT. 1)  I YSTART=1 

IF  (IYSTART.GT. IFPX)  IYSTART=IFPX 

I YEND= I YSTART+ I WINY 

IF  (IYEND.GT. IFPX)  IYEND=IFPX 

IF  (IYEND.LT. 1)  I YEND= 1 

DO  400  K=IYSTART, IYEND 

FPLANE (K, IWINPOSY2 (I) ) ="YY" 

RYSEQ ( I , I FR) =RYSEQ ( I , IFR) +RLOW ( K , IWINPOSY2 ( I ) ) 

IF  (IFR.NE.l)  GOTO  400 

YWNOISE (I) = RYSEQ (1,1) 


400 


CONTINUE 
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RTX= FLOAT ( ICTX) + ( IFR-1 ) *SHX 
RTY= FLOAT ( ICTY) + ( IFR-1 ) *SHY 
IX= INT (RTX  + .5) 

IY= INT (RTY+ .5) 

IF  (  (IX. EQ. ICTX)  .AND.  (IY.EQ. ICTY)  )  GOTO  500 
IF  ( ( IX.GT. IFPX) .OR. ( IY.GT. IFPY) )  GOTO  500 
IF  ( (IX.LT.l) .OR. (IY.LT.l) )  GOTO  500 
FPLANE ( IX , I Y ) =  "TT" 

500  CONTINUE 

TYPE "FOCAL  PLANE  BOUNDRIES  .." 

FPLANE (ICTX, ICTY) =  "**" 

WRITE (10,505)  (I, 1=1, IFPX) 

505  FORMAT (4X, 201 3) 

DO  600  1=1, IFPY 

WRITE (10,510)  I,  (FPLANE  ( J  ,  I )  ,J=1,IFPY) 

510  FORMAT (IX, 13, 2X, 20  (A1,2X) ) 

600  CONTINUE 

TYPE"FRAME :  ",IFR 

1000  CONTINUE 

IF  (IWPRT.NE.l)  GOTO  800 

DO  750  1=1,5 

flt  DO  725  J=1 , INUM 

XAVG ( J / =RXSEQ  ( I  ,  J) 

725  CONTINUE 

CALL  PLOT 10 (XAVG, INUM,I,X0,Y0 ,RESX) 

750  CONTINUE 

DO  780  1=1,5 

DO  760  J= 1 , INUM 

YAVG (J) =RYSEQ ( I , J ) 

760  CONTINUE 

ITEMP= 1+5 

CALL  PLOT 10 (YAVG, INUM, ITEMP , X0 , Y0 , RESY) 

780  CONTINUE 

c* ***********  ******************************************************** 

C 

C  NORMALIZE  DATA 

C 

c******************************************************************** 

800  DO  1020  1=1, NBND 

RXMAX=0 

RYMAX=0 

DO  1020  J= 1 , INUM 

IF  (RXSEQ ( I , J) . LE . RXMAX )  GOTO  1005 
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1005 

1015 

1020  CONTINUE 


RXMAX=RXSEQ ( I  ,  J ) 

IF  (RYSEQ (I , J) .LE.RYMAX)  GOTO  1015 
RYMAX=RYSEQ (I , J) 

XM (I) =RXMAX 
YM ( I ) =RYMAX 


C 

C 

C 

C 

C1100 

c 

c 

c 

c 

c 

C1200 


WRITE  (10,9000) 

WRITE (10,9005) 

DO  1100  1=1 ,NBND 

WRITE (10 , 9010)  ( RXSEQ ( I , J) ,J=1,INUM) 

CONTINUE 
WRITE (10,9005) 

WRITE (10,9015) 

WRITE (10,9005) 

DO  1200  1=1 ,NBND 

WRITE (10 , 9010)  (RYSEQ ( I , J) ,J=1,INUM) 

CONTINUE 


9000 

9005 

C9010 

9015 

9018 

9020 

9025 

9027 

9030 

9031 

9032 

9034 

9036 

9037 

9038 

9035 
9040 
9042 
9045 
9050 
9052 
9C55 


FORMAT ( 
FORMAT ( 
FORMAT 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 


IX," 


X  WINDOW  PROFILE:  ") 


") 


( IX , 2  0F4 . 0 ) 
IX,"  Y 

IX, "WINDOW: 
IX , 5F 1 2 • 6 ) 
IX , "PROFILE 


WINDOW 

",  I  3 ) 


PROFILE:  ") 


TARGET  VECTORS  FOR  FILE:  ",5A2) 


IX, "FILE  HEADER  INFORMATION:  ") 


IX, "FILE  CONTAINS: 

IX, "FOCAL  PLANE  SIZE: 
IX, "SUBPIXEL  SIZE: 

IX, "BLUR  DEVIATION: 

IX, "DRIFT  RATE: 

IX, "PLANE  RESOLUTION: 
IX, "TARGET  THRESHOLD: 
IX, "TARGET  DIRECTION: 
IX, "TARGET  VELOCITY: 
IX, "PIXEL  VELOCITY: 

IX, "FRAME  RATE: 

IX, "FOOTPRINT: 

IX , "CONV  WINDOW:  X: 

IX, "INITIAL  POSITION: 


" , 1 1 2  ,  "  FRAMES") 

"  ,  1 1 2  ,  "  DETECTORS") 

" , I 1 2 , "  SUB  PIXELS") 

"  ,  F 1 2 . 6  ) 

" , F 1 2 . 6 ,  "  METERS/S") 

" , F 1 2 . 6  ,  "  METERS") 

" ,112,"  GRAYLEVEL" ) 
",F12.6,"  RADIANS") 

" , F 1 2 . 6  ,  "  METERS/S") 

" , F 1 2 . 6  ,  "  PIX/FRAME") 

" , 1 1 2  ,  "  FRAMES/S") 

" , F 1 2 . 6  ,  "  METERS") 

"16,"  Y :  " , 1 6 ) 

” , 1 6 , 1  6  ,  "  DETECTOR (X,Y) ") 


IF  (IPR.NE.l)  GOTO  9250 


WRITE (12,9025) 
WRITE (12,9005) 
WRITE (12,9005) 
WRITE (12,9027) 
WRITE ( 12 , 9005) 
WRITE (12,9030) 
WRITE (12,9031) 


( INFILE (L) , L= 1 , 5 ) 


INUM 
I FPS 
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o  o  o  o  o 


WRITE (12,9032)  ISP 
WRITE (12,9045)  IFRATE 
WRITE (12,9050)  FPRINT 
WRITE (12,9034)  BLUR 
WRITE (12,9036)  DR 
WRITE (12 ,9005) 

WRITE (12,9055)  ICTX,ICTY 
WRITE (12,9040)  RMPS 
WRITE (12,9035)  RDIR 
WRITE (12,9037)  TRES 
WRITE (12,9038)  ITH 
WRITE (12,9042)  PXDST 
WRITE (12,9005) 

WRITE (12,505)  (I,I=1,IFPX) 

DO  620  1=1 , IFPY 

WRITE (12, 5 10)  I,  (FPLANE  ( J  ,  I )  ,J=1,IFPX) 

620  CONTINUE 

WRITE (12,9052)  IWINCX , IWINCY 
WRITE (12,9005) 

WRITE (12,9005) 

WRITE (12,9000) 

WRITE (1 2 ,9005) 

DO  9100  1=1, NBND 

WRITE (12 , 9018)  I 

WRITE  (12,9020)  (RXSEQ ( I , J ) , J= 1 , INUM ) 

9100  CONTINUE 

WRITE  (12,9005) 

WRITE (12,9015) 

WRITE (12,9005) 

DO  9200  1=1, NBND 

WRITE (12,9018)  I 

WRITE (12 ,9020)  (RYSEQ (I , J) , J=1 , INUM ) 

9200  CONTINUE 

***************>***★*****************★*******************★********* 

PROCESS  VECTOR  PROFILES  TO  EXTRACT  TARGET  INFORMATION 

******************************************************************* 


C  STRIP  OFF  DC  NOISE 

9250  DO  9300  1=1, NBND 

IXP=1000 
I YP= 1000 
XMAX=0 
YMAX=0 

DO  9300  K=1 , INUM 

RXSEQ ( I , K ) =RXSEQ ( I , K) -XWNOISE ( I ) 
IF  ( RXSEQ  ( I , K )  . LE . XMAX )  GOTO  9  2  30 
IF  (K.LT.IXP)  IXP=K 
XM (I) =RXSEQ(I ,K) 


151 


9230 


9260 


XMAX=XM (1 ) 

RYSEQ ( I ,K) =RYSEQ ( I , K ) -YWNOI SE ( I ) 
IF  (RYSEQ(I,K) .LE.YMAX)  GOTO  9260 
IF  (K.LT.IYP)  I YP=K 
YM ( 1 ) = RYSEQ  ( I ,  K ) 

YMAX  =  YM  (I) 

XP(I) =IXP 
YP  ( I ) = IYP 


9300  CONTINUE 


9310 


9400 


9450 


IF  (IPR.NE.l)  GOTO  9450 
WRITE (12,9005) 

WRITE (12,9310) 

FORMAT (IX,"  PROFILE  WITHOUT  DC") 

WRITE (12,9000) 

DO  9400  1=1 ,NBND 

WRITE (12,9018)  I 

WRITE (12,9020)  (RXSEQ(I,K) ,K=1,INUM) 


CONTINUE 
WF-ITE  (12,9015) 

DO  9450  1=1, NBND 

WRITE (12,9018)  I 

WRITE  (12,9020)  ( RYSEQ ( I , K )  ,K  =  1,INUM) 

CONTINUE 


C  NORMALIZE  CONVOLUTED  VECTORS 


C  DO  9500  1=1, NBND 

C  DO  9500  K  =  1 , I NUM 

c  RXSEQ ( I , K ) =RXSEQ ( I , K ) /XM ( I ) 

c  RYSEQ ( I , K ) =  RYSEQ (I , K ) / YM ( I ) 

C9500  CONTINUE 
C 

C  IF  (IPR.NE.l)  GOTO  9550 

C  WRITE (12,9005) 

C  WRITE (12,9510) 

C9510  FORMAT (IX,"  NORMALIZED  WITHOUT  DC") 

C  WRITE (12,9000) 

C  DO  9530  1=1, NBND 

C  WRITE (12,9018)  I 

C  WRITE (12,9020)  ( RXSEQ ( I , K ) , K= 1 , INUM ) 

C9530  CONTINUE 
C  WRITE (12 ,9015) 

C  DO  9550  1=1, NBND 

C  WRITE (12 ,9018)  I 

C  WRITE (12,9020)  (RYSEQ ( I , K)  ,  K  =  1 , INUM) 

C9550  CONTINUE 


C 

C  PROCESS  LOW  SNR  VECTORS 

C 
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IF  (ICON.NE.l)  GOTO  594 
DO  593  1=1,1 XLEN 

IF  (CONVX(I) .GT.X1MAX)  X1MAX=C0NVX ( I ) 

593  CONTINUE 

DO  594  1=1 , IYLEN 

IF  (CONVY (I) .GT.Y1MAX)  Y1MAX=C0NVY ( 1 ) 

594  CONTINUE 

DO  595  1=1, IXL 

RXVECT ( I ) =  FLOAT ( XVECT ( I )  ) / XMAX 

595  CONTINUE 

DO  596  1=1, IYL 

RYVECT ( I ) =FLOAT (YVECT ( I ) ) / YMAX 

596  CONTINUE 

IF  (ICON.NE.l)  GOTO  598 
DO  597  1=1 , I XLEN 

RCONVX ( I ) =FLOAT (CON’  .  ' ) /X1MAX 

597  CONTINUE 

DO  598  1=1, IYLEN 

RCONVY ( I ) =FLOAT (CONVY (I ) ) / Y1MAX 

598  CONTINUE 

WRITE  (12,310) 

WRITE (12,350)  IXNUM 

WRITE (12 ,599)  (RXVECT(I) ,1=1, IXL) 

WRITE (12,370)  IYNUM 

WRITE (12,599)  (RYVECT(l) ,1=1, IYL) 

IF  (ICON.NE.l)  GOTO  600 
WRITE (12,310) 

WRITE (12,599)  (RCONVX (I) ,I=1,IXLFN) 

WRITE (12,511) 

WRITE (12,599)  (RCONVY(I) ,1=1, IYLEN) 

599  FORMAT (1X,10F8.5) 

C 

C  GRAPH  OUTPUT 

C 

600  ACCEPT"  DO  YOU  WANT  TO  GENERATE  PLOTS  OF  OUTPUT?  " , IGRPH 
IF  (IGRPH.NE.l)  GOTO  750 

C  ACCEPT " PLOT  (1)  CONVOLVED  OR  (2)  NORMAL  DATA?  ",ICN 

C  ACCEPT "PLOT  (1)  X  AXIS  OR  (2)  Y  AXIS  ?  ",IAXIS 

C 

C  TYPE"ENTER  TITLE  OF  GRAPH:  " 

C  READ (11,3)  INFILE(l) 

C 

C  TYPE"FLIP  TERMINAL  SWITCH... AND  HIT  ANY  KEY..." 

C  PAUSE 

C  NG  =  1 

C  MD  =0 

C  YN  =  0 

C  YX  =  0 

C  I FS=  0 

C  IF  ( (ICN.NE.l) .AND. (IAXIS.NE.l) )  GOTO  760 

C  CALL  GR  PH2 (INFILE, NG , YVECT, XVECT, IYNUM, MD , YN , YX , IFS ) 

C60  IF  ( (ICN.NE.l) .AND. (IAXIS.EQ.l) )  GOTO  770 
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330 

FORMAT ( IX, "TARGET  WINDOW: 

"  ,  1 4 , "  BY  " ,14) 

340 

FORMAT (1X,"TARGET  RESOLUTION: 

" , F 1 0 . 6  ,  "  METERS" 

3  50 

FORMAT (IX, "X-VECTOR:  ",I4," 

ELEMENTS") 

360 

FORMAT (2014) 

370 

FORMAT (IX, "Y-VECTOR:  ",I4," 

ELEMENTS" ) 

380 

FORMAT (2014) 

c 

C  CONVOLVE  DATA 

C 

400  ACCEPT"DO  YOU  WANT  TO  CONVOLVE  DATA?  " , ICON 

IF  (ICON . NE . 1 )  GOTO  590 

ACCEPT " ENTER  WINDOW  SIZE->" , IWIN 
IXLEN= IXL+IWIN-1 
I YLEN  =  I YL+ 1 W I N - 1 
DO  500  1=1, IXLEN 

DO  500  J= 1 , IWIN 

IF  ( (I-J+l) .LT.l)  GOTO  500 
CONVX (I ) =CONVX (I) +XVECT (I-J+l) 

500  CONTINUE 

DO  505  1=1,1 YLEN 

DO  505  J=1 , IWIN 

IF( (I-J+l) .LT.l)  GOTO  505 
CONVY ( I ) =CONVY ( I ) + YVECT (I-J+l ) 

505  CONTINUE 

WRITE (10,510) 

WRITE (10,380)  (CONVX (I) ,1=1, IXLEN) 

510  FORMAT ( IX, "X  VECTOR  CONVOLUTION") 

WRITE (10,511) 

511  FORMAT (1X,"Y  VECTOR  CONVOLUTION") 

WRITE (10,380)  (CONVY  (I)  ,I  =  1,IYLEN) 

IF  (IPF.NE.l)  GOTO  575 
WRITE (12,510) 

WRITE (12,380)  (CONVX (I) ,1=1, IXLEN) 

WRITE (12,511) 

WRITE  (12,3  80)  (CONVY (I)  ,1  =  1, 1 YLEN) 

590  ACCEPT " PRINT  NORMALIZED  VECTORS?  " , IANS 

IF  (IANS.NE.l)  GOTO  600 

XMAX=0 
YM AX  =  0 
X1MAX=0 
Y1MAX=0 

DO  591  1=1, IXL 

IF  (XVECT (I ) .GT.XMAX)  XMAX  =  XVECT  ( I ) 
551  CONTINUE 

DO  592  1=1, IYL 

IF  (YVECT (I) .GT.YMAX)  YMAY= YVECT ( I ) 
592  CONTINUE 
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DO  67  1  =  1, IXL 

WRITE  BINARY ( 3 )  HIRES(I,J) 

67  CONTINUE 


55 


50 

60 


70 

100 


200 


300 

310 

320 


TYPE "CALCULATING  X-VECTOR  ..." 

DO  60  J=1 , IXL 

DO  50  1=1 ,NUMPOINTS 

IF  (TARGETX (I) .EQ. J)  XVECT ( J ) =XVECT ( J )  + 1 5 -TARGET I ( I ) 
IF  (TARGETX  (I)  .GT.l)  GOTO  50 
TYPE "POSSIBLE  NOISE  DATA  POINT:" 

TYPE "TARGET  POINT:  ", TARGETX ( I ), TARGETY ( I ),"=", TARGETI ( I ) 
CONTINUE 

IF  (XVECT (J) .GT.0)  IXNUM=IXNUM+1 

CONTINUE 

TYPE "NUMBER  X-ELEMENTS=  ",IXNUM 
TYPE "CALCULATING  Y-VECTOR:  " 

DO  100  J=1 , IYL 

DO  70  1=1 ,NUMPOINTS 

IF  (TARGETY (I) .EQ.J)  YVECT ( J ) =YVECT ( J ) + 1 5 -TARGET I ( I ) 
IF  (TARGETY  (I)  .GT.l)  GOTO  70 
IF  (TARGETY (I) .LT. (IYL-1) )  GOTO  70 
TYPE "POSSIBLE  NOISE  DATA  POINT:" 

TYPE "TARGET  POINT:  ", TARGETX ( I ), TARGETY ( I ),"=", TARGETI ( I ) 
CONTINUE 

IF  (YVECT (J) .GT.0)  I YNUM= I YNUM+ 1 

CONTINUE 

TYPE "NUMBER  OF  \ -ELEMENTS3  ",IYNUM 
IDIR  =  IFIX (180/3. 1 41 592  65  *FDIR+ .01) 

TYPE 

TYPE"DATA  FOR  TARGET  FILE:  " 

TYPE 

TYPE" FACING  DIRECTION:  " , IDIR 
TYPE"WINDOW :  ",IXL,"  BY  ",IYL 

TYPE "RESOLUTION:  ",RESOLU 

TYPE 

TYPE "X-VECTOR:  " 

WRITE (10,200)  (XVECT  (I)  ,1  =  1, IXL) 

TYPE "Y-VECTOR:  " 

WRITE (10,200)  ( YVECT  ( I )  ,1  =  1,1 YL ) 

FORMAT (2014) 

IF  (IPR.NE.l)  GOTO  400 


WRITE (12,310) 
WRITE (12,350) 
WRITE (12,360) 
WRITE (12 ,370) 
WRITE (12,380) 


IXNUM 

(XVECT (I) ,1=1, IXL) 
IYNUM 

(YVECT ( I ) ,1=1, IYL) 


FORMAT ( IX , "VECTOR  DATA  FOR  TARGET  MODEL  DATA  FILE:  ",S20) 
FORMAT '80X) 

FORMAT ( IX, "TARGET  DIRECTION:  ",I6,"  DEGREES") 


c 


OUTPUT  HEADER 


IF  (IPR.NE.l)  GOTO  42 

WRITE (12,300)  INFILE  (1) 

WRITE (12,310) 

WRITE  (12,320)  IDIR 
WRITE (12,330)  IXL,IYL 
WRITE (12,340)  RESOLU 
WRITE (12,310) 

42  TYPE "CLEARING  VECTOR  MATRICES..." 

DO  45  1=1,70 

XVECT (I ) =0 
CONVY(I) =0 
CONVX  (  I ) =0 
YVECT (I) =0 

45  CONTINUE 

ACCEPT "CHANGE  GREYSCALE?  ",IANS 
IF  (IANS.NE.l)  GOTO  575 

ACCEPT "ENTER  GREYSCALE  THRESHOLD->" , ITHRES 
DO  47  1=1 ,NUMPOINTS 

IF  (TARGETI (I) .GT. ITHRES)  TARGETI ( I )= ITHRES 
47  CONTINUE 

C  SET-UP  2-D  ARRAY 

575  DO  578  1=1,70 

DO  578  J=1 , 70 

HIRES  (I , J) =0 

578  CONTINUE 

DO  580  1=1 ,NUMPOINTS 

HIRES (TARGETX (I )  , TARGETY  ( I )  ) = 15-TARGET I  ( I ) 

580  CONTINUE 

ACCEPT "PRINT  2-  PLOT  OF  TARGET  FILE?  ",IANS 
IF  (IANS.NE.l)  GOTO  65 
WRITE (12,310) 

WRITE (12,584)  (INFILE(I)  ,1  =  1,5)  , ITHRES 
DO  582  J= 1 , I YL 

WRITE  (12 , 585)  (HI  RES ( I, J)  ,I  =  1,IXL) 

582  CONTINUE 

584  FORMAT (IX, "HIRES  PLOT  OF  FILE:  ",5A2,"  THRESHOLD:  ”,I3) 

585  FORMAT (IX, 3413) 

65  ACCEPT "GENERATE  BINARY  FILE?  ",IANS 

IF  (IANS.NE.l)  GOTO  55 
ACCEPT "ENTER  FILENAME  ->" 

READ (11,3)  INFILE ( 1 ) 

CALL  OPEN (3, INFILE, 3, IER) 

DO  67  J= 1 , I YL 
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TYPE 

ACCEPT "DO  YOU  WANT  ANOTHER  RUN?  " ,IANS 
IF  (IANS.EQ.l)  GOTO  15 

IF  (IPR.EQ.l)  CLOSE  3 
CALL  RELB(IBUF) 

END 


c******************************************************************** 

C 

C 

C  TARGET  EVALUATOR  : TARGEVAL  VERSION  2.0 

C 

C  FUNCTION:  THIS  PROGRAM  GENERATES  X  AND  Y  VECTOR  PROFILES 

C  FOR  TARGET  DATA  FILES,  AND  EITHER  SAVES  DATA: 

C  1)  LINE  PRINTER 

C  2)  DATABASE 

C 

DIMENSION  INFILE (20) 

INTEGER  NUMPOINTS,  TARGETX  (  9 00 )  , TARGETY ( 900 )  , TARGETI ( 9 00 ) 
INTEGER  XVECT ( 70 ) ,  YVECT(70),  CONVX ( 70 ) , HIRES ( 70 , 70 ) 

INTEGER  CONVY ( 70 ) 

REAL  RXVECT (70) ,RYVECT (70) ,RCONVX (70) ,RCONVY (70) 


1 


m 


3 


IXNUM=0 
I YNUM  =  0 
TYPE 

TYPE"ENTER  TARGET  FILENAME:  ’’ 
READ (11,3)  INFILE  ( 1 ) 

FORMAT (SI  3) 


ACCEPT" SEND  DATA  TO  PRINTER?  ",IPR 
TYPE "OPENING  FILE..." 

CALL  OPEN(4, INFILE, 1,IER) 

IF  (IER.NE.l)  TYPE  "CANT  OPEN  TARGET  FILE:  ",IER 
IF  (IER.NE.l)  STOP 


READ (4)  NUMPOINTS 
READ (4)  FDIR 

TYPE"#  OF  POINTS:  ", NUMPOINTS  ,  "  DIRECTION:  " , F  D I R 
READ (4)  RESOLU 
READ (4)  IXL 
READ (4)  IYL 

TYPE"WINDOW :  ", IXL, IYL 
DO  40  1=1, NUMPOINTS 

READ (4)  TARGETX  (I) 

READ (4)  TARGETY  (I) 

READ (4)  TARGETI  (I) 

40  CONTINUE 

CALL  CLOSE ( 4 , IER) 

TYPE "FINISHED  READING  DATA..." 
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I YEND= IBKOFFY+WINSI ZE-1 
ISCALE=DETS1ZE 

DO  200  J= I YSTART , I YEND , I  SCALE 

C  FILL  BUFFER  WITH  ISCALE  ROWS  OF  DETECTORS 

DO  175  L=l, ISCALE 
IROW= J+L-l 
ICOL= IXSTART 

CALL  GROW { I BUF , I ROW , I COL , WINS I ZE , I DETBUF ( 1 ,  L) ) 

175  CONTINUE 


INTEGRATE  ROWS 


IBND=WINSI ZE-1 
DO  190  1=1 , IBND, ISCALE 
IPIX=0 

DO  180  K=l, ISCALE 

DO  180  L=1 , ISCALE 
IROW=I+K-l 
ICOL=L 

IPIX=IPIX+1DETBUF (I ROW, ICOL) 

CONTINUE 

RPIX=FLOAT ( IPIX) /FLOAT (DETSIZE*DETSIZE) 

IF  (IPR.EQ.l)  WRITE ( 3 )  RPIX 
ICTR=I/ ISCALE+1 
RDETROW ( I CTR ) =RPIX 

CONTINUE 

WRITE (10,189)  (RDETROW (M) ,M=1 ,NDET) 

FORMAT (10X,20F3.0) 

CONTINUE 

IF  ( I  PAUSE . EQ . 1 )  PAUSE 

CALCULATE  NEW  TARGET  POSITION 

IF  (ITAR.NE.l)  GOTO  500 
DELTAX=SIN (RDIR) *PIXDIST 
DELTAY= (-1) *COS (RDIR) *PIXDIST 

CXPOS=CXPOS+DELTAX 

CYPOS=CYPOS+DELTAY 

DP  I  FT  SCENE 

INSX=IDSX+INSX 

CALL  PICIN (I BUF, BKGFILE , IHDR) 


CONTINUE 
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WRITE (3) 

RMPS 

WRITE  (3) 

RDIR 

WRITE (3) 

TRES 

WRITE  (3) 

ITGRAY 

;  PLANE  THRESHOLD 

WRITE (3) 

PIXDIST 

CONTINUE 

GENERATE  LOW  RESOLUTION  SCENE  SEQUENCES 


DO  500  ICNT= 1 , INUM 

TYPE "DISPLAYING  FRAME:  ",ICNT 

MASK  BUFFERED  BACKGROUND  WITH  TARGET 

IF  ( ITAR.NE . 1 )  GOTO  45 

IXPOS=IFIX (CXPOS) 

IYPOS=IFIX (CYPOS) 

TYPE "CURRENT  TARGET  POSITION  X:",IXPOS,"  Y:",IYPOS 
DO  45  1=1 ,NPTS 

IX=TARGETX (I) +IXPOS 
IY=TARGETY (I) +IYPOS 
IV= 15 -TARGET I  (1) 

IF  (IV.LT.l)  GOTO  45 
IF  ( IV. GT. 15)  IV=15 
IF  (ITGRAY.LT.O)  GOTO  40 
IF  (IV.LT. ITGRAY)  IV=ITGRAY 
CALL  PPNT ( IBUF , IY , IX , I V) 

CONTINUE 


CONVOLVE  WITH  POINT-SPREAD  FUNCTION 


IF  (NPSF.NE.l)  GOTO  192 
TYPE "CONVOLVING  PSF  ..." 

TYPE "CALCULATING  COEFF ,  M  =  "  , RMEAN , "  SIGMA=",SIG 
CALL  GCOEFF ( S IG , RMEAN , DETS IZE , COEFF) 

TYPE"CALCULATINE  OTF  ..." 

CALL  PSPREAD (DETSIZE, IBKOFFX, IBKOFFY, INSX , WINS IZE , IBUF, COEFF) 


INTEGRATE  HIRES  WINDOW  AND  DISPLAY  AS  FOCAL  PLANE  PIXELS 


I  XSTART= IBKOFFX  +  INSX 
IXEND  =  I BKOFFX  +WINSI ZE+ INSX 
I YST ART  =  IBKOFFY 


UUUfN  CJ  uuu 


RDIR=RDIR*0 .01745329 


ITGRAY=-1 

ACCEPT" DO  YOU  WANT  TO  CHANGE  GRAYSCALE  OF  TARGET?  " , IANS 
IF  (IANS.NE.l)  GOTO  20 

ACCEPT" ENTER  GRAYSCALE  FOR  TARGET  (0-15):  " , ITGRAY 

20  ACCEPT "ENTER  STARTING  POSITION  (X,Y):  ",IXST,JYST 

CXPOS= IXST+ IBKOFFX 
CYPGS= I YST+IBKOFFY 

ACCEPT "ENTER  FRAME  RATE:  ", IFRATE 
PIXDI ST=RMPS/ ( IFRATE*TRES) 

RLOWRES=DETS I ZE*TRES 


ESTABLISH  NOISE  PARAMETERS 


IDSX=IFIX ( RDRI FT/ IFRATE + . 5 ) 

IDSX= IFIX ( IDSX/TRES+ . 5 ) 

RDRIFT=IDSX*TRES*FLOAT (IFRATE) 

TYPE"DRIFT  RATE:  " ,RDRIFT, "METERS/SEC" 

TYPE "PIXELS/ FRA:  ",IDSX 

TIE  UP  LOOSE  ENDS 

ACCEPT "ENTER  NUMBER  OF  FRAMES  IN  THIS  SEQUENCE:  " ,INUM 
ACCEPT "PAUSE  AFTER  EACH  FRAME?  ’\IPAUSE 


OPEN  AND  SETUP  OUTPUT  DATA  FILE 


ACCEPT "DO  YOU  WANT  TO  GENERATE  A  SCENE  DATA  FILE?  ",IPR 
IF  (IPR.NE.l)  GOTO  30 

TYPE "ENTER  OUTPUT  FILENAME:  " 

READ (11,13)  OUTFILE ( 1 ) 

OPEN  3, OUTFILE 
NUMDET=NDET 

ILXST=IXST/DETSIZE+IXL/ (2*DETSIZE) +1 
ILYST=IYST/DETSIZE+IYL/ (2*DETSIZE) +1 

TYPE "FOCAL  PLANE  TARGET  COORD:  X:",ILXST,"  Y:",ILYST 

WRITE (3)  INUM  ?  NUMBER  OF  FRAMES 

WRITE (3)  NUMDET  ;  SIZE  OF  FOCAL  PLANE 

WRITE (3)  ISIZE  ;  SIZE  OF  SUB  PIXELS 

WRITE (3)  IFRATE  ;  FRAME  RATE 

WRITE (3)  SIG  ;  STD  DEV  OF  OPTICS  BLUR 

WRITE (3)  RDRIFT  ;  DRIFT  RATE 

WRITE (3)  RLOWRES  ;  FOOTPRINT 

WRITE (3)  ILXST  ;  INITIAL  TARGET  POSITION 

WRITE (3)  ILYST 
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OPEN  2 , TARFILE , ERR=1 01 
GOTO  7 

TYPE"CAN 1 T  OPEN  TARGET  FILE,  ERROR:", IE 
STOP 

ACCEPT " S I MULATE  SENSOR  DRIFT?  " ,IANS 
IF  (IANS.NE.l)  GOTO  8 

ACCEPT " ENTER  DRIFT  RATE  IN  METERS/SEC:  " , RDR1 FT 

ACCEPT "SIMULATE  OPTICAL  TRANSFER  FUNCTION?  " ,NP 
IF  (NPSF.NE.l)  GOTO  9 

ACCEPT "ENTER  SIGMA  OF  OPTICAL  BLUR:  ",S 
RMEAN=DETSIZE/ 2+1 

SET  UP  MEMORY  BUFFERS  FOR  BACKGROUND  VIDEO  FILE 

TYPE "BUILDING  BUFFERS..." 

CALL  PICFMT ( IBUF , BKGFILE , IHDR) 

CALL  MAKB(IBUF) 

CALL  PICIN (IBUF, BKGFILE, IHDR) 

READ  TARGET  DATA  FILE 

IF  ( ITAR.NE . 1 )  GOTO  25 
TYPE 

TYPE "READING  TARGET  DATA  FILE..." 

READ (2, 98)  NPTS 
READ (2,99)  FDIR 
READ (2,99)  TRES 
READ (2, 98)  IXL 
READ (2, 98)  IYL 
DO  10  1=1, NPTS 

READ (2, 98)  TARGETX ( I ) 

READ  (2, 98)  TARGETY ( I ) 

READ (2, 98)  TARGET I ( I ) 

CONTINUE 
FORMAT (15) 

FORMAT (FI 2. 6) 

CLOSE  2 

TYPE 

WRITE (10,11)  (TARFILE(I) ,1=1,3) 

FORMAT (IX,"  DATA  FOR  TARGET  FILE:  ",3A2) 

TYPE"NUMBER  OF  DATA  POINTS:  ",NPTS 

TYPE "TARGET  FACING:  " ,FDIR, "  DEGREES" 

TYPE "TARGET  RESOLUTION:  " , TRES , "  METERS  " 

TYPE "TARGET  SIZE:  ",IXL,"  BY  ",IYL 

TYPE 

GET  TARGET  SPECIFICATIONS 

ACCEPT "ENTER  TARGET  VELOCITY  (MPS):  ",RMPS 
ACCEPT " ENTER  TARGET  DIRECTION  (DEG):  " , RDI R 


PARAMETER  LORGX=155 ,LORGY=50 
PARAMETER  TX=600 ,TY=600 ,TI=600 
PARAMETER  IBKOFFX=30 , IBKOFFY=30 

PARAMETER  NDET=20 
PARAMETER  DETSIZE=9 
PARAMETER  WINSIZE=180 


DET  =  5 ,  BKG  =  78 
DET  =  9 ,  BKG=30 
DET  =  15 ,  BKG=  45 
DET=20 ,  BKG=28 


COMPILER  FREE 

DECLARE  VARIABLES 
INTEGER  XIN , YIN 

DECLARE  ARRAYS 

DIMENSION  BKGFILE (20) ,TARFILE (20) ,OUTFILE (20) 
INTEGER  IBUF (NBUFSZ ) 

INTEGER  TARGETX (TX) , TARGETY (TY) , TARGET I (TI) 
INTEGER  IDETBUF  (WINSIZE ,DETSIZE) 

REAL  COEFF (DETSI ZE , DETS I ZE) ,  RDETROW (WINSI ZE) 

INITIALIZE  VARIABLES  AND  GET  DATA  FILES 

NUMDET=NDET 

ISIZE=DETSIZE 

IBKGWIN=WINSI ZE 

RDRIFT=0 

IDSX=0 

INSX=0 

SIG=0 


VERSION  2.5" 
FOR  ECLIPSE" 


TYPE 

TYPE"  SCENE  GENERATOR  VERSION 

TYPE"  FOR  ECLI 

TYPE 

TYPE"ENTER  BACKGROUND  VIDEO  FILENAME  -> 
READ (11,13)  BKGFILE ( 1 ) 

FORMAT (S20) 

ACCEPT" SIMULATE  TARGET?  "  ,  ITAR 
IF  (ITAR.EQ. 1)  GOTO  6 
RDIR=0 
RMPS=0 
IFRATE=0 
I  XL=0 
I  YL=  0 
IXST=0 
IYST=0 
RLOWRES  =  0 
GOTO  7 

TYPE"ENTER  TARGET  DATA  FILENAME  ->  " 
READ (11,13)  TARFILE(l) 


I ROW= I +BK0FFY-1 

ICOL=BKOFFX+ INSX 

WRITE  (12,180)  (IDONE(L)  ,L=1,20) 

CALL  PROW (I BUF, I ROW, 1 COL, WINS I ZK, I DONE) 
CALL  GROW ( IBUF , I ROW , ICOL , WINS I ZE , I DONE ) 
WRITE (12,180)  (IDONE(L) ,L=1,20) 

CONTINUE 
RX=FLOAT (MX) 


NORMALIZE  PICTURE  TO  4  BIT  DISPLAY  (0-15) 

TYPE"NORMALI Z ING  FOV  ..." 

DO  400  J=1 , WINSIZE 

IROW= J+BKOFFY-1 
ICOL=BKOFFX+ INSX 

CALL  GROW ( I BUF , I ROW, ICOL,WINSIZE , I DONE ) 
RATIO=FLOAT (15) /RX 
DO  350  1  =  1 , WINS  I ZE 

I DONE ( I ) = I F I X ( FLOAT ( I DONE ( I ) ) *RATIO+.75) 
IF  ( I DONE ( I ) . GT .15)  IDONE(I)=15 

CONTINUE 

CALL  PROW (IBUF, IROW, ICOL, WINSIZE, IDONE) 

CONTINUE 

RETURN 


PROGRAM  SCENE 
AUGUST  15,  1984 


VERSION  2.5 
FOR  THE  ECLIPSE 


FUNCTION:  GENERATES  SIMULATED  FOCAL  PLANE  OUTPUT  ( 

TO  TEST  RESOLUTION  ENHANCEMENT  ROUTINES.  ( 

( 

INPUT:  1)  BACKGROUND  VIDEO  FILE  ( 

2)  TARGET  DESCRIPTION  FILE  (FROM  TGEN.FR)  ( 

3)  DESIRED  NOISE  STATISTICS  ( 

4)  DESIRED  TARGET  CHARACTERISTICS  ( 

(SPEED, DIRECTION, .. )  ( 

C 

OUTPUT:  SEQUENCE  OF  FOCAL  PLANE  FRAMES  CONTAINING  ( 

TARGET  AND  BACKGROUND  INFORMATION  ( 

( 

**************************************( 


DECLARE  CONSTANTS 


PARAMETER  NBUFSZ=300 
PARAMETER  HORGX=45,  HORGY=50 


n  n 


40 


CONTINUE 


TYPE"WIDTH  OF  PSF :  ", PNTSPD 
C  GET  FIRST  N-l  ROWS  TO  SET-UP  FOR  CONVOLUTION 

IBND=PNTSPD-1 
DO  50  1=1 , IBND 

IROW=BKOFFY- (PNTSPD/ 2 ) +1-1 
ICOL=BKOFFX+INSX- (PNTSPD/2) 

IWI DTH=WINS I ZE+PNTSPD-1 

CALL  GROW ( IBUF , I ROW, I COL, IWIDTH , I ARRAY (1,1) ) 

50  CONTINUE 

C  SET-UP  ARRAY  POINTER  FOR  PNTSPD  WIDTH  CONVOLUTION 

DO  100  1=1, PNTSPD 

IPTR(I) =1-1 

100  CONTINUE 


C 

0)  i 

150 

C 


COMPUTE  CONVOLUTION 
MX  =  0 

DO  300  1=1 ,WINSIZE 

COMPUTE  NEW  POINTERS 

DO  150  K=l, PNTSPD 

IPTR(K) =MOD ( IPTR (K) , PNTSPD) +1 

CONTINUE 
GET  ROW  N 


175 

C80 

200 


I ROW  = I + (PNTSPD/2) +BKOFFY-1 

CALL  GROW ( I BUF, IROW , ICOL, IWIDTH, I ARRAY (1, IPTR (PNTSPD)  )  ) 
TYPE "STARTING  ROW:  ",I 
DO  200  J=1,WINSIZE 
RVAL=0 

DO  175  K=l, PNTSPD 

DO  175  L=l, PNTSPD 
IX=K+ J-l 
IY=IPTR (L) 

RVAL=RVAL+COEFF ( K , L) *FLOAT 


1 ( I ARRAY ( IX , IY) ) 


IF  (RVAL.GT. 32000)  RVAL=32000 

CONTINUE 

I DONE (J) =IFIX (RVAL+.5) 

FORMAT ( IX , 2  51 5 ) 

IF  (IDONE(J) .GT.MX)  MX=IDONE(J) 


CONTINUE 


C  SAVE  NEWLY  COMPUTED  ROW 
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WRITE (3 ,9756)  RESX 
WRITE (3 ,9756)  RESY 
WRITE (3,9757)  IXL 
DO  9750  1=1, IXL 

INDEX=IXL-I+1 

WRITE  (3,9756)  XAVG(INDEX) 

9750  CONTINUE 

WRITE (3 ,9757)  IYL 
DO  9755  1=1, IYL 

WRITE (3,9756)  YAVG  ( I ) 

9755  CONTINUE 

9756  FORMAT (IX, F12. 6) 

9757  FORMAT (IX, 15) 

CLOSE  3 
9999  END 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  POINT-SPREAD 


FUNCTION:  CONVOLVES  ANY  2-D  PICTURE  INDEXED  BY  PICBUF  ARRAY 
WITH  A  COEFFICIENT  ARRAY  TO  SIMULATE  A  OPTICAL 
POINT-SPREAD  FUNCTION.  (CAN  BE  USED  FOR  ANY  SIZE 
2-D  CONVOLUTION) . 


VARIABLES:  DETSIZE 
BKOFFX 
BKOFFY 
INSX 
WINSIZE 
IBUF 
COEFF 
MX 


SIZE  OF  CONVOLUTION  MATRIX 

OFFSET  OF  SCENE  VIEWING  WINDOW  IN  X  DIR 

OFFSET  OF  SCENE  ”  "  "  Y  DIR 

ADDITIONAL  OFFSET  IN  X  DIRECTION  (DRIFT) 

SIZE  OF  VIEWING  WINDOW  (WINSIZE, WINSIZE) 

UTILITY  ACCESS  FOR  SCENE  (USED  BY  PICBUF) 

CONVOLUTION  COEFFICIENTS 

MAXIMUM  POINT  COMPUTED  (FOR  AVERAGING) 


SUBROUTINE  PSPREAD ( DETS I ZE , BKOFFX , BKOFFY , INSX , WINS I ZE , IBUF , COEFF ) 


INTEGER  IAPRAY (200,9) ,  IDONE(190) 

INTEGER  I PTR  (9)  , IBUF (300) 

REAL  COEFF (DETSIZE, DETSIZE) 

INTEGER  PNTSPD ,  DETSIZE,  BKOFFX, BKOFFY, WINSIZE 

C  INSURE  POINT-SPREAD  RANGE  IS  ODD 

IF  ( (DETSIZE/2-FLOAT (DETSIZE) /2) .NE.0)  GOTO  20 
PNTSPD=DETSI ZE-1 
GOTO  40 

20  PNTSPD=DETSI ZE 
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C  PREPROCESS  X-VECTORS 


MW I N  =  M AXW I N 
ITEST=1 

IF  (IPR.EQ.l)  WRITE (12,9000) 

CALL  VPROC (MWIN , IWPSX ,NSETX , MAXFR ,NBND ,NAVG , XWNOISE , RXSEQ/ XP , I  PR , XAVG , I 
CALL  DECON (XAVG , RREAL , IWINCX , IXL , ITEST , IPR , BLUR) 

I XL= I  XL- IWINCX  +  1 
DO  9555  1=1, MAXFR 

XAVG ( I ) = ABS (RREAL  ( I )  ) 

RREAL (I) =0 

9555  CONTINUE 

IF  (IWPRT.NE.2)  GOTO  9556 

CALL  PLOTIO ( XAVG , INUM ,1,X0,Y0, RESX ) 

9556  IF  (IPR.NE.l)  GOTO  9963 

WRITE (12,9005) 

WRITE (12,9562) 

9562  FORMAT (IX,"  X  TARGET  VECTOR") 

WRITE (12,9020)  (XAVG (I ), 1=1 , IXL) 

WRITE (12,9005) 

WRITE (12,9005) 

C  PREPROCESS  Y-VECTORS 

9963  IF  (IPR.EQ.l)  WRITE ( 12 , 9 0 15 ) 

CALL  VPROC (MWIN , IWPSY , NSETY , MAXFR ,NBND ,NAVG , YWNOI SE , RYSEQ , 

*YP , IPR , YAVG , IYL) 

CALL  DECON (YAVG , RREAL , I WINCY , I YL , ITEST , I  PR , BLUR) 

I YL= IYL-IWINCY+1 
DO  9565  1=1, MAXFR 
YAVG ( I ) = ABS ( RREAL ( I ) ) 

9565  CONTINUE 

IF  (IWPRT.NE.2)  GOTO  9566 

CALL  PLOTIO (YAVG, INUM, 10, X0,Y0,RESY) 

9566  IF  (IPR.NE.l)  GOTO  9980 

WRITE (12,9005) 

WRITE (12,9567) 

9567  FORMAT (IX,"  Y  TARGET  VECTOR") 

WRITE (12,9020)  (YAVG ( I ) , 1=1 , IYL) 

9580  ACCEPT "SAVE  TARGET  VECTORS?  ",ISVECT 

IF  (ISVECT.NE.l)  GOTO  9999 
OPEN  3, "COMP. DA" 


ACCEPT "ENTER  TARGET  DIRECTION:  " 
READ (11,9966)  IDIR 
9966  FORMAT ( A2 ) 

WRITE (3,9968)  IDIR 
9968  FORMAT ( A2 ) 

WRITE (3,9756)  RMPS 


c 

C70 

C 

C80 

C 

C 

C 

C 

750 


IF 


CALL  GR 
( I AXIS . NE . 1 ) 
CALL  GR 
CALL  GR 


PH 2 ( INFILE , NG , XVECT , YVECT , I XNUM , MD , YN ,  YX ,  IFS ) 
GOTO  780 

PH 2 ( INFILE , NG , CONVY , CONVX , I YLEN , MD, YN , YX , IFS) 
PH 2 ( INFILE , NG , CONVX , CONVY , I XLEN , MD ,  YN , YX , IFS ) 


GENERATE  DATABASE 

ACCEPT'1  DO  YOU  WANT  TO  GENERATE  A  DATA  BASE  FILE?  "  ,IANS 
IF  (IANS.NE.l)  GOTO  999 

INFILE (1) = "MD" 


405 


410 


420 


430 


TYPE 

TYPE" 

TYPE 

TYPE" 

TYPE" 

TYPE" 

TYPE" 

TYPE" 

TYPE 

TYPE" 

TYPE" 

TYPE 

ACCEPT" 


CURRENT  PLANE  INVENTORY:  " 

1- 6)  B  -  5  2  " 

2- 7)  KC-135" 

3- 8)  DC-10" 

4- 9)  F-15" 

5  -  10)  MIG-21" 

1ST  NUMBER  -  MODEL  DATABASE  FILE" 
2ND  NUMBER  -  COMPARE  TEST  FILE" 

WHICH  ONE?  " , I P I C K 


TYPE "ENTER  DIRECTION:  (MUST  BE  ’2'  DIGIT  INTEGER  -  00,  09,  45,  ETC) 
READ  (11,405)  INFILE ( 2 ) 

FORMAT ( A2 ) 

ACCEPT "PROCESS  IN  FORTRAN  FIVE  FORMAT?  ",IFIVE 

GO  TO  (410,420,430,440,450,460,470,480,490,495) ,IPICK 


INFILE (3) ="B5" 
INFILE (4) ="2." 
INFILE (5) = " DB " 
IF  ( 1  FIVE . EQ. 1 ) 
INFILE (6) =0 
GOTO  800 
INFILE (3) ="KC" 
INFILE (4) ="13" 
INFILE (5) = " 5 . " 
INFILE (6) = " DB" 
IF  (IFIVE.EQ.l) 
INFILE (7) =0 
GOTO  800 
INFILE (3) = " DC " 
INFILE (4) ="10" 
INFILE (5) ="X. " 
INFILE (6) ="DB" 
IF  (IFIVE.EQ.l) 


INFILE (5) = "DF" 


INFILE (6) = " DF " 


INFILE (6) ="DF" 
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INFILE  (7) =0 
GOTO  800 

440  INFILE (3) ="F1" 

INFILE ( 4 ) =  "5 .  ” 

INFILE (5) ="DB" 

IF  ( I  FIVE  .  EQ .  1 )  INFILE  (5)  =,,DF" 
INFILE (6) =0 
GOTO  800 

450  INFILE (3) ="MI" 

INFILE (4) = " G  2 " 

INFILE (5) ="1.  " 

INFILE (6) ="DB" 

IF  (IFIVE.EQ.l)  INFILE (6) ="DF" 
INFILE (7) =0 
GOTO  800 
460  GOTO  410 

470  GOTO  420 

480  GOTO  430 

490  GOTO  440 

495  GOTO  450 


800  TYPE"OPENING  TARGET  FILE _ " 

IF  ( I  PICK . GT . 5 )  TYPE"ENTER  FILENAME:  " 

IF  (IPICK.GT.5)  READ (11,3)  INFILE(l) 
tm  CALL  OPEN (2, INFILE, 3, IER) 

IF  (IER.NE.l)  STOP"CANT  OPEN  FILE...  SUCK  IT  UP!" 

IF  (IFIVE.EQ.l)  GOTO  850 

IF  (IPICK.LT. 6)  GOTO  602 

TYPE"ENTER  DIRECTION:  " 

READ (11,601)  IDIR 

601  FORMAT (A2) 

WRITE (2,603)  IDIR 
603  FORMAT (IX, A2) 

ACCEPT" ENTER  SPEED:  ",MPS 
WRITE (2)  MPS 
WRITE (2)  RESOLU 

602  WRITE (2) RESOLU 
WRITE (2)  IXNUM 
DO  650  1=1, IXL 

IF  (XVECT(I) .GT.0)  WRITE(2)  XVECT ( I ) 

650  CONTINUE 

WRITE (2)  IYNUM 
DO  700  J=1 , IYL 

IF  (YVECT (J) .GT.0)  WRITE(2)  YVECT(J) 

700  CONTINUE 

CALL  CLOSE  (2) 
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ACCEPT "ANOTHER  TARGET?  ?  ?:  " , IANS 
IF  (IANS.EQ.l)  GOTO  1 
GOTO  999 

850  WRITE (2,890)  RESOLU 

WRITE (2 , 895 )  IXNUM 
DO  860  1=1, IXL 

IF  (XVECT(I) .GT.O)  WRITE(2,895)  XVECT(I) 

860  CONTINUE 

WRITE (2,895)  IYNUM 
DO  870  J=1 , IYL 

IF  (YVECT (J) .GT.O)  WRITE(2,895)  YVECT(J) 

870  CONTINUE 

890  FORMAT (IX, F12. 6) 

895  FORMAT (IX, 15) 

CALL  CLOSE  (2) 

999  TYPE" ALL  DONE  ..." 

END 

c******************************************************************** 

c 

C  PROGRAM  TGEN  -  TARGET  GENERATOR  VERSION  2.0 

C  22  JULY  1984 

C  VERSION  1.0 

C  17  JULY  1984 


INTEGER  WIN,  XPOS , YPOS , IXL , IYL , SCALE , XIN , YIN , PIX ,NUMPOINTS ,CNT 
INTEGER  XSTART ,  XEND,  ASTART , AEND , TEST 
REAL  DIR , VERTLN , RESPIX , RTEMP 

INTEGER  ICT  (120)  ,DISP  (240)  , TARGETX ( 1000 )  , TARGETY ( 1000 )  , TARGET I i 
INTEGER  APIX  (  500 ) 

DIMENSION  INFILE (20) 


(1000) 


INITIALIZE  OCKTEK  DISPLAY 


IFIELD=0 
IX  =  0 
IXL=320 
I  Y=C 
IYL=240 
I ENWV =15 
IENWC= 1 5 
ICMAX= 1 1 
CANGLE=0 
GRID=0 
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CALL  SINTRO  ( ICT , 6 3K , IER) 

IF  (IER.EQ.l)  GOTO  51 

TYPE  "INTRO  ERROR  CODE:  ",IER 
STOP  "UNABLE  TO  INITIALIZE" 

CALL  OPEN ( 3 , " I ACMON . XB " , 2 , 1 ER ) 

IF  (IER.NE.l)  TYPE  "WARNING:  UNABLE  TO  ACCESS  IACMON.XB" 
IF  (IER.EQ.l)  CALL  LXB  (ICT, 3) 

CALL  MPRUN  (ICT, 0,1) 


OPEN  TARGET  DATA  FILE 

TYPE 

TYPE 

TYPE  "WELCOME  TO  TARGET  GENERATOR" 

TYPE 

TYPE  "ENTER  TARGET  FILENAME  (NAME_DI R . DA) :  " 

READ (11,5)  INFILE ( 1 ) 

FORMAT (S 20) 

CALL  OPEN (2, INFILE, 3, IER) 

IF  (IER.NE.l)  TYPE "CAN ' T  OPEN  TARGET  FILE" 

IF  (IER.NE.l)  STOP"ERROR  CODE:  ",IER 

INITIALIZE  VARIABLES 

ACCEPT" IS  THIS  A  TEST:  " , TEST 

WIN=240 

IXL=WIN 

XSTART=32 

XEND=WIN+XSTART-1 

I  YL=  1 

GET  INITIAL  DATA 
TYPE 

ACCEPT "ENTER  DIRECTION  TARGET  IS  FACING  (0-360):  ",DIR 

DIR=DIR*0. 017453 

TYPE 

ACCEPT "ENTER  VERTICAL  SCALE  LENGTH  OF  DISPLAY  (METERS):  " , VERTLN 
RESPIX=VERTLN/ 240 

TYPE "FOR  CURRENT  DISPLAY,  EACH  PIXEL  REPRESENTS ", RESPIX , "MTRS" 
TYPE 

ACCEPT" DO  YOU  WANT  TO  (1) SPECIFY  RESOLUTION  OR  (2) SCALE  ?",IANS 

IF  (IANS.EQ.2)  GOTO  50 

TYPE 

ACCEPT" ENTER  DESIRED  PIXEL  RESOLUTION:  " , RDESRES 
RTEMP=RDES RES /RESPIX 
SCALE=IFIX (RTEMP) 

IDESRES=IFIX (SC ALE* RESPIX) 

TYPE  SCALE, IDESRES 


non 


GOTO  60 


6  • 


TARGET  SCALING  LOOP 

50  ACCEPT "ENTER  SCALE  REDUCTION  (INTEGER):  ", SCALE 

TYPE 

60  IXL=IXL/ SCALE 

TYPE "I XL:  ",IXL 
DO  300  J  =  1 , WIN , SCALE 

TYPE "LINE:  ",J 

DO  200  I =XSTART , XEND , SCALE 

C  CALCULATE  X,Y  COORDINATES  OF  REDUCED  PIXEL 

RTEMP= (I-XSTART) /SCALE+1 
XPOS=IFIX (RTEMP) +1 
RTEMP=J/ SCALE 
YPOS= IFIX (RTEMP) +1 
PIX  =  0 

XIN=SCALE 

IF  (  (XIN+I)  .GT.XEND)  XIN=(XEND-I) 

YIN=SCALE 

IF  ( (YIN+J) .GT.YEND)  YIN= (YEND-J) 

IF  (TEST.EQ. 1)  TYPE "READING  IN  VIDEO  BLOCK:  " 

IF  (TEST.EQ. 1)  TYPE"X  :  "  ,  I , "  Y:",J 
IF  (TEST.EQ. 1)  TYPE "XL : " , XIN , "  YL:",YIN 
CALL  RVBLK ( I CT, API X , I, XIN, J, YIN) 

C  ADD  SCALE  BY  SCALE  PIXELS  INTO  ONE  VALUE,  PIX 

C  FIX  SUMMING  INDICES  DUE  TO  FAULT  IN  OCTEK  ROUTINES 

I SQ= SCALE* SCALE 
ASTART=2 
AEND= I SQ+ 1 

IF  (TEST.EQ. 1)  TYPE" SUMMING  ARRAY:  " , ASTART , AEND 
DO  100  K= AST ART , AEND 

PIX=PIX+APIX (K) 

100  CONTINUE 

C  SAVE  PIXEL  VALUE  IN  ARRAY 

RTEMP=PIX/ (SCALE*SCALE) 

IF  (TEST.EQ. 1 )  TYPE" PIX ( " , I , J , " =  ", RTEMP 
DISP (XPOS) =IFIX (RTEMP) 

200  CONTINUE 

C  WRITE  ARRAY  TO  DISPLAY  (ONE  VIDEO  LINE) 

CALL  WVBLK  ( ICT , DISP , XSTART , IXL,YPOS, IYL) 

300  CONTINUE 
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ono  jk  o  ono 


RTEMP=WIN/ SCALE 
WIN=IFIX (RTEMP) 

XEND=XSTART+WIN-1 

RESPIX=VERTLN/IXL 

TYPE 

TYPE" PRESENT  SCALE:  1  PIXEL= " , RESPIX , "METERS" 
ACCEPT "ANOTHER  REDUCTION?  " ,  IANS 
IF  (IANS.EQ.l)  GOTO  50 

READ  FINAL  TARGET  INTO  ARRAY 

TYPE 

TYPE "READING  TARGET  DISPLAY" 

CNT  =  1 
IY =WIN-1 
DO  400  J= 1 , IY 

DO  400  I =XSTART , XEND 
READ  PIXEL  I,J 

IVAL= IRDPI X ( ICT / I , J ) 

IF  ( IVAL.GT. 14 )  GOTO  400 

TARGETX (CNT ) =I-XSTART+1 
TARGETY (CNT) =J 
TARGET I (CNT) =IVAL 
CNT=CNT+1 

CONTINUE 

NUMPOINTS=CNT— 1 
IRANGE=XEND-XSTART 

TYPE  "NUMBER  OF  POINTS:  ", NUMPOINTS 
TYPE  "WINDOW:  " , IRANGE ,WIN 

OUTPUT  DATA  TO  DISK 

TYPE 

TYPE "WRITING  DATA  TO  DISK..." 

WRITE (2)  NUMPOINTS 
WRITE (2)  DIR 
WRITE (2)  RESPIX 
WRITE (2)  WIN 
WRITE (2)  WIN 


DO  500  1=1 , NUMPOINTS 

WRITE (2)  TARGETX (I) 
WRITE (2)  TARGETY (I) 
WRITE (2)  TARGET I ( I ) 

CONTINUE 


TYPE "ALL  DONE 


£********************************************************************* 

c 

C  SUBROUTINE  VECTOR-PREPROCESS  LANGUAGE:  FORTRAN  V 

C 

C  CALLED  BY:  PROFILE. FR 

C 

C  FUNCTION:  REDUCES  SAMPLED  WINDOW  DATA  TO  SINGLE  X-Y 

C  VECTORS  BY  SUBTRACTING  DC  NOISE  AND  COMBINING 

C  WINDOWS  BASED  ON  WINDOW  CONTRAST 

C 

c** ***************  *************************************************** 

SUBROUTINE  VPROC (MXW , I WPS ,  NSET , MXFR , NBND ,NAVG , RN , RSEQ , 

1MXPNT , I PR , AVG , ICNT ) 


(9j 


INTEGER  MXW 
INTEGER  MXFR 
INTEGER  NBND 
INTEGER  NAVG 
INTEGER  IPR 
INTEGER  ICNT 
INTEGER  IWPS 
INTEGER  NSET 


; MAX  NUMBER  OF  WINDOWS 
;MAX  NUMBER  OF  FRAMES 
;  NUMBER  OF  WINDOWS 
; NUMBER  OF  WINDOWS  TO  AVERAGE 
; PRINT  SWITCH 

; RETURN  LENGTH  OF  AVG  VECTOR 
.•WINDOWS  PER  SET 
; SETS  OF  WINDOWS 


INTEGER  MXPNT (MXW) 
REAL  RN (MXW) 

REAL  RSEQ (MXW, MXFR) 
REAL  AVG (MXFR) 
INTEGER  I SORT (10) 


; TABLE  OF  WINDOW  START  POSITIONS 
; DC  NOISE  OF  WINDOW 
; WINDOW  DATA 
; RETURN  VECTOR 
; INTERNAL  WINDOW  SORT  ARRAY 


DO  400  1=1, IWPS 


C  INITIALIZE  VARIABLES 


ICNT  =  0 
ILSTBND=0 
ICURBND= I 
RMINNOISE=1000 . 0 
RMAXNOISE=0 
ISHIFT=0 


DO  50  J=-1,MXFR 

AVG (J) =0 

50  CONTINUE 

DO  300  J= 1 , NSET 

C  FIND  NEXT  LOWEST  SNR  WINDOW 

RMINNOISE=1000 
DO  100  K=I, NBND, IWPS 

IF  (RN (K) . LE . RMAXNOI SE)  GOTO  100 
IF  (RN (K) .GT.RMINNOISE)  GOTO  100 
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I 


RMINNOISE=RN(K) 

ICURBND=K 

,  100  CONTINUE 

-  RMAXNOISE=RM INNOISE  ;  LOW  SNR  NOW  THRESHOLD 

C  PRINT  RESULTS 

WRITE (10,110)  I CURBND,RM INNOISE 
;  IF  (IPR.NE.l)  GOTO  120 

■  WRITE (12,110)  ICURBND , RM INNOISE 

110  FORMAT (IX," SELECTED  WIN:  "  ,14,"  WITH  DC 

1  NOISE:  " , F 1 2 . 6 ) 

C  SHIFT  BOUNDRIES  TO  ALIGN  FOR  AVERAGING 

1  120  IF  (J.LT.2)  GOTO  130  ;  DONT  SHIFT  FIRST  CHOICE 

ISHIFT=MXPNT (ICURBND) -ILSTBND+ISHIFT 
130  ILSTBND=MXPNT (ICURBND) 

C  AVERAGE  SELECTED  WINDOWS 

DO  200  K=1,MXFR 

INDEX =K+ I SHI FT 
IF  ( INDEX.LT. 1)  GOTO  200 
IF  (INDEX. GT.MXFR)  GOTO  200 
RTEMP=RSEQ ( ICURBND , INDEX ) 

IF  (RTEMP.GT. AVG (K) )  AVG ( K ) =RTEMP 

CONTINUE 

CONTINUE 

C  PUT  AVERAGED  VECTOR  IN  FIRST  WINDOW  SET 

IRCNT=1 

DO  325  J=1 , MXFR 

IF  (AVG(J) .LE.0)  GOTO  325 

IF  (IRCNT. GT.MXFR)  GOTO  325 
RSEQ ( I , IRCNT) =AVG (J) 

I RCNT= I RCNT+ 1 

325  CONTINUE 

TYPE"NUMBER  OF  POINTS  IN  VECTOR (  MINUS  1):  ", IRCNT 
DO  350  J  =  IRCNT , MXFR 
RSEQ (I , J) =0 

350  CONTINUE 

WRITE (10,620)  ( RSEQ ( I , J) ,J=1,MXFR) 

400  CONTINUE 

C  ORDER  WINDOWS  FOR  RECOMBINATION  BY  CHOOSING  LOWEST  FIRST  POINT 

RMAX=0 
RMIN=1000 
DO  410  I=1,IWPS 

I SORT ( I ) =1 
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410  CONTINUE 

DO  450  I=1,IWPS 

RMIN=1000 
DO  425  J=1,IWPS 

IF  (RSEQ(J,1) .LE.RMAX)  GOTO  425 

IF  (RS£Q(J,1) .GT.RMIN)  GOTO  425 

RMIN=RSEQ(J,1) 

I SORT (I) =J 

425  CONTINUE 

RMAX=RMIN 

450  CONTINUE 

DO  452  I=1,IWPS 

TYPE  1,1  SORT ( I ) 

452  CONTINUE 

C  COMBINE  WINDOW  SET  TO  FORM  CONVOLUTED  VECTOR 

ICNT=0 

RMAX=0 

DO  500  1=1 ,MXFR 

DO  500  J=1,IWPS 

INDEX= ( I — 1 ) *IWPS+J 
IF  (INDEX. GT.MXFR)  GOTO  500 
AVG ( INDEX) =RSEQ( I SORT (J) ,1) 

IF  (AVG (INDEX) .GT.0)  ICNT=ICNT+1 
IF  (AVG (INDEX) .GT.RMAX)  RMAX=AVG ( INDEX ) 

500  CONTINUE 

WRITE (10,620)  ( AVG ( I ) ,1=1,1 CNT ) 

WRITE (12,610) 

WRITE (12,620)  (AVG ( I ) , 1=1 , ICNT) 

DO  600  1=1 ,MXFR 

AVG ( I ) =AVG ( I ) /RMAX 

600  CONTINUE 

IF  (IPR.NE.l)  GOTO  999 
WRITE (12,610) 

WRITE (12,620)  ( AVG ( I )  ,  I  =  1 ,  ICNT ) 

610  FORMAT (IX,"  AVERAGE  CONVOLUTION:") 

620  FORMAT (1X,5F12. 6) 

999  RETURN 

END 
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